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Introduction 

In  this  final  year  of  our  award,  we  have  made  significant  progress  on  several  areas  of 
theoretical  research  on  the  prion  problem,  particularly  on  a  promising  new  model  for  the 
in  vitro  mammalian  prion  fibril  which  involves  a  new  suggestion  for  a  C-terminal  region 
left  handed  beta  helix.  We  have  also  made  progress  on  exploring  other  domain  swapped 
structures  for  oligomers  and  on  the  role  of  mutations  in  our  domain  swapped  prion  trimer 
model.  Overall  on  the  award,  we  have  accomplished  much:  (i)  published  seven  papers 
and  one  PhD  dissertation  (two  more  are  in  the  works  this  year),  (ii)  presented  27  invited 
lectures  at  universities  or  meetings  and  ten  contributed  talks  at  conferences;  (iii)  co¬ 
organized  two  major  international  workshop  on  amyloid  diseases  (the  last  being 
http://i2cam.org/i2camyloid);  (iv)  received  national  press  coverage  for  our  research  on 
prion  diseases  (see  the  attached  item  in  Appendix  J);  (v)  received  fellowship  honors  for 
the  participants  (a  Guggenheim  for  the  PI,  a  Kevli  fellowship  at  the  Kevli  Institute  for 
Theoretical  Physics  of  UC  Santa  Barbara  for  former  student  David  Mobley,  and  two 
fellowships  for  finishing  student  Jianping  Pan). 

We  anticipate  2-3  more  papers  emerging  from  the  final  year  of  research  on  this  grant  and 
we  detail  those  areas  of  progress  in  the  report  that  follows,  focusing  on  the  above 
mentioned  work  on  a  new  model  for  in  vitro  grown  fibrils,  new  domain  swapped 
oligomer  models,  and  models  of  the  influence  of  proline  mutations  on  the  kinetics  of 
formation  of  the  domain  swapped  prion  trimer. 


Key  Research  Accomplishments  During  the  Award  Period 

1.  Modeling  of  Oligomeric  Intermediates  in  prion  aggregation.  See  App.  B,  manuscript 

published  in  Biophysical  Journal. 

2.  Modeling  of  insertion  of  Abeta  peptides  in  model  membranes.  See  App.  C, 
manuscript  published  in  Biophysical  Journal. 

3.  Dissertation  on  modeling  of  prions  and  Abeta.  See  App.  D,  directions  to  online  copy 
of  dissertation  by  David  Mobley,  with  strong  review  chapter. 

4.  Proposing  a  new  structure  for  the  minimal  infectious  unit  of  prion  disease:  the 
Domain  Swapped  Prion  Trimer  (DSTP).  See  App.  E,  manuscript  of  paper  published 
in  FASEB  Journal. 

5.  Demonstration  of  the  necessity  of  the  membrane  in  exponential  in  vivo  growth  of 
prion  infectivity.  See  App.  F,  manuscript  published  in  the  Letters  section  of 
Biophysical  Journal. 

6.  Possible  mechanism  for  copper  inhibition  of  prion  conversion.  See  App.  G, 
manuscript  published  in  the  Letters  section  of  Biophysical  Journal. 
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7.  Detailed  modeling  of  in  vitro  yeast  prion  aggregation  See  App.  H,  manuscript  of 
paper  published  in  Phys.  Rev.  E. 


8.  Review  of  amyloidogenic  proteins  from  a  materials  perspective.  See  App.  I,  published 
in  Materials  Research  Bulletin. 


Research  Accomplishments  during  final  term  of  support 


Modeling  of  in  vitro  fibrils  with  left  handed  beta  helices. 


Figure  1 :  Ribbon  representations  of  the  nearly  four  turn  C-terminal  LH(3H 
model  (CTP4)  and  the  nearly  two  turn  N-terminal  LH(3H  (NTP2)  model  . 
Notable  in  CTP4  is  the  disulfide  bond  (cysteines  highlighted),  and  in  NTP2  the 
removed  hydrophobic  loop  of  dominant  glycine  and  alanine  character. 

Figures  drawn  with  RASTOP. 


Recent  work1  on  in  vitro  grown  mammalian  prion  fibrils  in  low  pH  shows  electron 
microscopy  evidence  for  prion  fibrils  with  twinned  crosslinked  filaments,  possessing 
repeat  units  with  lengths  consistent  with  twelve  beta  strands  perpendicular  to  the  fibrillar 
axis,  no  alpha  helical  structure,  and  no  strong  evidence  for  inter-monomer  disulfide  bonds 
between  cysteines.  The  cross  sectional  diameter  of  the  filaments  is  of  order  25-35 
angstroms,  and  the  linkers  between  the  filaments  are  of  length  50-60  angstroms.  There  is 
a  gap  of  order  8-15  angstroms  between  these  repeat  units  which  is  thus  likely  devoid  of 
beta  sheet  bonding.  The  cross  linked  regions  are  electron  rich  and  approximately  twice 
as  wide  as  the  "empty”  regions  between  filaments.  These  facts  alone  allow  us  to  begin 
to  build  a  powerful  empirical  case  for  LH|]H  structure  to  the  filaments.  We  note  the 
following  observations: 
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•  With  the  inclusion  of  side  chains  pointing  to  the  outside,  it  is  easy  to  see  that  the 
cross  sectional  size  of  LH(3H  motifs  are  of  the  order  of  25-35  angstroms. 

•  The  two  versions  of  LH|iH  trimers  proposed  already  for  the  structures  of  prion 
oligomers2'3  contain  LH(3H  structures  in  each  monomer  on  the  N-terminal  side  of  the 
truncated  prion  protein  with  3  beta  helical  turns  in  residues  90-144. 

•  The  remnant  alpha  helical  portions  of  the  original  LH|iH  trimer  model  has  63 
residues  between  the  structure  disrupting  proline  at  residue  165  and  the  terminal 
glutamine  at  residue  228,  easily  enough  to  accommodate  three  full  LH(3H  turns. 

•  The  two  cysteines  with  a  disulfide  bond  are  at  residues  179  and  214  respectively, 
differing  by  35  residues  or  almost  two  full  turns  of  a  LH|]H.  Indeed,  assuming  the 
cysteines  are  at  LH|]H  corners,  there  is  no  difficulty  in  accommodating  the  intra¬ 
monomer  disulfide  bond.  This  requirement  fixes  the  three  turn  C  LH|iH  threading  to 
within  1-2  residues. 

•  The  21  residue  linking  region  between  residues  145  and  165  is,  fully  extended, 
approximately  60  angstroms  and  thus  of  sufficient  length  to  explain  the  observed  cross 
links. 

This  suggests  that  we  can  build  a  model  in  which  each  filamentary  repeat  is  composed  of 
two  N-terminal  and  two  C-terminal  derived  LH|iH  s  such  that  they  add  to  give  six  turns. 
Fig.  (1)  illustrates  an  example  of  such  LH(3H  s,  in  which  the  C-  LH|iH  has  almost  four 
turns  (CTP4)  and  the  N-  LH(3H  has  just  over  two  turns  (NTP2).  In  order  to  provide 
enhanced  stability  relative  to  the  threading  of  the  prion  trimer  models  mentioned  earlier, 
we  have  pulled  out  a  predominant  G-A  loop  which  contains  the  putative  transmembrane 
segment  of  the  prion  protein(4).  We  have  also  considered  fibril  models  built  from  three 
turns  at  the  C-terminal  (CTP3)  and  from  taking  the  three  turn  LH|iH  of  the  miniprion 
trimer  model  for  the  N-terminal  LH|iH  (NTP3). 

There  is  considerable  latitude  in  building  a  suitable  model.  Figure  (2)  gives  an 
illustration,  which  shows  how  we  can  get  into  the  right  general  region  of  model  space  to 
describe  the  fibrils  by  stacking  in  an  overall  prion  tetramer  the  C  and  N  terminal  LH[]H  s 
(C  LH|iH  and  N  LH(3H)  with  the  arrow  denoting  the  n-to-c  terminal  sense  of  the  helix. 
The  individual  monomers  are  color  coded  in  the  schematic  diagram  and  in  the  ribbon 
diagram.  The  difference  between  the  two  schematics  is  the  sense  of  the  filaments;  on  the 
left,  they  stack  in  the  same  sense,  while  on  the  right  they  stack  in  opposite  sense. 


Figure  2:  Sample  models  of  mammalian 
prion  fibril  repeat  units,  (a)  shows  same 
sense  filaments  with  color  coding  indicating 
the  partner  CLH(3H/NLHpH  structures  from 
a  given  monomer;  (b)  shows  antisense 
filaments;  (c)  is  a  ribbon  rendering  of 
model  (a)  generated  within  VMD  using 
CTP3  and  NTP3  models  for  the  C  and  N 
terminal  LH(3Hs.  Image  of  (c)  created  with 
VMD. 
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A  separate  degree  of  freedom  is  the  number  of  turns  per  helix;  we  have  considered  both  3 
and  4  turn  CLH(3H  (abbreviated  CTP3  [residues  178-226  of  human  PrP],  CTP4  [residues 
166-216  of  human  PrP  ])  and  2  and  3  turn  NLH(3H  (abbreviated  NTP2  [residues  90-145 
of  human  PrP],  NTP3  [residues  90-145  of  human  PrP]).  NPT3  is  the  three  turn  region 
taken  from  the  proposed  LH(3H  trimer  model,  while  NPT2  is  a  new  model  we  have 
produced  in  which  we  extract  a  loop  [residues  114-126  of  human  PrP]  of  predominantly 
glycine  and  alanine  residues  known  to  belong  to  the  section  of  the  prion  which  can  cross 
the  membrane  [residues  111-134  of  human  PrP]4.  We  find  that  the  NPT3  LH(3H  by  itself 
is  very  unstable,  in  large  part  due  to  the  poor  packing  of  the  middle  layer  which  contains 
substantial  glycine  and  alanine  content,  and  due  to  poor  side-chain  to  side-chain 
hydrogen  bonding.  As  discussed  below  the  NPT2  LH(3H  is,  apart  from  the  fluctuating 
loop  region,  as  stable  as  the  two  CLHpH  models,  and  contains  reasonable  side-chain  to 
side-chain  hydrogen  bonding. 


Fig  3:  Embedding  of  model 

tetramer  composed  of  NTP2, 
CTP4  LH(3Hs  plus  connecting 
loops  in  experimental  density  data 
of  Ref.  1  Image  created  with 
Chimera. 


Fig.  (3)  shows  a  model  in  which  we  employed  approximate  mirror  symmetry  to  agree 
with  the  densities  of  the  1200  angstrom  helix  of  Ref.  (1)  If  we  make  the  assumption  that 
the  N-terminal  region  of  the  NLH(3H  mirrors  across  the  gap  of  the  density,  then  we  can 
see  that  the  clumping  of  the  cross  links  gives  rise  to  a  low/high  density  variation  similar 
to  that  seen  in  experiment.  An  attractive  feature  of  this  model  is  that  previously  proposed 
domain  swapping  using  the  proline  containing  loop  (K101-P105)(see  Ref.  3)  across  the 
gap  is  possible  and  can  be  a  stabilizing  feature. 

The  space  of  exploration  is  relatively  large  even  within  this  simplified  model  in  that  we 
can  vary 

•  over  four  relative  sense  orientations  in  each  filament  of  the  tetramer  (in  addition 
to  that  shown,  there  are  CtN|N|C|,  C|N|NtC|,  and  N|CtC |N  j). 

•  The  filaments  can  run  in  antisense  to  one  another  as  shown  in  the  figure. 
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•  In  addition  to  the  long  filament  domain  swapping  associated  with  loops  between 
NLHpHs  and  CLHfiHs,  there  can  be  domain  swapping  between  adjacent  NLH|)Hs  of  the 
type  discussed  in  Ref.(3) 

•  The  orientation  of  the  triangular  cross  section  of  the  beta  helix  is  not  readily 
determined  by  the  available  data  though  it  is  constrained  by  the  aspect  ratio  of  the  density 
cuts. 

•  There  is  flexibility  on  how  many  turns  to  include  in  the  CLH(3H  and  NLH(3H. 

An  important  key  to  assessing  the  model  on  theoretical  grounds  is  the  issue  of  the 
stability  of  the  CLH(3H.  We  show  in  Fig.  (4)  the  RMSDs  from  selected  three  turn 
segments  of  five  LH|3Hs  from  the  PDB  together  with  the  CTP3,  CTP4,  and  NTP2  models 
of  LH(3Hs  from  our  work.  The  short  time  frame  for  the  simulation  is  not  intended  to 
confirm  stability  in  any  real  sense,  but  rather  to  be  comparative  with  known  structures, 
for  which  the  RMSD  ranges  from  1-3  Angstroms.  (The  NTP3  model  drawn  from 
existing  prion  trimer  structure  models  It  is  clear  that  these  preliminary,  unoptimized  prion 
models  are,  in  this  measure,  within  the  range  of  viability  suggested  by  known  stable 
structures.  In  particular,  the  CTP3  model  looks  quite  promising.  We  suspect  this  is 
facilitated  by  (i)  excellent  side  chain-to-side  chain  hydrogen  bonding,  and  (ii)  added 
stability  from  the  covalent  disulfide  bond  between  the  cysteines. 

Our  plan  is  to  examine  those  models  which  have  the  best  overall  chance  of  describing  the 
fibril  data,  and  then  optimize  the  factors  above  to  give  the  best  agreement  with 
experiment  using  construction  programs,  electron  density  map  generation  with  EMAN 
and  electron  density  map  optimization  and  visualization  via  tools  in  the  CHIMERA 
package.  We  will  then  prepare  all  atom  AMBER  MD  runs  for  octamers  of  prion  protein 
to  ensure  that  we  build  in  the  gap  region  and  validate  the  stability  with  relatively  low 
RMSD  from  the  starting  structure  (1-2  angstroms). 


0  100  200  300  400  500  600  700  800  900  1000 

Time  (ps) 


Figure  4:  RMSD  from  starting 
structure  in  all  atom  Amber  MD  runs 
(npt  at  300K)  for  five  LH(3Hs  from  the 
PDB,  and  for  the  3  and  4  turn  C- 
terminal  (CTP3,CTP4)  and  2  turn  N- 
terminal  (NTP2)  LH(3H  models 
discussed  here.  For  the  PDB 
proteins,  starting  structures  were  the 
posted  coordinate  files. 


We  note  that  the  conformational  diversity  obvious  within  these  fibril  models  may  have 
some  relevance  to  the  conformational  diversity  observed  both  in  vitro5  and  in  the  prion 
strains.  We  are  particularly  intrigued  by  the  notable  stability  of  the  CTP3  and  CTP4 
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models  in  the  following  sense:  small  amounts  of  these  structures  could  play  significant 
templating  roles  both  in  vivo  and  in  vitro.  To  satisfy  the  general  constraint  about  the 
amount  of  alpha  helical  content  in  infectious  material,  which  must  derive  from  the  C- 
terminal  region,  we  cannot  have  large  amounts  present  from  purified  brain  extracts.  On 
the  other  hand,  the  relatively  small  amount  of  proteinase  resistant  material  found  in  many 
strains  by  conformation  dependent  immunoassay6  may  point  to  a  critical  role  for  these 
CLH[iHs  in  disease  propagation,  since  it  is  C-terminal  stretches  that  provide  the 
proteinase  resistant  core  in  synthesized  fibrils  '  . 

Stabilization  of  LHj3Hs.  A  fundamental  question  we  need  to  address  is  "what  stabilizes 
the  beta-helices?”  There  is  an  obvious  cost  to  these  structures  as  making  the  beta  helices 
requires  bends  or  turns  in  the  beta  strands,  so  they  can  wrap  on  top  of  each  other  allowing 
the  formation  of  backbone  Hydrogen  bonds.  These  turns  can  be  gradual  as  in  Perutz's 
proposal  of  a  circular  beta-helix9,  or,  they  can  be  sharp,  as  in  the  triangular  beta-helices 
commonly  found  in  the  protein  databases1011.  In  either  case,  these  deviations  from  planar 
structure  cost  torsional  energy  and  can  reduce  the  strength  of  the  backbone  hydrogen 
bonds.  Thus,  one  needs  other  factors  that  can  compensate  for  this  increase  in  energy  and 
help  stabilize  the  beta-helical  structure. 


Figure  5  :  Side  chain-to- 
side  chain  hydrogen  bonds 
per  turn  for  a  number  of 
LH(3Hs  from  the  Protein  Data 
Bank,  together  with  our  C- 
terminus  and  N-terminus 
LHfSH  models 
(CTP4,CTP3,NTP3,NTP2) 


□  [Frustration  lndex| 
■  Packing  Fraction 


Figure  6:  Frustration  index 
(CTP4,CTP3,NTP3,NTP2) 
|Nm|  and  packing  fraction  pf 
for  several  LH(3Hs  from  the 
PDB  along  with  the  three 
turn  C-terminal  LH(3H  model 
discussed  in  this  proposal. 
AVG  is  the  average  of  these 
quantities  for  the  PDB 
proteins. 


Studying  the  beta-helical  structures  in  the  protein  databases,  we  have  found  three 
important  stabilizing  effects: 

1.  Side-chain  Hydrogen  Bonding:  The  beta-helix  structure  allows  side-chain  to  side-chain 
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as  well  as  side-chain  to  backbone  Hydrogen  bonding,  for  suitably  placed  residues.  We 
show  in  Fig.  (5)  that  the  CTP3,  CTP4  models  have  superior  side  chain  bonding  near  the 
best  of  selected  LH(3Hs  from  the  PDB,  while  the  as-threaded  NTP3  used  in  modeling 
prion  oligomers  has  essentially  no  side  chain-to-side  chain  hydrogen  bonding  which  may 
account  for  its  relative  instability.  This  situation  is  substantially  improved  in  the  NTP2 
model  with  the  extruded  loop. 

2.  Good  packing  of  the  hydrophobic  interior:  The  known  beta-helices  pack  the  interior  of 
the  helix  rather  well  and  are  very  compact.  This  can  also  lead  to  a  substantially  lower  free 
energy  for  the  structure.  This  is  illustrated  quantitatively  in  Fig.  (6)  for  seven  LH|iHs 
from  the  PDB  and  the  CTP3  model.  This  shows  that  the  CTP3  model  is  at  the  low  end  of 
packing  fractions  for  known  stable  LH(3Hs.  The  situation  is  visualized  in  Fig.  (7)  which 
shows  a  comparison  of  the  packing  in  three  turns  of  the  beta  helix  labeled  1KGQ  on  the 
PDB,  and  the  CTP3  and  NTP3  models,  where  we  have  placed  one  angstrom  spheres 
around  each  atom.  Clearly  1KGQ  has  better  packing,  but  the  CTP3  and  NTP3  models  are 
not  substantially  worse. 

3.  Good  matching  of  hydrophobic  residues  in  the  interior  and  hydrophilic  residues  in  the 
exterior  of  the  beta-helix:  By  reducing  any  mismatches  of  hydrophobic  and  hydrophilic 
residues  the  beta-helix  structure  can  gain  substantial  stability.  Fig.  (6)  shows  a 
comparison  of  seven  LH(3Hs  from  the  PDB  together  with  the  CTP3  model.  On  this 
account,  the  CTP3  model  is  in  very  reasonable  accord  with  known  stable  structures. 


Figure  7:  Top  view  of  internal  packing  of 
three  turn  segments  of  LH(3Hs.  Three  turns 
have  been  taken  from  1KGQ  on  the  PDB,  and 
the  three  turn  CLH(3H  and  NLH(3H  models  for 
the  prion  protein  are  shown.  In  each  case,  one 
angstrom  spheres  are  centered  on  each  atom. 
Drawings  rendered  by  RASTOP. 


These  observations  are  less  pertinent  to  a  homo-peptide  such  as  poly-glutamine,  but  for  a 
hetero-peptide  such  as  part  of  the  prion  protein  or  A-beta  peptide,  this  suggests  that  any 
proposal  for  a  threading  of  a  sequence  onto  a  beta-helical  structure  should  take  these 
factors  into  consideration.  Here,  we  propose  to  develop  a  scheme  for  finding  stable  beta- 
helical  structures,  for  different  sequences,  guided  by  the  known  structures  in  the 
databases. 

In  threading  a  peptide  sequence  onto  a  beta-helical  structure,  there  is  some  flexibility  as 
to  which  peptide  is  placed  in  which  location.  This  flexibility  is  further  increased,  when 
loops  are  set  aside  that  can  take  on  a  non-beta  structure.  We  plan  to  use  the  following 
strategy  to  obtain  stable  beta-helix  structures  given  a  peptide  sequence.  First,  we  will  use 
standard  programs  to  thread  the  peptide  sequence  onto  the  beta-helical  structure, 
including  selecting  the  loop  regions  which  will  have  non  beta-strand  dihedral  angles.  For 
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this  trial  structure,  after  some  energy  minimization  within,  e.g.,  the  Amber  package,  we 
will  measure  the  number  of  side-chain  Hydrogen  bonds  Ns,  the  packing  fraction  pf 
(defined  as  the  ratio  of  the  enclosed  side  chain  volume  to  the  volume  of  the  subtended 
triangular  prism),  and  the  number  of  satisfied  Hydrophilic-Hydrophobic  matches  Nm  ,  in 
the  beta-helical  region.  (Nm  is  defined  as  the  number  of  polar/charged  residues  pointing 
out  of  the  helix  minus  those  pointing  into  the  helix  plus  the  number  of  hydrophobic 
residues  pointing  into  the  helix  minus  the  number  pointing  out  of  the  helix.)  We  define  a 
dimensionless  scoring  function  p,  with  <NS>,  e.g.,  the  mean  of  Ns  from  LH|)Hs  of  the 
PDB, 


p=  -(Ns/<Ns>  +  pf/<pm>+Nm/<Nm>)/3 

By  normalizing  each  of  the  three  variables  to  their  mean  values,  we  equally  weight  their 
importance.  Obviously  we  can  explore  other  weighting  procedures  as  we  consider 
simulations  on  LH|)H  structures.  The  lower  the  score,  the  more  stable  the  structure.  Note 
that  the  sign  is  positive  for  Nm  in  the  score  as  lower  Nm  will  be  more  stable.  The  mean  p 
score  for  the  structures  in  the  database  is  -1,  and  its  variance  is,  assuming  no  correlations 
between  Ns,  Nm,  and  pf,  and  with,  e.g.,  ANm  the  variance  of  Nm 

Ap2  =  (ANs2/<Ns>2  +Apf2/<pf>2+  ANm2/<Nm>2)/9. 

If  we  find  that  the  score  for  a  trial  threaded  structure  is  within  Ap  of  -1  ,  it  is  a  good 
candidate  beta-helical  structure,  and  will  be  explored  more  closely  through  a  full 
Molecular  Dynamics  simulation.  If  it  has  a  higher  score,  it  is  unlikely  to  be  stable  and 
will  be  rejected.  We  plan  to  use  this  scheme  to  obtain  stable  beta-helical  threading  for 
various  known  amyloid  proteins.  We  will  test  the  scheme  for  (i)  independence  of  the 
main  scoring  variables,  (ii)  correlations  between  particular  residue  frequencies  in  the 
PDB  and  good  values  of  p,  and  (iii)  intra-molecular  correlations  among  residues. 

With  the  specific  examination  of  the  NLH|)H.  if  we  are  unable  to  stabilize  via 
rethreading,  we  will  explore  the  possibility  that  extrinsic  perturbations,  especially  in  the 
form  of  a  membrane,  might  provide  the  necessary  stability.  In  this  case  we  will  attempt 
first  to  model  the  membrane  as  a  mean  field  with  mildly  negative  surface  charge  drawing 
upon  our  experience  in  modeling  A(340-42  insertion  in  membranes12.  In  particular,  we 
will  investigate  how  well  the  G-A  loop  in  the  NTP2  model  can  insert  into  the  membrane, 
and  will  update  the  simulation  model,  which  is  biased  towards  alpha  helices,  to  allow  for 
formation  of  beta  structure  outside. 

Disposition  to  prion  disease  by  point  mutations.  With  our  domain  swapped  prion  trimer 
model,  the  short  loop  K1 01 -PI 05  has  special  significance  in  terms  of  gluing  the  oligomer 
together.  This  is  also  potentially  relevant  to  fibril  growth  as  discussed  earlier  in  this 
section.  Two  of  the  point  mutations  which  dispose  individuals  to  get  the  inherited  GSS 
form  of  prion  disease  involve  the  mutation  of  prolines  for  leucines  on  this  loop  (P102L 
and  P105L),  and  one  reported  inherited  form  of  Creutzfeldt-Jakob  Disease  involves  the 
mutation  P105T  .  We  note  that  these  mutations  can  promote  enhanced  conversion  from 
cellular  prion  to  scrapies  form  by  at  least  two  effects:  (i)  any  substitution  for  proline  in 
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these  loops  will  make  them  more  flexible  and  hence  more  likely  to  convert  rapidly  if  this 
conversion  is  important  to  stabilizing  the  aggregate,  which  is  true  of  both  the  domain 
swapped  oligomers  and  fibrils  and  (ii)  these  substitutions  as  well  as  artificial  ones  studied 
in  transgenic  mice  by  Nazor  et  alu  are  relatively  hydrophobic  to  proline,  which  means 
they  are  inherently  destabilizing  to  the  nominal  random  coil  of  the  wild  type  cellular 
form,  while  if  there  is  favorable  water  exclusion  by  hydrophobic  clustering  in  the  core  of 
the  domain  swapped  oligomers  then  the  energy  surface  from  cellular  prion  form  to 
aggregate  will  be  ramped  further.  We  have  made  a  preliminary  examination  of  the 
influence  of  mutants  on  the  loop  flexibility  as  measured  by  the  associated  B 
factors, defined  for  a  given  atom  as  B  =  87t<u  >/3. 

Fig.  (8)  illustrates  schematically  within  a  simple  two-state  model  how  the  simultaneous 
hydrophobic  destabilization  of  the  wild  type  protein  and  hydrophobic  stabilization  of  the 
domain  swapped  trimer  can  ramp  the  energy  surface  to  lower  the  energy  barrier  and 
speed  the  conversion  reaction.  Obviously  this  naive  picture  ignores  intermediates,  but  we 
anticipate  that  even  in  this  case  the  partial  hydrophobic  protection  which  might  arise  in 
the  domain  swapped  model  can  enhance  the  conversion  rate.  Fig.  (8)  also  illustrates  how 
the  relative  softening  of  the  loop  can  enhance  the  rate  by  again  reducing  the  energy 
barrier.  This  assumes  that  the  bending  stiffness  of  the  loop  with  proline  or  mutant 
participate  in  the  reaction  coordinate  involved  in  conversion. 

Fig.  (9)  below  shows  a  preliminary  examination  of  the  B-factor  for  the  domain  swapped 
prion  trimer  model.  Uniform  substitutions  at  the  three  proline  containing  loops  were 
made,  and  the  LH(3Hs  were  restrained  while  the  loops  were  left  free  to  fluctuate.  While 
there  is  no  systematic  favoring  of  one  mutant  residue  over  another,  it  is  clear  that  the  WT 
DSTP  produces  a  B-factor  at  or  near  the  bottom  for  each  of  the  three  loops.  This 
supports  the  notion  that  softer  loops  (with  proline  removed)  can  search  more  quickly. 

We  will  study  the  relative  hydrophobic  stabilization  of  the  core  domain  swapped  loop 
region  by  two  approaches:  (i)  examination  of  explicit  water  exclusion  in  the  core,  and  (ii) 
energy  minimization  in  the  presence  of  salvation.  If  we  find  significant  energy  lowering 
in  the  latter  case,  we  will  attempt  to  carry  out  free  energy  perturbation  analysis  in  all 
atom  MD. 

Within  the  fibril  model  framework  discussed  earlier  in  the  report,  it  is  of  interest  to  study 
the  fatal  familial  insomnia  mutation  D178N,  which  requires  also  methionine  in  codon 
129.  We  note  that  for  the  fibril  arrangement  N|C|C that  there  is  potential  proximity 
between  residue  178  in  one  monomer  and  residue  129  in  another.  With  the  asparagine 
substitution,  the  formation  of  a  stabilizing  hydrogen  bond  between  the 
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Fig.  8:  Two-state  model  to  motivate  perturbation 
study  of  mutations  on  domain  swapped  hinge 
regions  of  the  model  prion  trimer.  Upper  panel 
illustrates  the  effect  of  mutation  to  relatively 
hydrophobic  residues  from  proline.  The  wild  type 
three  monomer  surface  is  modestly  destabilized 
(the  hinge  regions  are  in  the  solvent  exposed 
random  coil)  while  they  may  be  relatively 
hydrophobically  protected  in  the  DSTP  structure. 
As  a  result  the  energy  surface  ramping  reduces 
the  kinetic  barrier  to  oligomerization.  The  lower 
panel  illustrates  the  generic  effect  of  softening  by 
removing  the  proline  residue.  With  the 
assumption  that  the  librational  coordinates  of  the 
loop  are  part  of  the  reaction  coordinate,  the 
softening  flattens  each  surface  (wild  type  and 
infectious  oligomer)  and  thus  reduces  the  kinetic 
barrier  between  them. 
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Figure  9:  B  factor  of  domain  swapped 
prion  trimer  from  1  ns  all  atom 
simulations  of  mouse  prion.  The  main 
secondary  structures  are  restrained  as 
the  proline  bearing  loops  are  allowed  to 
fully  fluctuate.  In  collaboration  with  S. 
Yang. 


methionine  sulfur  and  the  terminal  side  chain  amide  of  the  asparagines  is  possible,  which 


14 


is  not  the  case  with  aspartic  acid.  The  conjecture  is  that  this  stabilizing  link  will  hold  for 
an  intermediate  on  the  path  to  conversion  with  templating  via  the  CLHfiH  of  NLH  H 
growth.  It  is  interesting  to  note  that  such  models  within  our  scheme  for  constructing 

fibrils  will  be  less  prone  to  fibril  formation.  We  will  construct  such  models  to  examine 

whether  they  can  account  for  enhanced  conversion  kinetics. 

Examining  the  possible  role  of  the  LH/3H  in  prion  conversion.  We  are  interested  in  two 
models  here.  First,  as  shown  in  the  figure  below,  the  novel  tri-octarepeat  motif  first 
identified  by  Zahn15  and  illustrated  in  Fig.  (10)  has  a  cross  section  that  matches 

beautifully  to  the  LH[)H  as  illustrated  in  Fig.  (11).  Attachment  to  the  beta  helix  may  be 

facilitated  by  hydrogen  bonds  between  glutamine  side  chains.  The  goal  of  the  modeling 
will  be  to  see  if  (a)  the  tri-octrepeat  attachment  can  indeed  occur  to  any  putative  LH(3H 
structures  at  the  N  or  C  terminus,  and  if  so  (b)  can  this  transmit  the  template  information 
to  the  N-terminal  region.  We  will  study  the  former  with  all  atom  MD  using  soft  restraints 
for  the  tri-octarepeat  and  hard  restraints  for  the  different  LH(3H  motifs;  we  will  examine 
bonding  for  all  orientations  of  the  LH(3H  structures  and  the  tri-octarepeat.  A  second 
possibility  we  wish  to  explore  is  that  the  CLH(3H  templates  NLH(3H  conversion  from 
either  fibrils  or  proteinase  resistant  oligomers.  In  this  picture  the  CLH|)H  anchored 
structures  helps  to  stabilize  less  proteinase  resistant  oligomers. 

Modeling  of  prion  strains  with  domain  swapped  oligomers.  We  will  construct  models 
and  explore  the  possibility  that  prion  strain  is  encoded  in  (i)  number  of  monomers  in  a 
domain  swapped  oligomer,  (ii)  splay  of  LHpHs  in  a  domain  swapped  oligomer,  and  (iii) 
arrangements  of  CLH|)Hs  and  NLH(3Hs  in  fibril  repeat  units  which  because  of  the 
terminal  properties  do  not  promote  fibril  growth. 

We  have  already  tested  the  stability  of  a  domain  swapped  dimer  of  the  PrP106  protein  in 
direct  analogy  to  our  earlier  work  on  the  domain  swapped  trimer  (see  Fig.  (12) ).  As  with 
the  trimer  there  is  enhanced  overall  stability  due  to  entanglement  via  domain  swapping. 
This  suffers,  however,  from  the  same  nonoptimal  threading  of  the  NLHpH  noted  earlier 
in  the  proposal.  Accordingly,  assuming  we  find  the  most  optimal  NLH(3H  threading,  we 
will  study  other  oligomers  with  different  numbers  of  monomers  to  assess  generic  stability 
of  such  domain  swapped  oligomers. 
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Figure  10:  Ribbon  structure  of  tri- 
octarepeat  of  Ref.  (  15).  Glutamines  are 
highlighted  in  yellow.  Produced  by  VMD 


Figure  11:  Top  view  of  tri- 
octrepeat  motif  laid  on  CTP3  model 
from  the  prion  protein.  Produced  in 
RASTOP. 


Top  (N-terminus) 


Figure  12:  Ribbon  view  of  model 
domain  swapped  PrP106  dimer 
(alpha  helices  not  shown).  Produced 
in  RASTOP. 


Notice  in  the  side  view  that  the  individual  LH(3Hs  are  splayed.  This  opens  the  possibility 
that  there  are  stable  or  metastable  orientations  of  LH(3H  in  oligomers  which  are  splayed 
relative  to  the  overall  oligomer  axis,  which  provides  another  possible  conformation 
encoding  for  strains.  To  study  this,  we  will  (i)  carry  out  unbiased  long  time  MD  studies 
of  potential  oligomer  configurations  and  look  for  significant  residency  time  in  splayed 
states,  and  (ii)  employ  soft  restraints  together  with  a  systematic  sweep  of  splay  angles  to 
hunt  for  relatively  stable  starting  configurations  for  unrestrained  MD  runs. 

Reportable  Outcomes:  Awards,  Activities,  Presentations  and  Publications 

1.  Awards,  Activities,  and  Presentations  (Appendix  A) 

2.  Modelling  of  Oligomeric  Intermediates  in  prion  aggregation,  manuscript 

published  in  Biophysical  Journal.  (Appendix  B) 

3.  Modeling  of  insertion  of  Abeta  peptides  in  model  membranes,  manuscript 
published  in  Biophysical  Journal.  ( Appendix  C) 

4.  Dissertation  on  modeling  of  prions  and  Abeta.  Directions  to  online  copy  of 
dissertation  by  David  Mobley,  with  strong  review  chapter.  (Appendix  D) 

5.  Proposing  a  new  structure  for  the  minimal  infectious  unit  of  prion  disease:  the 

Domain  Swapped  Prion  Trimer  (DSTP).  manuscript  of  paper  published  in 
FASEB  Journal.  (Appendix  E) 

6.  Demonstration  of  the  necessity  of  the  membrane  in  exponential  in  vivo  growth 

of  prion  infectivity.  Manuscript  published  in  the  Letters  section  of  Biophysical 
Journal.  (Appendix  F) 

7.  Possible  mechanism  for  copper  inhibition  of  prion  conversion.  Manuscript 
published  in  the  Letters  section  of  Biophysical  Journal.  (Appendix  G) 

8.  Detailed  modeling  of  in  vitro  y east  prion  aggregation  Manuscript  of  paper 
published  in  Phys.  Rev.  E.  ( Appendix  H) 
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9.  Review  of  amyloidogenic  proteins  from  a  materials  perspective.  Manuscript 
published  in  Materials  Research  Bulletin.  (Appendix  I) 

10.  Yahoo!  News  Search  Results  for  Cox  Singh  "mad  cow"  (Appendix  J) 


Conclusions 


In  summary,  through  the  course  of  this  research  we  have  most  significantly  (i)  produced 
new  and  more  realistic  models  for  prion  membrane  associated  aggregation  kinetics;  (ii) 
modeled  insertion  of  the  A-beta  peptide  into  model  membranes  and  noted  a  possible 
correlation  with  inherited  forms  of  Alzheimer’s  disease  (while  eyeing  possible  relevance 
to  insertion  of  the  putative  transmembrane  segment  of  the  prion  protein);  (iii)  created  a 
new  domain  swapped  prion  trimer  model  which  rescues  the  beta  helical  trimer  model 
from  instability;  (iv)  shown  that  the  membrane  not  only  mediates  toxicity  of  the  prion 
infectious  units  but  also  exponential  growth;  (v)  noted  that  copper  can  protect  the  prion 
protein  against  conversion  to  the  infectious  form;  (vi)  demonstrated  that  a  possible  C- 
terminal  beta  helix  structure  can  potentially  explain  in  vitro  fibril  structures  coupled  with 
N-terminal  beta  helix  formation  and  novel  domain  swapping. 
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Simulations  of  Oligomeric  Intermediates  in  Prion  Diseases 


David  L.  Mobley,  Daniel  L.  Cox,  Rajiv  R.  P.  Singh,  Rahul  V.  Kulkarni,  and  Alexander  Slepoy* 
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ABSTRACT  We  extend  our  previous  stochastic  cellular  automata-based  model  for  two-dimensional  (areal)  aggregation  of 
prion  proteins  on  neuronal  surfaces.  The  new  anisotropic  model  allows  us  to  simulate  both  strong  /3-sheet  and  weaker 
attachment  bonds  between  proteins.  Constraining  binding  directions  allows  us  to  generate  aggregate  structures  with  the 
hexagonal  lattice  symmetry  found  in  recently  observed  in  vitro  experiments.  We  argue  that  these  constraints  on  rules  may 
correspond  to  underlying  steric  constraints  on  the  aggregation  process.  We  find  that  monomer-dominated  growth  of  the  areal 
aggregate  is  too  slow  to  account  for  some  observed  doubling-time-to-incubation-time  ratios  inferred  from  data,  and  so  consider 
aggregation  dominated  by  relatively  stable  but  noninfectious  oligomeric  intermediates.  We  compare  a  kinetic  theory  analysis  of 
oligomeric  aggregation  to  spatially  explicit  simulations  of  the  process.  We  find  that  with  suitable  rules  for  misfolding  of  oligomers, 
possibly  due  to  water  exclusion  by  the  surrounding  aggregate,  the  resulting  oligomeric  aggregation  model  maps  onto  our 
previous  monomer  aggregation  model.  Therefore  it  can  produce  some  of  the  same  attractive  features  for  the  description  of  prion 
incubation  time  data.  We  propose  experiments  to  test  the  oligomeric  aggregation  model. 


INTRODUCTION 

Prion  diseases  are  a  group  of  neurodegenerative  diseases 
including  bovine  spongiform  encephalopathy  (BSE)  in 
cattle,  scrapie  in  sheep  and  goats,  chronic  wasting  disease 
in  deer  and  elk,  and  kuru  and  Creutzfeldt-Jakob  disease 
(CJD)  in  humans.  These  diseases  came  to  the  forefront  after 
BSE  reached  epidemic  proportions  in  Great  Britain  in  the 
early  1990s,  and  it  was  later  shown  that  transmission  of  BSE 
to  humans  can  lead  to  new  variant  CJD  (vCJD)  in  humans 
(Bruce  et  al.,  1997;  Hill  et  al„  1997;  Scott  et  al.,  1999). 

Prion  diseases  are  unusual  in  that  they  appear  to  be  caused 
by  infection  with  some  minimal  infectious  “seed”  of 
misfolded  prion  protein,  which  alone  may  be  able  to  cause 
disease  by  catalyzing  further  misfolding  and,  in  many  cases, 
aggregation  of  the  prion  protein.  These  aggregates  are 
typically  amyloidlike  fibrils  or  amyloid  plaques  (Caughey, 
2000).  The  infectious  agent  is  unusually  hard  to  eliminate  by 
various  methods  including  ultraviolet  irradiation,  suggesting 
it  contains  no  nucleic  acid  and  rather  only  protein,  the  so- 
called  “protein-only”  hypothesis  in  prion  diseases  (Weiss- 
mann  et  al.,  2002). 

In  the  case  of  CJD,  a  sporadic  form  of  the  diseases  also 
exists,  occurring  more  or  less  randomly  worldwide  with  an 
incidence  of  about  one  in  a  million  people  per  year.  It  has 
been  suggested  that  this  incidence  is  due  to  the  very  rare 
event  of  nucleating  the  minimal  infectious  seed  by  chance  in 
a  healthy  individual  (Come  et  al.,  1993). 

Developing  an  understanding  of  these  diseases  is 
important  because,  for  one,  they  are  invariably  fatal.  To 
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date,  no  treatment  exists.  Additionally,  it  is  not  yet  clear  how 
large  the  vCJD  epidemic  in  humans  will  be;  an  understand¬ 
ing  of  the  disease  process  is  important  to  be  able  to  guide  the 
search  for  treatment  ideas. 

In  many  cases,  prion  diseases  result  in  large,  up-to- 
micron-scale  plaques  in  the  brains  of  people  and  animals 
with  these  diseases.  They  also  involve  vacuolization  or 
spongiform  change  in  the  brain  due  to  death  of  neurons 
(Scott  et  al.,  1996).  Additionally,  the  normal  fomi  of  the 
prion  protein  (known  as  PrPc)  has  long  been  known  to 
misfold  and  aggregate  in  vitro  when  catalyzed  by  the 
presence  of  a  misfolded  prion  protein  (PrPSc)  seed  (Come 
et  al.,  1993).  Together,  these  observations  have  suggested  to 
some  that  the  aggregation  process  itself  may  be  important  in 
these  diseases  (Come  et  al.,  1993;  Masel  et  al.,  1999).  It  has 
also  been  suggested  that  the  rate-limiting  step  in  aggregation 
is  nucleation  of  an  appropriate  seed,  thus  the  rapid 
aggregation  in  the  seeded  case  described  above  (Come 
et  al.,  1993). 

Another  fact  which  may  be  important  to  this  issue  is  that 
the  prion  protein  is  normally  GPI-anchored  to  the  cell 
surface.  Aggregation  in  vitro  as  mentioned  above  is  observed 
in  solution  rather  than  in  the  presence  of  the  GPI  anchor  on 
a  cell  surface,  leaving  the  possibility  that  the  aggregation 
process  in  vivo  is  different. 

Aggregation  models  developed  to  explore  the  aggregation 
process  in  prion  disease  include  one-dimensional,  fibrillar 
aggregation-and-fission  models  (Masel  et  al.,  1999;  Slepoy 
et  al.,  2001),  since  aggregates  grown  in  vitro  are  typically 
seen  to  be  fibrillar.  Additionally,  our  earlier  work  suggested 
that  an  areal  aggregation  model  could  explain  certain  other 
properties  of  the  diseases  (Slepoy  et  al.,  2001).  By  areal 
aggregation,  we  mean  two-dimensional  aggregation  in 
a  relatively  regular  array,  probably  on  the  cell  surface  due 
to  GPI  anchoring,  in  contrast  to  the  one-dimensional,  fibril¬ 
lar  aggregation  observed  in  vitro,  and  also  in  contrast  to 
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two-dimensional  plaques  of  crossing  fibrils  which  can 
be  observed  in  vivo.  This  earlier  model  is  attractive  in  that 
it  can  provide  a  simple  explanation  for  the  long  lag  phase 
which  is  sometimes  observed  in  growth  of  the  amount  of 
infectious  material  in  the  brain.  This  lag  phase  of  little  or 
no  growth  is  followed  by  a  doubling  phase  with  a  short 
characteristic  doubling  time.  Additionally,  our  earlier  model 
provides  a  possible  explanation  of  some  of  the  difference 
between  infectious  and  sporadic  forms  of  CJD  (Slepoy  et  al., 
2001).  In  later  work,  we  used  this  model  to  explain  and  fit 
experimental  dose  incubation  curves  (Kulkami  et  al.,  2003). 

However,  there  were  drawbacks  to  the  earlier  aggregation 
model  we  proposed.  First,  no  such  areal  aggregates  had  so  far 
been  observed.  Second,  the  fissioning  essential  to  the  model 
would  involve  breaking  of  strong  bonds  between  the  proteins, 
probably  bonds  between  /3-sheets  (Serag  et  al.,  2002). 

More  recent  experimental  work  found  two-dimensional 
areal  aggregates  of  prion  protein  produced  during  the 
purification  process.  These  aggregates  were  examined  under 
electron  microscope  and  found  to  consist  of  trimeric  or 
hexameric  subunits.  These  subunits  are  linked  together  in 
a  regular  array,  possibly  by  their  N-terminal  sugars  or  a  weak 
protein-protein  interaction  (Wille  et  al.,  2002). 

This  suggested  we  should  modify  our  earlier  model  and 
attempt  to  reproduce  this  aggregate  morphology.  We  thought 
of  two  basic  schemes  for  growing  aggregates  of  this  sort: 

1 .  Growing  the  aggregate  outward,  monomer  by  monomer, 
from  an  initial  seed,  or 

2.  Oligomeric  intermediates  (possibly  very  flexible  and  of 
unstable  shape),  which  form  on  their  own  in  solution  and 
are  only  catalyzed  into  stably  misfolding  in  the  presence 
of  an  existing  misfolded  seed. 

Some  evidence  in  favor  of  case  2  has  already  been 
produced.  Monomers  of  yeast  prion  can  form  intermediates 
if  left  to  stand,  which  allows  aggregation  to  proceed  at  an 
initial  faster  rate  when  catalyzed  by  addition  of  a  seed  (Serio 
et  al.,  2000).  Additionally,  the  conformation-dependent  im¬ 
munoassay  developed  by  Safar  et  al.  (2002)  detects  both 
protease-sensitive  and  protease-resistant  PrPSc.  In  hamster 
brains,  sensitive  PrPSc  is  observed  earlier,  followed  by 
resistant  PrPSc.  This  could  correspond  to  case  2  above,  where 
the  sensitive  PrPSc  is  the  intermediates  that  are  not  yet  stably 
misfolded  and  the  resistant  PrPSc  is  stably  misfolded  inter¬ 
mediates. 

Work  here  has  been  done  to  further  explore  these  two 
potential  modifications  of  our  earlier  model  to  examine 
whether  they  retain  the  same  features  and  if  additional 
insight  can  be  gained. 

It  is  important  to  note  that  even  if  areal  aggregation  is  not 
important  to  the  time  course  of  these  diseases,  the  aggregates 
observed  by  Wille  and  co-workers  have  already  provided 
insight  into  the  structure  of  the  misfolded  prion  protein 
(Wille  et  al.,  2002).  Theoretical  modeling  may  be  able  to 
place  further  constraints  on  the  protein  or  subunit  structure 


necessary  to  reproduce  these  aggregates,  and  hence  provide 
valuable  information  because  these  aggregates  can  form, 
even  if  they  are  not  important  to  the  disease  progression. 

BASICS  OF  OUR  MODEL 

Here  we  explore  the  two  basic  schemes  suggested  above  for 
growing  aggregates  like  those  observed  by  Wille  et  al. 
(2002).  To  do  so,  we  use  a  modification  of  our  earlier  model. 
Therefore  a  recap  of  common  features  of  these  models  is 
useful. 

These  models  are  stochastic  cellular  automata  models, 
meaning  that  they  take  place  on  a  lattice  with  probabilistic 
interaction  and  diffusion  rules  governing  the  progression  of 
the  system.  In  this  case,  sites  on  the  lattice  are  either 
occupied  by  individual  prion  proteins,  or  water  (empty,  in 
the  simulation).  The  protein  form  at  a  site  can  also  vary  from 
PrPc  to  PrPSc. 

Rules  vary  depending  on  the  model  being  explored,  but  the 
basic  procedure  is  the  same.  For  every  simulation  step,  which 
represents  a  small  amount  of  time,  we  allow  proteins  and  any 
aggregates  to  diffuse  a  small  amount  on  the  lattice  (each 
object  has  a  probability  1  /(size)112  of  moving  one  lattice  site 
in  a  given  step).  Then  we  look  at  every  protein  in  the  lattice 
and  update  its  state  according  to  the  rules.  For  example,  in  our 
original  model,  the  conformation  of  an  individual  prion 
protein  is  determined  solely  by  its  number  of  neighboring 
prion  proteins,  and  this  can  vary  from  step  to  step.  After 
doing  this,  we  add  more  normal  prion  monomers  to  replace 
any  that  converted  to  PrPSe.  This  is  due  to  the  assumption  that 
this  process  would  be  taking  place  in  a  small  area  on  a  cell, 
and  the  normal  prion  monomers  would  be  added  by  the  cell 
or  diffuse  in  from  other  locations  on  the  cell  surface  to  keep 
the  monomer  concentration  relatively  constant. 

GROWTH  VIA  MONOMER  ADDITION 

First,  case  1  from  above  was  explored.  Simple  rules  were 
developed  (Fig.  1)  which  can  reproduce  aggregates  similar  to 
those  observed  by  Wille  et  al.  (2002).  It  is  important  to  note 
that  although  the  rules  were  designed  to  reproduce  such 
aggregates,  most  modifications  of  these  rules  could  not  do 
so.  This  means  that  the  rules  provide  some  constraints  on 
the  protein-protein  interactions  necessary  to  reproduce  such 
aggregates.  Also,  for  the  purposes  of  this  model,  we  are 
assuming  the  subunits  are  hexameric,  but  the  corresponding 
model  for  trimeric  intermediates  is  actually  much  simpler 
than  this  model  and  will  produce  similar  results.  Details  of 
the  algorithm  for  this  model  are  covered  in  Fig.  2. 

The  rules  are  as  follows.  The  simulation  begins  with 
a  single  hexagonal  subunit  consisting  of  six  misfolded 
monomers  (light  gray  hexagons  in  Fig.  1)  which  stick  some 
of  their  residues  into  an  adjacent  site,  excluding  anything 
else  from  occupying  that  site  (black).  Healthy  monomers 
(light  gray  spheres )  can  then  attach  via  a  sugar-bond  or  other 
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FIGURE  1  Simple  rules  for  monomer-by-monomer  growth  of  aggregates  like  those  observed.  Some  possible  rules  can  be  excluded,  thus  these  rules  give 
insight  into  how  the  proteins  involved  must  be  interacting  with  one  another,  (a)  The  initial  seed  consists  of  six  misfolded  monomers  (light  gray  hexagons ) 
surrounding  a  central  region  (black)  which  is  occupied  by  some  residues  sticking  into  it  from  the  adjacent  six  sites,  (b)  A  healthy  monomer  (light  gray  sphere)  can 
move  adjacent  to  a  misfolded  one  and  attach  via  a  sugar  bond  or  other  weak  protein-protein  interaction  (proteins  sugar-bonded  are  colored  dark  gray).  This 
cannot  happen  if  the  monomer  moves  into  the  site  between  two  misfolded  proteins,  (c)  Subsequent  monomers  can  move  next  to  the  attached  one  and  misfold  and 
begin  to  form  a  new  hexamer.  Residues  from  the  two  stick  into  the  black  region,  preventing  anything  else  from  moving  there.  This  cannot  happen  if  the  second 
monomer  is  adjacent  to  the  existing  hexamer;  this  would  produce  irregular  aggregates  unlike  those  observed  by  Wille  and  co-workers  (Wille  et  al.,  2002).  (d) 
The  forming  hexamer  can  grow  and  finish  via  subsequent  monomer  addition,  (e)  Continue  a-d  for  a  long  time,  and  an  aggregate  like  the  one  shown  can  form. 


weak  protein-protein  interaction  to  this  subunit  (dark  gray 
spheres  to  dark  gray  hexagons)  but  only  radially  outward 
from  a  monomer  in  the  initial  hexamer.  Additional  mono¬ 
mers  moving  adjacent  to  the  attached  monomer  can,  together 
with  it,  misfold  but  only  if  the  second  monomer  does  not  also 
neighbor  the  original  hexamer.  Then  additional  monomers 
can  attach  to  this  forming  hexamer,  allowing  it  to  complete. 
Repeating  this  process  many  times  can  produce  mostly  reg¬ 
ular  aggregates  with  some  holes,  similar  to  those  observed. 

The  rules  are  also  probabilistic:  above,  “can”  means  that 
some  fraction  of  the  time  the  event  occurs.  These  prob¬ 
abilities  can  be  changed  in  the  simulation  and  give  different 
growth  rates,  but  the  same  essential  features  and  scaling  as 
described  below. 

If  this  is  in  fact  how  these  aggregates  are  forming,  we  find 
out  about  the  orientation  of  monomers  within  a  hexagonal 
subunit.  We  find,  as  mentioned  in  the  discussion  of  the  rules 
above,  that  the  N-terminal  sugars  or  attachment  sites  must 
stick  radially  outward  from  each  monomer  in  a  hexagonal 
subunit  (Fig.  1  h).  This  is  in  agreement  with  the  hexagonal 
structure  proposed  by  Wille  and  co-workers  (Wille  et  al., 
2002).  Additionally,  we  find  that  no  such  regular  aggregates 
can  be  produced  unless  the  monomer  attaching  to  a  pre¬ 
viously  attached  monomer  (Fig.  1  c)  can  only  attach  if  it  is 
not  adjacent  to  an  existing  hexamer.  This  seems  to  indicate 
that  the  other  spaces  must  be  occupied  by  residues  from  the 
existing  hexamer,  preventing  attachment  in  those  sites. 

This  model  can  also  reproduce  gaps  in  aggregates  as 
observed.  In  this  model  gaps  are  due  to  variations  of  the 
growth  rate  from  average  for  part  of  the  aggregate,  causing 
several  parts  of  the  aggregate  to  grow  apart  and  then  rejoin 
after  leaving  a  gap. 

One  reason  for  developing  this  model  was  to  see  if  it 
would  capture  the  same  features  of  the  disease  as  our  original 


model.  Our  original  model  explained  the  difference  be¬ 
tween  the  lag  phase  and  the  doubling  phase  by  suggesting 
that  the  doubling  phase  is  initiated  when  aggregates  begin 
to  fission,  then  regrow  to  a  certain  fissioning  size  and 
break  again.  Key  to  this  explanation  is  our  result  that 
aggregation  speeds  up,  so  that  the  time  for  an  aggregate  to 
double  in  size  from  half  its  fission  size  to  its  fission  size  is 
much  less  than  the  time  for  it  to  get  from  its  initial  size  to 
its  fissioning  size. 

To  see  if  this  model  could  produce  the  same  separation  of 
lag  and  doubling  phases,  we  examined  the  aggregate  growth 
rate  as  a  function  of  size  in  this  model  (Fig.  3)  and  found  it 
speeds  up  only  slowly.  Naively,  one  would  expect  the 
growth  rate  to  be  roughly  proportional  to  the  square  root 
of  the  size,  as  the  growth  rate  is  proportional  to  the 
circumference  of  the  aggregate,  which,  assuming  a  circular 
aggregate,  is  2m\  The  size  of  the  aggregate  is  proportional  to 
the  area,  iti-2,  so  the  radius  is  proportional  to  the  square  root 
of  the  size  and  thus  the  rate  proportional  to  the  square  root 
of  the  size.  To  a  good  approximation,  the  growth  rate  ob¬ 
served  here  is  well-fit  by  an  offset  plus  a  term  proportional  to 
(size)112,  as  expected. 

In  this  simple  picture,  one  can  calculate  the  ratio  of  the 
doubling  time  to  the  lag  time.  The  lag  time  is  the  time  to  go 
from  the  initial  size,  say  size  0  for  simplicity,  to  size  n\  the 
doubling  time  from  size  n/2  to  size  n.  Integrating  the  rate  to 
get  the  times  and  taking  the  ratio  we  find  tdoub/fiag  = 
1  -  1/V2  or  —0.293.  This  means  that  this  model  cannot 
produce  such  a  large  separation  between  lag  and  doubling 
times  as  our  earlier  model  could,  at  least  not  without  further 
modification. 

This  also  indicates  that  if  there  is  a  lag  phase  and  if  the 
difference  between  it  and  the  doubling  phase  is  due  to 
acceleration  of  aggregation,  this  picture  is  not  sufficient  and 
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FIGURE  2  Flow  chart  of  sim¬ 
ulation  for  monomer  addition 
model.  We  typically  use  P3  = 
0.2;  we  tried  a  variety  of  differ¬ 
ent  values  for  this  and  values 
near  0.2  seem  to  produce  the 
most  regular  aggregates.  We 
also  typically  use  P$  —  0.9.  This 
is  not  important  and  roughly  sets 
the  simulation  timescale.  Also, 
for  our  statistics,  we  typically 
average  >1000  such  runs  as  the 
one  described  here. 


something  more  like  case  2,  growth  from  intermediates,  may 
be  a  better  representation  of  the  disease  process. 

GROWTH  VIA  INTERMEDIATES 

In  this  case,  aggregation  is  assumed  to  be  the  assembly  of 
independent  hexameric  intermediates  into  a  larger  areal 


aggregate.  The  intermediates  themselves  are  not  misfolded 
but  only  misfold,  in  this  model,  when  they  either  aggre¬ 
gate  with  an  existing  misfolded  seed,  or  come  together  in 
such  a  way  that  they  can  misfold  and  form  a  new  stable 
seed.  In  this  way,  the  model  works  essentially  just  like  the 
model  of  Slepoy  et  al.  (2001),  except  now  hexameric  inter¬ 
mediates  are  playing  the  role  of  monomers  (Fig.  4).  As 
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FIGURE  3  Growth  rate  (change  in  aggregate  size  per  step)  as  a  function 
of  size  for  seeded  areal  aggregation  in  the  monomer  growth  model.  Growth 
rate  goes  as  the  square  root  of  the  size  with  an  offset,  which  was  as  expected 
for  this  model. 


mentioned  above,  there  is  some  evidence  that  intermediates 
greatly  increase  aggregation  rate  in  studies  of  yeast  prions,  so 
this  emphasis  on  the  importance  of  intermediates  may  be 
reasonable. 

To  be  able  to  map  this  model  back  into  our  old  model, 
though,  we  need  to  know  how  the  intermediate  concentration 
depends  on  monomer  concentration.  And  this  is  not  obvious. 
So  a  simulation  was  developed  to  explore  how  the  con¬ 
centration  of  hypothetical  hexameric  intermediates  would 
depend  on  monomer  concentration.  Again,  here  we  are 
assuming  the  intermediates  are  hexameric  but  we  can  easily 
modify  the  model  to  accommodate  trimers. 


To  get  at  the  concentration  of  intermediates,  it  was  as¬ 
sumed  that  two  monomers  have  a  probability  Px  of  begin¬ 
ning  a  new  hexameric  subunit  when  they  come  into  contact 
(see  Fig.  5).  This  new  subunit  can  grow  by  addition  of 
monomers  when  they  move  into  appropriate  positions 
(changing  this  probability  does  not  affect  the  outcome  of 
the  simulation,  only  the  timescale,  so  it  was  set  to  1). 
However,  this  growth  process  competes  with  a  “dissolving” 
process  by  which  a  monomer  that  is  part  of  an  intermediate 
but  only  has  one  neighboring  monomer  can  break  off  with 
a  probability  P3.  Thus  the  end  destiny  of  any  intermediate 
that  begins  is  either  to  form  a  complete  hexameric  in¬ 
termediate,  in  which  case  it  can  persist,  or  to  dissolve  com¬ 
pletely.  Details  of  the  algorithm  for  this  model  are  shown  in 
Fig.  6. 

This  dissolving,  or  reversibility,  was  included  because  it 
was  not  obvious  that  at  low  monomer  concentrations,  one 
would  expect  a  reasonable  formation  rate  of  intermediates  via 
this  mechanism.  It  was  initially  thought  that  at  concentrations 
below  something  on  the  order  of  P3,  breaking  would 
dominate  and  the  formation  rate  of  intermediates  would  be 
almost  zero.  First,  the  simulation  that  was  developed  was 
used  to  examine  the  dependence  of  time  for  intermediate 
formation  as  a  function  of  monomer  concentration  (Fig.  7).  It 
was  found  that  at  high  monomer  concentration,  the  time  to 
form  an  intermediate  scales  between  1/c  and  1/c2  (c  is 
concentration).  This  is  because  the  likelihood  of  starting  an 
intermediate  scales  as  the  dimer  concentration  ( 1/c2),  whereas 
the  time  to  add  monomers  to  it  scales  as  1/c.  On  the  other 
hand,  at  very  low  monomer  concentration,  the  time 
asymptotically  approaches  1/c6.  This  is  due  to  the  fact  that 
at  these  concentrations,  dissolving  dominates  and  it  is  only  in 
the  very  rare  event  that  six  monomers  are  in  the  same  place  at 


FIGURE  4  ( a )  As  in  Slepoy’s  model  (Slepoy  et  al.,  2001),  subunits  were 
healthy  monomers  ( light  gray  spheres )  aggregating  with  misfolded 
monomers  ( dark  gray  hexagons)',  (b),  subunits  are  hexagonal  intermediates 
( light  gray/dark  gray)  aggregating  with  misfolded  hexagonal  structures 
(< medium  gray /dark  gray).  In  both  cases,  the  aggregation  process  and  kinetics 
ought  to  be,  and  indeed  are,  similar. 
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FIGURE  5  Rules  for  the  formation  of  intermediates.  Note  that  growth  and 
dissolving  compete,  so  that  any  intermediate  eventually  either  becomes 
a  complete,  stable  hexagon  or  dissolves  back  into  monomers.  ( a )  Two 
monomers  have  a  probability  P\  of  joining  to  begin  a  new  intermediate, 
which  is  not  yet  stably  misfolded.  Black  represents  a  region  blocked  by  some 
of  their  residues.  ( b )  This  can  grow  by  addition  of  monomers  to  either 
“end.”  After  attaching,  the  monomer  sandwiched  between  the  other  two  has 
two  neighbors  and  is  not  allowed  to  break  off,  whereas  the  ones  with  only 
one  neighbor  can.  (c)  A  monomer  with  only  one  neighboring  monomer  has 
a  probability  P3  of  breaking  off  in  a  given  step.  This  competes  with  the 
growth  process,  (d )  Continuing  addition  of  monomers  can  result  in  a  finished 
hexameric  intermediate  where  every  monomer  has  two  neighbors  and  is  safe 
from  breaking  off. 
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FIGURE  6  Flow  chart  of  sim¬ 
ulation  for  the  formation  of 
intermediates.  Note  that  Px  we 
vary  for  different  runs,  P2  we 
typically  set  to  1  (it  sets  the 
simulation  timescale  and  is  un¬ 
important),  and  P3  we  also  vary. 
Finished  intermediates  are  re¬ 
moved  so  that  we  can  run  to  a 
larger  number  of  finished  inter¬ 
mediates  without  the  lattice  get¬ 
ting  clogged.  Here,  also,  we 
typically  average  >1000  trials 
for  good  statistics. 


almost  the  same  time  that  an  intermediate  can  finish.  The 
probability  of  that  scales  as  1/c6. 

It  is  interesting  to  note  that  the  beginning  of  the  transition 
between  high  concentration  behavior,  where  most  inter¬ 
mediates  successfully  become  complete,  and  low  concen¬ 
tration  behavior,  where  only  a  lucky  few  do,  begins  at 
a  concentration  on  the  order  of  the  breaking  probability,  P 3. 
This  suggests  that  if  the  strength  of  bonds  between  in¬ 
termediates  could  be  weakened  somehow,  the  biological 
number  of  intermediates  could  be  drastically  decreased  by 
pushing  biological  monomer  concentrations  into  the  1/c6 
regime. 

The  goal,  however,  was  to  determine  the  dependence 
of  the  intermediate  concentration  on  monomer  concentra¬ 
tion.  This  just  provided  a  formation  rate,  and  the  functional 


form  was  uncertain.  So  another  sort  of  result  was  exam¬ 
ined,  wherein  we  began  examining  behavior  of  the  system 
as  a  function  of  time,  and  measured  the  number  of  differ¬ 
ent  partial  intermediates  (two  monomers,  . .  .five  mono¬ 
mers,  hexameric  intermediates).  We  first  examined  the  case 
with  no  breaking  (P3  =  0)  to  check  our  results,  because 
it  is  relatively  easy  to  work  out  kinetics  in  that  case.  A 
sample  of  one  of  these  plots  is  shown  in  Fig.  8,  with 
symbols  as  data  points  and  solid  lines  as  approximate 
kinetics  fits.  It  is  important  to  note  that  in  this  case,  and  in 
the  case  of  nonzero  breaking  probability,  the  number  of 
dimers,  trimers,  tetramers,  and  pentamers  reaches  equilib¬ 
rium  relatively  quickly  and  then  the  hexamer  number 
begins  to  grow  linearly  at  a  rate  equal  to  the  rate  of  dimer 
formation. 
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Soft  hexamer  aggregation:  Steps  vs.  cone 


p3=0,01;  view  of  transition  region 


FIGURE  7  Simulation  steps  (time)  to  form  a  hexameric  intermediate  as 
a  function  of  monomer  concentration.  Log-log  scale.  Note  the  broad 
transition  to  dissolving-dominated  behavior  at  low  concentration.  The 
transition  actually  continues  to  even  lower  concentration  than  can  be  seen 
here.  At  very  low  concentration  the  time  eventually  scales  as  1/c6.  Standard 
deviations  fall  within  the  size  of  the  data  points  on  this  plot. 


Number  of  each  size  versus  time 

To  1000  finished  hexes;  pl=0.1,  p3=0.001,  c=0.25% 


FIGURE  9  Number  of  each  size  as  a  function  of  time  (simulation  steps) 
with  nonzero  breaking.  Compare  to  Fig.  6;  note  that  the  number  of 
intermediates  reaches  equilibrium  faster  and  at  smaller  numbers,  but  that  the 
hexamer  number  still  grows  linearly  at  long  times. 


Sample  results  with  nonzero  breaking  are  shown  in  Fig.  9. 
These  results  are  qualitatively  similar,  except  the  number  of 
pre-intermediates  that  persists  is  much  lower.  In  the  high- 
breaking  limit,  the  very  low  level  of  intermediates  dem¬ 
onstrates  that  either  a  potential  intermediate  gets  “lucky” 
and  quickly  forms  an  intermediate,  or  it  dissolves  back  to 
monomers,  leaving  few  dimers,  trimers,  and  so  on. 

The  kinetics  equations  we  can  write  down  to  describe  this 
simulation  are  relatively  simple.  With  rnm  as  the  rate  constant 
for  forming  m-mers  from  n-mers,  and  bnm  as  the  rate  of 
breaking  n-mers  into  m-mers  plus  monomers,  we  can  write: 


Number  of  each  size  versus  time 

To  200  finished  hexes;  pi  =0.1 ,  p3=0,  c=0. 1  % 


FIGURE  8  Number  of  each  size  as  a  function  of  time  (simulation  steps), 
with  zero  breaking.  Note  that,  at  long  times,  intermediates  reach  equilibrium 
and  the  hexamer  number  begins  growing  linearly  with  time.  Points  are 
simulation  data  points;  solid  lines  (mostly  overlapping  points)  are 
approximate  kinetics  results. 
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- =  t'45 [c  1  ] [c4]  -  /'seN [c5]  ~  b54[c5\  (5) 

^T  =  rX[cl][c5].  (6) 


Since  we  know  that  the  hexamer  number  grows  linearly  at 
steady  state  and  all  of  the  other  concentrations  are  un¬ 
changing,  we  can  greatly  simplify  the  above  kinetics  by 
looking  at  the  steady  state  only.  We  can  work  backward  from 
the  steady-state  behavior  of  the  hexamers  to  find  the 
dependence  of  the  steady-state  rate  of  hexamer  formation 
on  the  different  kinetic  parameters  and  ultimately  on  the 
monomer  concentration. 

This  straightforward  kinetics  analysis  produces  the 
equilibrium  result 
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where  m  is  the  slope  at  equilibrium  of  the  hexamer  formation 
rate. 

The  constants  in  our  simple  result  for  m ,  above,  can  be 
measured  from  our  simulation.  However,  our  simulation 
does  not  necessarily  reproduce  what  these  constants  would 
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be  in  a  biological  system.  So  it  is  difficult  to  say  exactly  what 
the  rate  of  intermediate  formation,  m ,  would  be  in  a  real 
system.  However,  it  is  nevertheless  useful  to  know  the 
functional  form  of  its  dependence  on  the  monomer  con¬ 
centration. 

The  result  that  the  hexamer  number  begins  growing 
linearly  eventually  is  independent  of  monomer  concentra¬ 
tion.  This  is  important  because  it  means  some  hexamers  can 
form  given  these  simple  rules  even  if  breaking  dominates. 
Given  that  result,  it  seems  safe  to  assume  that  if  hexameric 
intermediates  are  stable,  some  will  form  in  biological 
systems. 

In  our  model,  the  hexamer  number  grows  linearly 
indefinitely,  which  is  obviously  unrealistic  biologically. 
The  reason  for  this  is  that  we  include  no  mechanism  to 
remove  finished  hexamers.  Realistically,  they  would  be 
cleared  from  the  body  somehow.  They  could  be  endocytosed 
from  the  cell  surface  and  degraded  via  the  proteasome 
mechanism  or  some  other  pathway.  Additionally,  any 
hexamers  being  taken  up  into  aggregates  would  reduce  this 
number.  Regardless,  realistically  the  number  should  stabilize 
at  some  fixed  value  determined  by  the  balance  of  the 
clearance  rate  and  the  formation  rate. 

With  the  result  that  some  hexamers  fomi  even  at  low 
monomer  concentrations  (and  more  would  form  if  they  are 
trimers),  a  model  was  developed  where  now  hexameric 
intermediates  occupy  a  single  cell  on  the  lattice  (equiva¬ 
lently,  these  could  be  trimeric  intermediates).  This  model, 
described  below,  largely  maintains  the  same  attractive 
features  of  the  original,  showing  that  if  areal  aggregation  is 
the  explanation  for  these  features,  as  we  suggested,  this 
aggregation  could  be  of  hexameric  intermediates. 

Part  of  our  basis  for  this  model  is  the  observation  that  the 
intermediates  are  not  yet  stably  misfolded  since  formation  of 
intermediates  in  studies  of  yeast  prions  does  not  lead  to 
a  change  in  circular  dichroism  results;  it  is  only  when  they 
aggregate  with  a  seed  that  they  stably  misfold  (Serio  et  al., 
2000).  This  also  is  justified  by  observing  that  if  intermediates 
were  stably  misfolded,  they  could  act  as  seeds  on  their  own, 
without  the  necessity  of  an  external  seed  initiating  the 
infection,  and  thus  there  would  be  no  difference  between 
sporadic  and  infectious  CJD.  Therefore,  for  aggregates 
consisting  of  misfolded  oligomers  like  those  observed  by 
Wille  et  al.  (2002),  intermediate  misfolding  must  be 
catalyzed  by  existing  aggregates  or  few-hexamer  misfolded 
oligomers.  We  hypothesize  that  the  mechanism  for  this  is 
intermediates  forming  bonds  to  an  existing  seed.  When 
solvent  is  excluded  locally  around  these  oligomers  and  their 
neighbors  include  a  misfolded  oligomer  or  aggregate,  they 
misfold.  The  important  point  is  that  it  is  solvent  exclusion 
around  an  intermediate  that  can  cause  it  to  misfold,  making 
this  a  very  rare  sporadic  event.  But  a  misfolded  seed  can  help 
this  process  by  providing  a  place  where  intermediates  bond, 
helping  the  solvent-exclusion  process.  These  rules  make  this 
model  essentially  identical  in  terms  of  kinetics  to  our  original 


model.  Details  of  the  algorithm  for  this  model  and  mapping 
are  shown  in  Fig.  10. 

However,  from  our  old  model  we  estimated  the  sporadic 
form  of  the  disease  could  have  a  peak  at  ~  1000  years,  given 
a  biological  concentration  of  10  In  our  new  model  we 
find  that  it  is  very  difficult  to  estimate  this  number  as  the 
scaling  of  the  time  as  a  function  of  monomer  concentration  is 
complicated.  It  was  hoped  that  this  model  would  give  a  result 
for  the  onset  of  sporadic  disease  that  could  be  compared  with 
the  time  for  onset  of  the  infectious  form  to  see  if  the  results 
were  consistent  with  the  roughly  l-in-106  incidence  of 
sporadic  CJD  that  we  earlier  pointed  out.  Unfortunately,  it  is 
difficult  for  our  model  to  give  a  concrete  answer  at  this  time 
as  the  answer  depends  too  much  on  the  value  of  the 
biological  monomer  concentration.  We  do  find,  however, 
that  the  power  law  used  previously  to  scale  the  sporadic  data, 
c-3,  is  a  lower  bound  on  the  separation.  That  is,  the  actual 
exponent  should  be  larger,  meaning  that  we  previously 
underestimated  the  separation  of  timescales.  Thus  although 
we  cannot  say  exactly  what  the  separation  of  timescales  here 
will  be,  we  can  say  that  it  will  be  greater  than  the  two  orders 
of  magnitude  that  we  previously  estimated. 

This  work  suggests  that  a  model  like  our  earlier  one, 
modified  to  involve  areal  aggregation  of  hexameric  or 
trimeric  intermediates,  could  maintain  the  same  attractive 
features  of  our  earlier  model  in  explaining  certain  aspects  of 
the  diseases.  However,  without  precise  knowledge  of  the 
biological  monomer  concentration  and  a  way  to  measure 
relevant  rate  constants,  it  is  difficult  to  make  numerical 
predictions  from  this  model. 

DISCUSSION 

Our  work  has  shown  that  both  in  the  case  of  monomer 
addition  to  a  seed,  and  in  the  case  of  growth  via  in¬ 
termediates,  it  is  possible  to  produce  aggregates  like  those 
observed  by  Wille  et  al.  (2002).  This  leaves  the  question  of 
how  such  aggregates  actually  grew.  If  areal  aggregation  is 
the  cause,  or  part  of  the  cause,  of  the  difference  between  lag 
and  doubling  times,  as  suggested  by  Slepoy  et  al.  (2001), 
then  our  work  suggests  that  intermediates  are  already  present 
in  vivo  before  aggregation. 

Our  work  has  also  shown  that  a  model  can  be  developed 
which,  with  suitable  parameters,  can  reproduce  areal  ag¬ 
gregates  like  those  actually  observed  while  maintaining  the 
same  features  of  our  original  model. 

Whether  or  not  areal  aggregation  is  actually  important  in 
these  diseases,  we  can  gain  insight  from  this  model.  If  the 
aggregates  observed  are  growing  via  monomer  addition,  we 
gain  some  constraints  on  the  structure  simply  from  our  rules. 
On  the  other  hand,  if  intermediates  are  important  to 
aggregation,  then  our  results  indicate  the  intermediate 
concentration  can  be  quite  important.  At  high  intermediate 
concentrations,  intermediates  form  relatively  fast.  However, 
at  low  intermediate  concentrations,  intermediate  formation 
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FIGURE  10  Flow  chart  for 
simulation  mapping  back  into 
our  original  model.  Here  we 
basically  have  free  monomers 
{ fMs ),  attached  monomers  that 
are  not  yet  stably  misfolded 
(i aM ),  and  monomers  that  have 
stably  misfolded  and  aggregated 
(//).  We  have  some  choice  of  a 
parameter,  Qbc  —  m.  This  model 
will  capture  the  features  of  our 
original  model  for  m  between  3 
and  6,  and  the  simulation  will 
proceed  in  exactly  the  same 
way.  We  compute  Nbc,  the  bond 
coordination  number,  with  Nbc 
=  nm  +  n.dM  +  (m  -  1)  X  nH, 
where  the  n  is  the  number  of 
neighboring  JMs,  and  so  on.  A^c, 
the  “hardening”  or  aggregating 
coordination  number,  is  given 
by  Nhc  =  /ifM  +  «aM  +  «H.  We 
refer  to  Qbc  as  the  bonding 
critical  coordination  number  and 
Qbc  as  the  “hardening”  critical 
coordination  number. 


timescales  as  1/c6.  This  result  is  exciting  because  it  suggests 
the  intermediates  as  a  target  to  prevent  aggregation.  Simply 
reducing  the  monomer  concentration  by  a  factor  of  2  would 
decrease  the  number  of  intermediates  by  a  factor  of  26  or  64. 
Within  our  model,  this  would  certainly  increase  the  aggrega¬ 
tion  time,  and  thus  slow  down  the  disease,  by  at  least  the 


same  factor.  For  a  disease  which  can  typically  incubate  for 
years,  this  obviously  would  be  a  great  advantage. 

In  this  case,  the  location  of  the  transition  between  low 
concentration  behavior  and  high  concentration  behavior  is, 
roughly  speaking,  set  by  the  probability  of  monomers 
breaking  off  from  an  intermediate  before  it  becomes  a  sta- 
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ble  hexamer.  Thus  if  this  probability  could  be  increased 
slightly — that  is,  the  bonds  between  monomers  could  be 
weakened  slightly — it  would  have  the  result  described 
above.  This  could  provide  an  explanation  for  one  experi¬ 
mental  observation.  Humans  have  a  methionine/valine 
polymorphism  at  codon  129  of  the  gene  for  the  prion 
protein.  To  date,  everyone  affected  by  vCJD  has  been 
methionine/methionine  homozygous.  This  effect  was  also 
seen  in  the  prion  disease  Kuru,  where  the  methionine/ 
methionine  genotype  was  associated  with  increased  sus¬ 
ceptibility  and  the  shortest  incubation  time  (Goldfarb, 
2002).  If  replacing  methionine  with  valine  weakened  the 
monomer-monomer  bonds  within  a  forming  intermediate 
and  reduced  intermediate  concentration,  this  could  have 
exactly  the  effect  described  above.  This  is,  however,  highly 
speculative,  but  as  Wille  et  al.  (2002)  refine  their  model  of 
the  oligomer  structures,  it  will  be  interesting  to  see  if  this 
residue  falls  in  the  region  important  to  bonding  between 
monomers. 

In  all,  our  work  shows  that  our  earlier  model  can  be 
extended  to  produce  aggregates  like  those  observed  in  vitro 
while  still  maintaining  its  attractive  features.  Our  work  also 
suggests  possible  mechanisms  for  formation  of  these 
aggregates.  If  the  aggregates  form  by  monomer  addition,  it 
constrains  protein  structure.  If  they  form  by  addition  of 
intermediates,  it  highlights  the  importance  of  bonds  within 
the  intermediates  as  a  target  for  possible  treatment  strategies. 
Our  model  suggests  that  an  experiment  to  measure  the 
biological  intermediate  concentration,  if  there  is  such 
a  concentration,  would  be  very  useful.  That  would  indicate 
whether  such  intermediates  are  present  at  a  high  enough 
concentration  to  be  important  biologically.  Additionally,  this 
work  suggests  that  experimentalists  should  check  and  see 
whether  reasonably-sized  aggregates  of  prion  protein  can  be 
found  in  vivo  on  the  cell  surface.  This  confinement  to  the  cell 
surface  conceivably  could  make  the  difference  between  the 
one-dimensional  fibrillar  aggregates  typically  observed  in 
vitro  and  two-dimensional  areal  aggregates  like  those 
suggested  by  the  model  of  Slepoy  et  al.  (2001).  Direct 
measurements,  or  detailed  simulations,  giving  the  strengths 
of  /3-bonds  between  monomers  compared  to  bonds  between 
subunits  would  be  very  useful. 

One  simple  way  to  experimentally  discern  between 
growth  via  monomers  or  intermediates  may  be  to  look  at 
high  resolution  at  the  boundary  of  actual  areal  aggregates.  If 
growth  is  by  monomers,  aggregates  will  form  with  mono¬ 
mer-scale  roughness  at  their  boundaries  (Fig.  1  e)  while  if 
growth  is  by  intermediates,  there  will  be  no  such  roughness 
(Fig.  4  b).  Experimentally,  the  absence  of  such  roughness 
would  not  prove  the  growth  via  intermediate  hypothesis 
because  incomplete  oligomers  at  the  edge  of  the  aggregate 
could  be  removed  in  the  purification  process,  possibly  by 
proteinase  K  digestion.  However,  the  presence  of  such 
roughness  would  certainly  suggest  that  monomer  growth  is 
important. 


A  more  general  scheme  for  experimentally  testing  the 
possible  role  of  intermediates  and  estimating  their  concen¬ 
tration  is  via  spin  labeling  (Hubbell  et  al.,  1998;  Columbus 
and  Hubbell,  2002).  Briefly,  a  small  molecule  with  a  free 
spin  can  preferentially  react  and  attach  to  cysteine  residues 
in  proteins.  Frequently,  these  residues  are  moved  around 
a  protein  via  mutagenesis  to  then  map  out  structures,  but  for 
these  purposes  a  less  refined  approach  is  required.  Since  the 
PrP  protein  already  possesses  cysteine  residues  at  the 
position  of  the  disulfide  bond,  the  spin  labels  can  attach 
there  (and  will  not  disrupt  the  disulfide  bond).  Then  the  spin- 
spin  interactions  will  produce  a  different  characteristic 
spectrum  for  monomers,  incomplete  intermediates,  and 
complete  intermediates,  in  particular,  with  a  progressive 
broadening  upon  moving  from  monomers  to  complete 
intermediates.  Since  the  spins  can  have  interactions  with 
other  spins  within  a  3-nm  sphere,  we  do  not  doubt  that  the 
broadening  will  be  observable.  Of  course,  since  the  spin 
labels  will  react  with  any  cysteines  present,  it  is  important 
to  carry  this  out  first  by  in  vitro  aggregation  experiments 
with  purified  prion  extracts.  This  will  help  to  identify 
the  conditions  which  can  lead  to  areal  aggregation  as 
observed  by  Wille  et  al.  (2000),  and  serve  as  an  existence 
proof  at  least  for  significant  oligomeric  intermediate  con¬ 
centrations. 
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ABSTRACT  Inspired  by  recent  suggestions  that  the  Alzheimer's  amyloid  (3  peptide  (A/3)  can  insert  into  cell  membranes  and 
form  harmful  ion  channels,  we  model  insertion  of  the  40-  and  42-residue  forms  of  the  peptide  into  cell  membranes  using  a  Monte 
Carlo  code  which  is  specific  at  the  amino  acid  level.  We  examine  insertion  of  the  regular  A/3  peptide  as  well  as  mutants  causing 
familial  Alzheimer’s  disease,  and  find  that  all  but  one  of  the  mutants  change  the  insertion  behavior  by  causing  the  peptide  to 
spend  more  simulation  steps  in  only  one  leaflet  of  the  bilayer.  We  also  find  that  A/342,  because  of  the  extra  hydrophobic 
residues  relative  to  A/340,  is  more  likely  to  adopt  this  conformation  than  A/340  in  both  wild-type  and  mutant  forms.  We  argue 
qualitatively  why  these  effects  happen.  Here,  we  present  our  results  and  develop  the  hypothesis  that  this  partial  insertion 
increases  the  probability  of  harmful  channel  formation.  This  hypothesis  can  partly  explain  why  these  mutations  are  neurotoxic 
simply  due  to  peptide  insertion  behavior.  We  further  apply  this  model  to  various  artificial  A/3  mutants  which  have  been  examined 
experimentally,  and  offer  testable  experimental  predictions  contrasting  the  roles  of  aggregation  and  insertion  with  regard  to 
toxicity  of  A/3  mutants.  These  can  be  used  through  further  experiments  to  test  our  hypothesis. 


INTRODUCTION 

Scientific  and  public  interest  in  Alzheimer’s  disease  has 
surged  in  the  last  several  decades.  The  reason  for  this  is 
simple:  with  increasing  life  expectancy,  Alzheimer’s  disease 
has  emerged  as  the  most  prevalent  form  of  late-life  mental 
failure  in  humans  (Selkoe,  2001). 

Alzheimer’s  disease  (AD)  is  a  neurodegenerative  disease 
involving  progressive  memory  impairment,  altered  behavior, 
decline  in  language  function,  disordered  cognitive  function, 
eventual  decline  in  motor  function,  and,  finally,  death 
(Selkoe,  2001).  In  AD,  the  brain  is  typically  marked  by  le¬ 
sions  (Selkoe,  2001),  neuronal  damage,  and  vascular  damage 
(Durell  et  al.,  1994).  These  lesions  are  typically  associated 
with  extracellular  plaques,  called  amyloid  plaques,  and 
intraneuronal  fibrillar  tangles  (Durell  et  ah,  1994;  Selkoe, 

2001) .  The  tangles  are  composed  of  a  protein  called  Tan  and 
are  called  Tau  tangles,  whereas  the  extracellular  plaques  are 
largely  composed  of  amyloid  /3  peptide  (A/3)  in  40-  and  42- 
residue  forms  (Selkoe,  2001)  (denoted  A/340  and  A/342, 
respectively).  These  insoluble  amyloid  plaques  composed  of 
A/3  are  considered  a  hallmark  of  AD.  However,  they  are  not 
specific  to  AD  (Dickson  and  Vickers,  2001)  and  have  been 
observed  in  older  patients  free  from  AD  symptoms  (Jarrett 
et  ah,  1993).  It  has  been  pointed  out  that  correlations 
between  amyloid  plaque  density  and  severity  of  dementia 
are  weak,  whereas  there  are  stronger  correlations  between 
soluble  A/3  levels  and  severity  of  dementia  (Walsh  et  ah, 

2002) .  This  is  one  reason  for  the  suggestion  that  oligomers  of 
A/3  may  be  more  important  to  toxicity  than  large  insoluble 
aggregates  or  plaques.  Evidence  for  this  idea  has  been 
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provided  in  vivo  (Walsh  et  ah,  2002)  and  in  vitro  (Hartley 
et  ah,  1999;  Lambert  et  ah,  1998). 

One  mechanism  by  which  oligomers  can  damage  cells  is 
formation  of  pores  or  ion  channels  through  the  cell 
membrane.  Early  work  in  this  area  showed  that  A/3  can 
insert  into  planar  lipid  bilayers  and  allow  a  calcium  current 
upon  insertion,  and  further  that  these  channels  can  be 
blocked  (Arispe  et  ah,  1993),  suggesting  that  the  calcium 
current  is  really  due  to  channel  formation,  not  just  bilayer 
permeabilization  by  the  peptide.  Theoretical  modeling  based 
on  predicted  secondary  structures  for  membrane-bound  A/3 
has  suggested  that  the  A/3  peptide  can  form  channels  with 
four  or  six  A/3  subunits  in  each  leaflet  of  the  bilayer  (for 
a  total  of  8  or  12  per  channel;  Durell  et  ah,  1994).  More 
recent  work  has  been  done  using  atomic  force  microscopy  to 
look  at  the  structure  of  A/3  inserted  in  planar  lipid  bilayers 
and  has  found  what  appear  to  be  channels  consisting  of  four 
or  six  visible  subunits  around  a  central  pore,  consistent  with 
the  theoretical  picture  described  above.  The  monomers 
oligomerize  after  insertion  into  the  bilayer.  Furthermore,  in 
the  presence  of  these  oligomers,  current  can  flow  (Lin  et  ah, 
2001).  Lin  et  ah  (2001)  also  show  that,  under  similar 
conditions,  A/342  induces  neuritic  degeneration  and  death  in 
cell  culture  and  that  this  toxicity  is  calcium-dependent  and 
blocked  by  zinc.  Imaging  work  by  another  group  has  also 
shown  that  A/340  oligomers  with  the  E22G  mutation  (where 
glutamate,  E,  at  residue  22  is  replaced  with  glycine,  G), 
which  causes  a  form  of  familial  AD,  can  form  pore-like 
structures  (Lashuel  et  ah,  2002).  These  pore-like  structures 
actually  could  be  intermediates  which,  when  not  membrane- 
bound,  build  up  into  the  amyloid  plaques  observed  in  the 
brain  of  AD  patients  (Lashuel  et  ah,  2002). 

Based  on  these  suggestions,  and  the  observation  of  Lin 
et  ah  (2001)  that  oligomers  in  the  membrane  form  after 
insertion  of  monomers,  we  model  insertion  of  the  A/3  peptide 
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into  the  cell  membrane.  We  first  examine  the  regular  A/340 
and  A/342  peptides,  then  the  40-  and  42-residue  versions  of 
all  of  the  mutations  in  the  A/3  peptide  that  are  known  to  cause 
familial  AD  (FAD)  and  reduce  the  average  age  of  onset  for 
the  disease  compared  to  people  with  sporadic  AD  (Selkoe, 
2001 ).  We  believe  FAD  mutants  provide  a  tool  for  assessing 
proposed  toxicity  mechanisms,  in  that  the  biological  toxicity 
mechanism  should  explain  why  these  mutants  cause  FAD. 
Our  reasoning  in  looking  at  these  mutants  is  that  if  the 
insertion  behavior  of  the  FAD  mutant  peptides  is  different, 
this  could  make  a  difference  in  the  prevalence  of  oligomers 
in  the  membrane  and  thus  have  an  effect  on  toxicity,  if 
membrane-associated  oligomers  are  indeed  important  for 
toxicity  in  vivo.  Although  some  earlier  modeling  work  has 
dealt  with  the  structure  of  A/340  in  a  lipid  bilayer  (Pellegrini- 
Calace  et  al.,  2003),  we  believe  this  work  is  the  first  to 
compare  insertion  of  FAD  mutants. 

This  system  is  modeled  using  a  Monte  Carlo  (MC)  code 
which  has  been  developed  to  study  insertion  behavior  of 
peptides  into  lipid  bilayers.  This  model,  which  is  specific  at 
the  amino  acid  level,  allows  us  to  simulate  larger  peptides 
and  longer  timescales  than  traditional  molecular  dynamics 
simulation  studies.  The  configurational  steps  are  sufficiently 
small  that  it  has  been  used  successfully  to  suggest  insertion 
mechanisms,  as  well  as  to  describe  insertion  conformations 
for  some  peptides  (Maddox  and  Longo,  2002a, b).  Here,  we 
find  that  in  all  cases  the  peptide  inserts  relatively  easily. 
However,  we  find  differences  in  the  conformations  the 
peptide  adopts  once  inserted.  These  differences  in  the 
prevalence  of  conformations  are  our  central  result.  Relative 
to  the  normal  A/3  peptide,  most  of  the  FAD  mutant  peptides 
are  more  likely  to  insert  only  partially  in  the  bilayer.  We 
point  out  similarities  between  this  partially  inserted  confor¬ 
mation  and  the  predicted  channel  structures  (Durell  et  al., 
1994).  Thus  we  suggest  that  FAD  mutants  may,  in  this  way, 
facilitate  formation  of  harmful  channels.  Moreover,  the 
A/342  peptide,  with  additional  hydrophobic  residues,  has 
a  greater  tendency  than  A/340  to  hang  up  in  this  confor¬ 
mation,  and  this  may  correlate  with  the  increased  toxicity 
of  A/342. 

INTRODUCTION  TO  THE  MODEL  AND  METHOD 
Model  energy  function 

The  Monte  Carlo  model  used  here  has  been  described  in  detail  in  an  earlier 
publication  in  this  journal  (Maddox  and  Longo,  2002a).  Accordingly,  we 
give  a  brief  overview  of  the  essentials  here  and  direct  the  interested  reader  to 
the  earlier  reference  for  greater  detail. 

The  model  follows  previous  work  from  the  past  decade,  most  notably  and 
closely  that  of  Milik  and  Skolnick  (1992,  1993)  and  Baumgaertner  (1996). 
Each  amino  acid  residue  is  treated  as  a  sphere  of  identical  1.5  A  radius. 
There  are  three  contributions  to  the  potential  energy  function  which  are 
residue-independent:  1),  Us,  which  is  a  hard-core  steric  interaction  prevent¬ 
ing  residue-residue  overlap;  2),  Ut,  an  energy  measuring  the  cost  of  rotating 
the  peptide  planes  of  successive  residues,  which  is  periodic  in  the  torsional 
angle  </>  between  successive  residues  and  has  a  shallow  minimum  at  (J>  — 


52.1°;  and  3),  f/A,  characterizing  the  energy  of  distortion  of  the  angle  6 
between  adjacent  bonds,  with  a  shallow  minimum  at  6  =  89.5°. 

The  lipid  bilayer  with  surrounding  water  is  modeled  as  a  medium  without 
molecular  specificity,  but  with  three  different  spatial  regions.  The  bilayer’s 
normal  is  taken  to  lie  along  the  z  axis,  so  these  regions  are  invariant  in  the 
x—y  plane.  The  bilayer  has  overall  thickness  2(zo  +  Zh),  where  zq  is  the  length 
of  the  acyl  chains  of  a  given  leaflet,  and  zh  is  the  width  of  the  headgroup 
region.  We  have  used  Zh  =  4.5  A  and  zq  —  13.5  A.  We  have  also  tried 
different  chain  lengths  z0  but  we  do  not  present  the  results  here  as  they  were 
not  significantly  different  except  for  reductions  in  the  amount  of  the 
transbilayer  conformation  of  the  inserted  peptide  when  the  membrane  is 
sufficiently  thick,  as  we  discuss  in  Results. 

The  water-lipid  medium  is  characterized  by  three  dimensionless 
functions,  two  of  which  couple  linearly  to  residue  specific  parameters  we 
discuss  in  the  next  paragraph.  These  functions  are  (as  shown  in  Fig.  1): 

1.  w(z),  which  measures  the  fractional  water  content;  this  is  modeled  as 
a  step  function  with  exponentially  rounded  edges  (decay  length  of  2  A) 
that  is  zero  in  the  hydrophobic  acyl  chains,  one  in  the  water  region,  and 
varies  smoothly  through  the  head  region. 

2.  A  polarity  function  p(z),  also  exponentially  rounded  with  the  same 
decay  length,  and  chosen  to  be  one  in  the  lipid  head  regions  and  water 
while  falling  to  some  small  value  l—fq  (where  /q,  the  polarity  factor, 
determines  the  polarity  of  the  tail  region,  with  larger /q  corresponding  to 
a  less  polar  tail  region)  after  approximately  one  residue  diameter  into  the 
tail  region. 

3.  A  hydrophobicity  function  y(z),  which  is  the  sum  of  two  exponentially 
rounded  step  functions:  one  which  is  zero  in  the  water  region  and 
saturates  in  the  head  region,  proportional  to  the  total  gain  of  hy¬ 
drophobic  energy  in  the  head  region,  and  a  second  which  saturates  in 
the  tail  region  after  approximately  one  residue  diameter  and  accounts  for 
the  hydrophobic  energy  gained  for  residues  penetrating  the  acyl  tail 
region.The  water  content  w(z)  couples  linearly  to  the  external  hydrogen 
bonding  energy  of  each  residue,  which  is  residue  independent  in  form. 
The  net  hydrogen  bonding  energy  is  taken  as  Uu  given  by  a  sum  over 
residues  i  as 

Un  =  X(w(z Off0  +  ( 1  -  w(zi))//int(/)) ,  (1) 


z-coordinate 


FIGURE  1  Functions  characterizing  bilayer  properties.  The  lipid  bilayer 
is  described  by  three  functions,  w(z)  for  the  fractional  water  content,  p{z )  for 
the  polarity,  and  y(z)  for  the  hydrophobicity.  Here  the  z  axis  is  perpendicular 
to  the  plane  of  the  bilayer,  and  the  functions  (and  the  bilayer)  are  symmetric 
at  ~z  =  0. 
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where  H0  —  —  6.12  kcal/mol  is  the  transfer  energy  of  an  unbonded  peptide 
group  to  water  and  Hmt(i)  the  internal  hydrogen  bonding  energy  associated 
with  a -helix  formation  given  by 

*.(0=3  I  VH(|rn  +  i  -  ri|).  (2) 

4  n=— 4,— 3,3,4 

VH  is  a  nearly  hard-core  function  of  the  separation  between  residue  i  and  the 
potential  helical  hydrogen  bonding  partners  along  the  peptide  chain,  as 
proposed  by  Milik  and  Skolnick  (1993). 

Residue  specificity  is  included  in  two  energies  associated  with  polarity 
and  hydrophobicity.  First,  a  potential  energy  term  Uq  —  Xi#o(*)p(zi)  *s 
included,  where  q0(i )  is  the  residue-specific  polar  energy  associated  with 
charged  or  partially  charged  functional  groups.  Second,  a  hydrophobic 
energy  Ub  =  is  included,  where 

B(i)  =  y(zi)b0(i)  +  (1  -  y(z ;))^ 

X  X  Wn  +  i-ril).  (3) 

n=— 4,— 3,3,4 

Here,  bo(i )  is  residue-specific  and  measures  the  water-to-alkane  Gibbs 
hydrophobic  transfer  energy  for  residue  i,  and  b\(i)  is  the  maximum 
reduction  in  hydrophobic  energy  due  to  helical  folding.  The  value  bo(i)  is 
taken  to  be  proportional  to  the  stochastic  accessible  area  of  the  residue.  The 
helical-folding  related  term  derives  from  the  loss  of  accessible  surface  area 
associated  with  helix  formation.  Values  for  qo(i),  ^o(0>  and  b\(i)  for  all 
residues  are  tabulated  in  the  previous  work  (Maddox  and  Longo,  2002a). 
Note  that  because  the  model  treats  the  lipids  and  water  only  as  media,  the 
hydrophobic  energy  must  be  included  explicitly  in  our  model  energy 
function. 

Because  of  the  way  the  hydrophobic  and  hydrogen  bonding  energies  are 
calculated — simply  based  on  local  helicity — the  simulation  is  biased  toward 
o;-helices,  as  /3-structure  involves  longer  range  interactions  and  is  not  taken 
into  account  by  the  model.  Therefore,  the  model  will  not  accurately  describe 
insertion  behavior  of  any  peptide  that  inserts  while  in  a  conformation  rich  in 
/3-structure.  Fortunately,  the  monomeric  A/3  peptide  is  predicted,  based  on 
secondary  structure,  to  be  a -helical  between  residues  15  and  40  or  42  (Durell 
et  al.,  1994;  Pellegrini-Calace  et  ah,  2003)  when  membrane-bound. 
Experimental  NMR  work  in  aqueous  sodium  dodecyl  sulfate  micelles, 
which  to  some  extent  resemble  a  water-membrane  medium,  confirms  this  for 
A/340  (Coles  et  al,  1998).  Thus  the  model’s  bias  away  from  /3-structure 
should  not  play  a  significant  role  here.  Indeed,  we  find  that  the  region 
mentioned  above  inserts  into  the  membrane  in  a  largely  helical  structure,  as 
described  below. 

Our  total  energy  is  then  taken  as  the  sum  U  =  Us  +  Ur  +  U\  +  Ub  + 
Uq  +  Un .  All  of  our  modeling  presented  below  is  done  at  pH  7.0  with 
temperature  305  K,  and  uses  a  polarity  factor  fq  —  0.85,  corresponding  to 
a  polarity  between  that  of  octanol  and  hexadecane.  The  choice  of  this  value 
is  based  on  experimental  studies  (Griffith  et  al.,  1974;  Roseman,  1988)  and 
earlier  simulation  work  (Maddox  and  Longo,  2002a). 


Monte  Carlo  simulation  details 

The  simulation  method  is  the  canonical  MC  method.  We  use  periodic 
boundary  conditions  in  all  three  directions,  and  in  the  case  where  the  peptide 
runs  across  the  boundary,  interactions  are  calculated  using  the  minimum 
separation  between  the  two  residues  in  question  (the  minimum  image 
convention).  New  peptide  conformations  are  generated  using  three  different 
sorts  of  moves: 


1.  Peptide  translation:  The  whole  peptide  is  randomly  translated  a  small 
distance  (between  0  and  0.2  A)  along  each  Cartesian  axis. 

2.  Spike  move  (two  sorts):  (a).  For  an  end  residue,  the  virtual  bond 
connecting  it  to  the  chain  is  rotated  slightly,  first  in  the  x— y  plane,  and 
then  in  the  y—z  plane.  The  angle  of  rotation  is  random,  between  0°  and 
20°.  (b).  For  a  central  residue,  the  residue  is  rotated  a  random  (between 
0°  and  20°)  amount  around  an  axis  joining  the  centers  of  its  nearest 
neighbors,  while  keeping  all  virtual  bond  lengths  fixed. 

3.  Slide  move:  A  random  virtual  bond  is  selected  and  all  residues  on  one 
side  (selected  randomly,  either  up  or  down  the  chain)  of  it  are  moved 
a  small,  random  amount  (between  0  and  0.2  A),  while  remaining  fixed 
relative  to  one  another.  The  move  leaves  the  initial  virtual  bond  the 
same  length  but  rotated  relative  to  the  residues  it  connects. 

One  MC  step  consists  of  one  modification  of  each  type  1, 2,  and  3,  where  the 
choice  of  residue  is  random  for  moves  2  and  3.  Modifications  are  accepted  or 
rejected  with  a  probability  given  by  the  usual  Boltzmann  factor  p  —  e~ rt  so 
that  favorable  moves,  with  a  negative  A U,  are  always  accepted,  and  some 
unfavorable  moves  are  accepted. 

In  our  work,  we  wanted  to  capture  insertion  behavior  without  biasing 
results  by  initial  peptide  conformations.  We  have  done  two  groups  of 
simulations  to  accomplish  this.  First,  we  have  started  the  peptide  outside  the 
bilayer  in  the  aqueous  phase  in  a  random  conformation.  Second,  we  have 
started  the  peptide  in  an  initially  helical,  fully  inserted  conformation. 

MC  simulations  can  be  used  to  investigate  nonequilibrium  properties 
(e.g.,  insertion  mechanisms)  or  equilibrium  properties  (e.g.,  inserted 
conformations)  of  a  system.  When  these  simulations  are  used  to  study 
equilibrium  properties,  it  is  important  to  ensure  that  the  system  has  fully 
equilibrated  before  data  collection  begins.  If  this  is  not  done  carefully,  one 
consequence  is  that  the  so-called  equilibrium  state  might  depend  on  the 
initial  conditions.  To  establish  an  appropriate  period  of  equilibration,  we 
monitored  the  average  energy  of  the  peptide  as  a  function  of  simulation 
steps.  As  the  peptide  equilibrates  (reaching  its  energetically  preferred 
conformation(s)),  the  average  energy  decreases  from  an  initially  higher 
value.  Thus  we  can  get  a  reasonable  idea  how  many  simulation  steps  it  takes 
for  this  to  happen  simply  by  plotting  energy  versus  step  number. 

Using  this  method,  we  find  that  for  insertion  from  an  initial  conformation 
outside  the  bilayer,  a  30-million-step  equilibration  period  is  usually 
sufficient,  whereas  for  an  initially  inserted  and  helical  conformation,  30 
million  steps  is  always  sufficient.  Although  both  initial  conformations 
eventually  produce  the  same  equilibrated  state,  the  inserted  helical  con¬ 
formation  converges  more  rapidly  and  is  used,  with  an  equilibration  period 
of  50  million  steps,  in  all  our  simulations  (unless  otherwise  noted). 

As  a  further  test  that  our  equilibration  period  is  sufficient,  we  have  also 
used  conformations  from  peptides  at  the  end  of  an  entire  simulation  run  as 
starting  points  for  new  simulations,  and  the  results  at  the  end  of  both 
simulations  are  within  our  error  bars  of  one  another. 

It  is  worth  pointing  out  that  no  equilibration  would  be  required  if  the 
insertion  mechanism  is  being  studied.  However,  in  this  work,  we  find  that  in 
every  case  insertion  is  fairly  easy,  as  we  discuss  in  Results.  Thus  we  focus  on 
peptide  conformations  at  equilibrium. 


Data  collection 

Every  MC  simulation  run  employs  a  unique  set  of  random  numbers, 
resulting  in  a  slightly  different  result  each  trial  (similar  to  the  way  in  which 
no  two  experimental  measurements  are  identical).  More  accurate  data  are 
therefore  generated  by  averaging  multiple  runs.  Here,  we  have  run 
a  minimum  of  10  trials  for  every  peptide:  five  beginning  in  initially  helical 
and  inserted  conformations,  and  five  with  initially  random  conformations 
outside  the  bilayer.  Following  equilibration  (discussed  above),  data  is 
collected  for  50  million  steps.  In  this  article  we  report  the  results  of  the 
initially  inserted  conformations,  thus  our  results  are  averaged  over 
a  minimum  of  five  such  trials.  However,  as  we  will  discuss  in  the  appendix 
on  data  analysis,  in  some  cases  we  use  more  trials. 
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RESULTS 

Overview 

Overall,  we  simply  input  the  sequence  of  A/340  and  A [342 
and  various  mutants  and  run  the  simulations.  As  discussed 
above,  we  do  multiple  trials  for  each  peptide  and  average  the 
results.  What  we  find,  briefly,  is  that  the  peptides  all  insert 
relatively  easily  into  the  bilayer  if  they  begin  initially  outside 
the  bilayer,  and  we  do  not  find  that  the  FAD  mutations 
significantly  effect  this  ease  of  insertion.  However,  we  find 
that  the  mutations  do  influence  the  conformation  the  peptide 
adopts  once  inserted  into  the  bilayer.  Therefore,  our  focus  in 
this  work  is  not  on  details  of  how  the  peptide  inserts,  though 
we  believe  it  is  likely  that  insertion  proceeds  via  one  of  the 
two  main  insertion  mechanisms  described  previously 
(Baumgaertner,  1996;  Maddox  and  Longo,  2002a). 

One  of  our  fundamental  results  is  that  the  peptide  appears  to 
exhibit  multiple  possible  inserted  conformations  which  have 
nearly  the  same  energies,  thus  allowing  the  peptide  to  switch 
between  conformations  often  in  the  course  of  a  simulation. 
We  have  previously  described  such  behavior  as  conforma¬ 
tional  partitioning  (Maddox  and  Longo,  2002b).  We  find  that 
that  the  A/3  peptide  and  its  mutants  can  always  adopt  the  same 
small  set  of  conformations.  However,  the  mutations  alter  the 
number  of  MC  steps  the  peptide  spends  in  each  of  these 
conformations  (which,  in  a  real  system,  would  correspond  to 
the  number  of  inserted  peptides  in  each  conformation).  These 
conformations  are  shown  in  Figs.  2  and  3. 

Essentially,  it  is  quite  easy  to  distinguish  three  conforma¬ 
tions:  the  first  has  the  last  several  residues  anchored  in  the 
lower  head  region.  We  call  this  conformation  transbilayer. 
The  second  conformation  is  similar,  but  usually  much  more 
prevalent  since  the  last  residues  are  hydrophobic  and  prefer 


FIGURE  2  Primary  inserted  conformations  of  the  A/3  peptide.  We  find 
that  in  every  case,  the  inserted  peptides  can  adopt  essentially  three  different 
conformations.  Mutations  appear  to  alter  the  percentage  of  steps  the  peptide 
spends  in  each  conformation  but  do  not  fundamentally  change  the 
conformations.  ( a )  Transbilayer.  The  peptide  inserts  with  the  last  several 
residues  near  the  C-terminus  in  the  lower  lipid  head  region;  the  portion 
crossing  the  bilayer  is  roughly  helical.  ( b )  Fully  inserted.  Just  like  a,  except 
the  last  several  residues  are  not  anchored  in  the  lower  head  region,  meaning 
that  the  conformation  is  fairly  flexible,  (c)  Partially  inserted.  Like  b,  except 
now  much  more  of  the  peptide  is  tethered  to  the  upper  head  region  by  the 
polar  residues  22-23  and  26-28,  whereas  before  only  residues  1-15  or  so 
were  in  the  upper  head  region.  The  conformations  shown  are  for  A/340,  but 
A/342  has  similar  conformations  with  two  additional  residues  (isoleucine 
and  alanine)  at  the  C-terminus. 


to  remain  in  the  tail  region.  This  is  especially  true  in  the  case 
of  A/342.  In  this  conformation,  the  essential  difference  is  that 
the  tail  of  the  peptide  is  not  anchored  and  thus  is  fairly  floppy 
and  able  to  change  the  angle  it  makes  with  the  z  axis  easily. 
We  call  this  conformation  fully  inserted.  The  third  con¬ 
formation  is  different  in  that  in  the  first  two,  only  residues 
1-15  or  so  remain  in  the  upper  lipid  head  region,  while  in  the 
third  conformation,  the  polar  residues  22-23  and/or  26-28 
also  remain  in  the  upper  head  region  (along  with  some  of 
their  neighbors),  where  there  is  still  some  water  content.  As 
a  result,  the  C-terminus  (residue  40  or  42)  does  not  stick 
down  into  the  lipid  tail  region  nearly  as  far  as  in  the  other  two 
conformations.  We  call  this  conformation  partially  inserted. 
In  this  case  all,  or  almost  all,  of  the  peptide  is  only  in  the 
upper  leaflet  of  the  bilayer. 

This  third,  partially  inserted  conformation  can  be  divided 
into  two  conformations  simply  by  distinguishing  whether  it 
is  residues  22-23  that  remain  in  the  upper  head  region,  or 
residues  26-28.  This  separation  of  conformations  is  shown 
in  Fig.  3.  Some  of  our  analysis  is  done  grouping  these 
together,  and  some  by  separating  them,  as  we  will  discuss 
below. 

Briefly,  we  find  that  when  the  FAD  mutations  have  an 
effect  on  the  insertion  behavior,  it  is  usually  by  causing  the 
mutant  peptides  to  favor  the  partially  inserted  conformation 
more  than  wild-type.  The  E22G  mutant  is  an  exception,  as  it 
essentially  eliminates  this  conformation.  We  refer  the  reader 
to  the  Appendix  for  detailed  discussion  of  our  data  analysis 
procedure.  Overall,  however,  the  basic  output  of  the 
simulation  is  the  number  of  steps,  or  percentage  of  steps, 
each  residue  in  the  peptide  spends  at  each  z-coordinate.  We 
can  plot  this  for  all  residues  (Fig.  4)  or  particular  residues 
(Figs.  5-7).  With  some  analysis  of  these  (the  data  analysis  is 
explained  in  the  Appendix)  we  are  able  to  get  accurate 
measurements  of  increases  or  decreases  in  the  number  of 
steps  the  peptide  spends  in  a  given  conformation,  relative  to 
wild-type.  We  are  able  to  do  this  whether  we  choose  to 
separate  the  peptide’s  conformations  into  three  or  four 
groups.  Results  obtained  using  these  methods  are  presented 
in  Table  1  and  Table  2. 


FIGURE  3  Subconfigurations  associated  with  locations  of  the  polar 
residues.  We  can  further  break  up  the  partially  inserted  conformation,  from 
Fig.  2  c,  into  two  conformations:  (a)  a  conformation  where  only  residues  22- 
23  remain  in  the  upper  head  region,  and  ( b )  a  conformation  where  both 
residues  22—23  and  26-28  remain  in  the  upper  head  region. 
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FIGURE  4  Number  of  steps  at  each  z-coordinate  (vertical  axis  is  z  axis; 
darker  means  more  steps)  plotted  versus  residue  number  (sequence  shown), 
from  1-40  for  A/340,  a ,  and  from  \-A2  for  A/342,  b.  Note  that  on  residue  30, 
there  appear  three  dark  regions,  corresponding  to  three  peaks,  whereas 
residues  26  and  28,  for  example,  have  four  peaks.  This  result  is  important  for 
our  data  analysis.  Note  also  that  the  transbilayer  conformation  for  A/342  is 
less  common  than  for  A/340  (compare  the  darkness,  or  number  of  steps,  of 
the  lowest  peak  on  residue  40  for  A/340  and  A/342).  The  sequence  of  A/3 
has  polar  residues  shaded,  hydrophobic  residues  underlined,  and  charged 
residues  with  charges  indicated. 


Normal  A/3  peptide  insertion 

As  detailed  in  the  Appendix,  some  results  for  A/340,  A/342, 
and  various  mutants,  are  shown  in  Figs.  5-7.  We  have  also 
extracted  the  number  of  steps  in  each  conformation  and 
present  these  (measured  as  percentages  of  total  steps)  in 
Table  1,  along  with  standard  deviations.  As  we  discuss  in  the 
Appendix,  the  percentages  for  conformations  ( b )  and  (c)  are 
not  completely  accurate  as  absolute  measures  of  the  number 
of  steps  in  those  conformations,  but  the  change  relative  to 
wild-type  is  accurate. 

It  is  important  to  note  that  there  is  a  fundamental 
difference  between  A/340  and  A/342.  While  both  can  adopt 
all  the  conformations  of  Fig.  2,  A/342  spends  many  fewer 
steps  in  conformation  (a).  This  is  because  the  last  two 
residues  of  A/342  add  significantly  to  the  hydrophobicity  of 
the  C-terminus  (see  sequence  in  Fig.  4),  as  both  are  nonpolar 
and  isoleucine  is  strongly  hydrophobic.  One  might  well  ask, 
however,  why  A/340  has  its  last  three  residues  in  the  lower 
head  region  at  all,  especially  the  final  two  valines,  which  are 
nonpolar  and  strongly  hydrophobic.  To  understand  this 
effect,  it  is  worth  noting  that  in  this  conformation,  the 
transbilayer  helix  begins  after,  or  around,  residue  16,  lysine, 
which  is  the  final  charged  and  polar  residue  before  a  number 
of  nonpolar,  hydrophobic  residues,  beginning  with  leucine. 
To  remove  residues  39  and  40  from  the  lower  head  region 
would  require  either  the  peptide  to  insert  at  a  more  shallow 
angle  so  that  these  residues  do  not  make  contact  with  the 
lower  head  region,  or  putting  a  kink  in  the  helix  (as  in  Fig.  2 
b)  so  that  they  do  not  make  contact.  In  the  first  case, 
changing  the  insertion  angle  would  move  lysine,  and 
possibly  residues  after  it  (it  is  immediately  followed  by 


/a)  A-Beta  40 

Residue  40 


(b)  A-Beta  42 

Residue  40 


FIGURE  5  Binned  numbers  of  steps  spent  at  each  z-coordinate  for  residue 
40  of  ( a )  A/340  and  ( b )  A/842.  Each  plot  also  shows  various  FAD  mutations. 
It  is  easily  apparent  that  A/342  spends  significantly  fewer  steps  in  the  fully 
inserted  conformation  ( leftmost  peak)  compared  to  A/340.  It  is  difficult  to  tell 
much  about  the  other  conformations  by  looking  at  this  distribution  for 
residue  40,  but  the  prevalence  of  these  can  be  extracted  from  other  residues. 
Note  that  the  lipid  head  regions  are  from  z  —  13.5  A  to  z  —  18  A  (and 
similarly  for  negative  z). 


valine),  into  the  upper  head  region  where  there  is  still  some 
water.  This  would  be  costly  energetically.  So  is  putting  a  kink 
into  the  helix.  Thus  we  believe  that,  in  the  case  of  A/340,  the 
energy  cost  of  either  putting  a  kink  in  the  helix,  or  imbedding 
K16  and  V17  in  the  upper  head  region,  is  comparable  to  the 
energy  cost  of  imbedding  V39  and  V40  in  the  lower  head 
region,  thus  the  conformation  in  Fig.  2  a  does  occasionally 
happen.  On  the  other  hand,  for  A/342,  the  additional  141 
makes  this  conformation  so  costly  that  it  almost  never 
happens.  This  reduction  in  the  transmembrane  conformation 
results  in  an  increase  in  the  prevalence  of  the  other  two 
conformations  relative  to  A/340  which,  if  correct,  could  help 
explain  why  A/342  is  typically  more  toxic. 

It  is  probably  important  to  point  out  that  our  results 
do  not  necessarily  mean  A/340  would  be  biologically 
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FIGURE  6  Binned  number  of  steps  spent  at  each  z-coordinate  for  (a) 
residue  30  of  A/340  and  ( b )  residue  35  of  A/342.  These  are  the  residues  we 
picked  which  best  distinguish  three  groups  of  conformations.  The  leftmost 
group  on  both,  which  is  very  small  for  A/342,  is  in  conformations  that  appear 
nearly  transbilayer.  The  middle  group  is  in  conformations  that  are  inserted 
and  fairly  floppy,  as  in  Fig.  2  b,  and  the  rightmost  group  is  in  the  partially 
inserted  conformation.  It  can  be  clearly  seen  that  for  A/340,  all  of  the 
mutations  but  E22G  and  A21G  result  in  an  increase  in  this  last  peak  relative 
to  wild-type,  and  for  A/342  all  of  them  except  E22G  do,  as  well. 


transmembrane.  We  have  tried  making  the  bilayer  slightly 
thicker  (by  5  A)  and  the  transbilayer  conformation  essen¬ 
tially  disappears.  This  means  the  transmembrane  insertion  is 
due  to  the  limited  space  between  head  regions,  as  we  argued 
above.  Biologically,  the  membrane  could  be  slightly  thicker 
than  our  model,  or  the  peptide  might  cause  a  small  bulge  in 
the  membrane  to  accommodate  a  fully  inserted,  rather  than 
transbilayer,  conformation.  So,  although  the  simulation  does 
produce  this  transbilayer  conformation,  we  have  no  reason 
to  believe  that  this  conformation  would  be  distinct  from 
the  fully  inserted  conformation  in  biological  systems.  On  the 
other  hand,  this  conformations  does  exist  in  our  model.  Here, 
it  is  the  reduction  in  the  prevalence  of  this  conformations  that 
makes  the  other  two  conformations  more  prevalent  for 
A/342  than  Aj340. 


(a) 


A-Beta  40 

Residue  26 


(b)  A-Beta  42 

Residue  26 


FIGURE  7  Binned  number  of  steps  spent  at  each  z-coordinate  for  residue 
26  of  (a)  A/340  and  ( b )  A/342.  For  this  residue,  there  are  four  apparent  peaks, 
corresponding  to  the  conformations  of  Fig.  2,  a  and  b,  and  Fig.  3.  For  A/342, 
the  transbilayer  conformation  is  so  small  that  there  is  no  apparent  peak.  For 
A/340,  it  is  the  leftmost  peak,  followed  by  the  fully  inserted  peak,  then  the 
partially  inserted  peaks:  the  peak  of  Fig.  3  a,  and  then  the  peak  of  Fig.  3  b. 
Notice  that  for  both  A/340  and  A/342,  the  FAD  mutants  increase  the  weight 
of  the  rightmost  conformations  (those  of  Fig.  3)  relative  to  wild-type,  except 
for  the  E22G  mutant. 


FAD  A/3  peptide  insertion 

FAD  mutations 

There  are  a  number  of  known  FAD  mutations,  including 
some  involving  Aft  (as  well  as  others  involved  in  other 
aspects  of  the  disease  including  Aft  production).  These  are 
named  by  the  populations  they  were  first  found  in  and  include 
Flemish  (A21G),  Arctic  (E22G),  and  Iowa  (D23N)  (Mur¬ 
akami  et  al.,  2002).  Murakami  et  al.  also  include  Dutch 
(E22Q)  and  Italian  (E22K)  but  there  is  some  dispute  about 
whether  these  are  properly  to  be  considered  AD  mutations 
(Melchor  et  al.,  2000;  Nilsberth  et  al.,  2001;  Wattendorff 
et  al.,  1995).  To  understand  this,  it  is  important  to  note  that 
AD  is  often  accompanied  by  cerebral  amyloid  angiopathy 
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TABLE  1  Frequency  of  each  conformation  for  native  and  FAD  A/3  peptides 

Three-peak  analysis  Four-peak  analysis 


Peptide  form  Onset  %  Trans.  %  Fully  ins.  %  Partially  ins.  A%  Upper(a)  A%  Upper(/?) 


WT  A/340 

72.8 

15.6 

-E 

3.8 

48 

+ 

5.2 

36.3 

-E 

6.8 

— 

- 

- 

— 

A/342 

2.1 

-E 

0.5 

58.5 

-E 

5.6 

39.4 

-E 

5.8 

- 

- 

- 

- 

A21G  A/340 

52 

17.6 

-E 

3.0 

48.3 

-E 

6.8 

34.1 

-E 

8.0 

0.5 

-E 

1.8 

2.9 

±  8.4 

A/342 

1.9 

-E 

0.7 

44.4 

-E 

5.4 

53.7 

-E 

5.4 

-2.8 

-E 

3.2 

14.8 

±  6.9 

E22G  A/340 

57 

27.4 

-E 

7.5 

67.4 

-E 

7.0 

5.13 

-E 

7.0 

-11.8 

-E 

1.6 

-13.2 

±  6.2 

A/342 

2.3 

-E 

1.7 

94.8 

-E 

7.0 

1.29 

-E 

1.0 

-16.8 

-E 

3.7 

-9.6 

±  5.6 

E22Q  A/340 

? 

13.6 

-E 

4.6 

43.6 

+ 

5.5 

42.8 

-E 

5.5 

6.4 

-E 

3.5 

2.8 

±  11.0 

A/342 

2.0 

-E 

0.3 

55.8 

+ 

6.9 

42.2 

-E 

7.0 

4.5 

-E 

3.4 

-1.8 

±  8.7 

E22K  A/340 

? 

15.5 

-E 

9.6 

37.9 

-E 

5.2 

46.6 

-E 

6.0 

9.9 

-E 

2.3 

2.2 

±  10.0 

A/342 

1.3 

-E 

0.8 

34.2 

+ 

10.5 

64.5 

-E 

9.4 

9.4 

-E 

8.0 

14.6 

±  17.1 

D22N  A/340 

69 

13.1 

-E 

9.2 

30.0 

-E 

8.8 

56.8 

-E 

10.2 

8.6 

-E 

7.2 

10.8 

±  16.2 

A/342 

1.2 

-1- 

0.5 

41.1 

-E 

16.3 

57.7 

-E 

16.7 

0.6 

-E 

5.8 

16.0 

±  11.4 

Ages  of  onset  (where  known)  (Grabowski  et  al.,  2001;  Li  et  al.,  2002;  Nilsberth  et  al.,  2001;  Roks  et  al.,  2000)  and  simulation  results:  percentage  transbilayer, 
percentages  in  the  fully  inserted  (Fig.  2  b)  and  partially  inserted  (Fig.  2  c)  conformations  as  calculated  from  the  three-peak  analysis  on  residue  30.  The  last 
two  columns  measure  the  change  in  percentage  of  the  conformations  of  Fig.  3,  a  and  b,  relative  to  wild-type  (from  the  four-peak  analysis).  For  each  mutation, 
one  row  represents  A/340  and  the  next  A/342.  Onset  age  is  not  specific  to  the  40-  or  42-residue  forms.  Here,  the  only  consistent  trend  that  we  find  is  that  most 
of  the  FAD  mutants,  with  the  exception  of  E22G,  appear  to  increase  the  percentage  of  the  partially  inserted  conformation  relative  to  wild-type.  This  appears 
to  be  true  for  both  A/340  and  A/342.  By  way  of  comparison,  the  FAD  mutations  increase  fibrillar  aggregation  of  A/342  in  vitro  in  every  case  except  the  A21G 
mutant  (Murakami  et  al.,  2002)  and  increase  soluble  A/3  levels  in  every  case  but  the  E22G  mutant  (Lashuel  et  al.,  2002;  Nilsberth  et  al.,  2001).  To  improve 
statistics,  wild-type  A/340  results  are  the  average  of  25  trials  rather  than  the  usual  5;  wt  A/342  are  the  average  of  20;  A21G  A/340  20  trials;  A21G  A/342  10 
trials,  and  A/342  E22Q  and  D23N  10  trials  each.  All  the  rest  are  five  trials,  as  described  in  Appendix:  Data  Analysis. 


(CAA),  deposition  of  A/3  in  blood  vessels  of  the  brain 
potentially  leading  to  vessel  rupture  and  stroke,  especially  in 
FAD  cases  (Murakami  et  al.,  2002).  Wattendorff  et  al.  (1995) 
point  out  that  classic  Alzheimer’s  plaques  are  rarely  found  in 
the  Dutch  CAA  case  and  dementia  and  death  are  due  to 
cerebral  hemorrhage  involving  damage  to  blood  vessels,  as  is 
also  the  case  in  the  Italian  E22K  mutant  (Melchor  et  al., 
2000).  But  the  fact  that  AD  also  involves  amyloid  angiopathy 
leaves  open  the  possibility  that  the  Dutch  and  Italian  forms 
are  vascular  forms  of  AD  (Wattendorff  et  al.,  1995). 

Here,  we  set  aside  the  issue  of  whether  or  not  the  Dutch 
and  Italian  A/3  mutations  are  actually  AD  mutations  or 
whether  they  should  be  regarded  as  something  different  and 
simply  model  insertion  of  these  peptides  into  cell  mem¬ 
branes.  It  is  known  that  even  the  Dutch  E22Q  and  Italian 
E22K  mutant  peptides  interact  with  cell  surfaces  of  cerebral 
vascular  smooth  muscle  cells  and  cause  cell  death  in  vitro 
(Melchor  et  al.,  2000).  Thus,  it  is  conceivable  that  the 
mechanism  may  be  similar  to  that  described  by  Lin  et  al. 
(2001)  and  the  type  of  cells  being  damaged  may  simply  be 
different — either  cerebral  vascular  smooth  muscle  cells 
(Melchor  et  al.,  2000)  or  brain  pericytes  (Verbeek  et  al., 
1997)  rather  than  neurons.  For  simplicity,  we  will  call  all  of 
these  FAD  mutations. 

It  is  also  important  to  note  that  these  A/3  mutations  are 
autosomal-dominant  (Nilsberth  et  al.,  2001)  and,  in  addition 
to  lowering  the  age  of  onset  for  AD  compared  to  sporadic 
cases,  cause  AD  in  all  subjects  with  the  mutations  who  live 
long  enough.  Many  of  these  mutations  also  lead  to  increased 


A/3  levels,  but  that  cannot  be  the  sole  cause,  as  at  least  the 
Arctic  mutation  (Lashuel  et  al.,  2002)  leads  to  decreased 
levels.  The  cause  of  disease  also  cannot  be  simply  the 
increased  propensity  to  form  fibrils  or  aggregates,  because  at 
least  the  Flemish  mutation  does  not  increase  fibril  formation 
(Murakami  et  al.,  2002). 

FAD  mutation  results 

Results  for  these  mutations  are  shown  in  Figs.  5-7,  and 
presented  in  reduced  fomi  in  Table  1.  We  find  that  the  most 
consistent  difference  in  results  relative  to  wild-type,  for  both 
A/340  and  A/342,  is  that  four  of  the  five  FAD  mutants 
increase  the  number  of  steps  the  peptide  is  in  a  conformation 
like  that  of  Fig.  2  c.  The  E22G  mutant,  however,  results  in 
a  huge  decrease  of  this  conformation. 

To  understand  our  results,  consider  the  changes  in  polarity 
and  hydrophobicity  the  mutations  involve.  The  polarity  and 
hydrophobicity  values  used  in  this  model  have  been  tabulated 
before  (Maddox  and  Longo,  2002a)  and  will  not  be  repeated 
here,  but  it  is  worthwhile  to  discuss  the  mutations  briefly.  The 
A21G  mutant  involves  a  decrease  in  hydrophobicity,  hence 
should  cause  residue  21  to  more  strongly  prefer  to  be  in  the 
upper  head  region  and  thus  result  in  an  increase  in  the 
prevalence  of  the  partially  inserted  conformation  (Fig.  2  c). 
This  is  what  we  find,  at  least  for  A/342  (for  A/340,  the 
percentage  stays  about  the  same).  The  E22Q  mutant  involves 
a  fairly  large  increase  in  polarity  and  slight  increase  in 
hydrophobicity,  thus  the  polarity  should  dominate  and  cause 
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TABLE  2  Frequency  of  each  conformation  for  artificial  A/3  mutants 

Three-peak  analysis  Four-peak  analysis 


Peptide  form  %  Trans.  %  Fully  ins.  %  Partially  ins.  A%  Upper(tf)  A%  Upper(/?) 


WT  A/340 

15.6 

-t- 

3.8 

48 

+ 

5.2 

36.3 

-F 

6.8 

- 

— 

- 

— 

A/342 

2.1 

+ 

0.5 

58.5 

+ 

5.6 

39.4 

-F 

5.8 

- 

- 

- 

- 

E22QD23N  A/340 

15.5 

-t- 

5.0 

30.0 

+ 

6.4 

54.6 

-F 

6.9 

4.4 

-F 

6.3 

13.9 

±  16.2 

A/342 

0.9 

+ 

0.7 

24.6 

+ 

11.4 

74.4 

-F 

11.4 

4.0 

-F 

10.5 

30.9 

±  19.0 

E22A  A/340 

23.1 

-t- 

7.6 

71.7 

-F 

7.1 

5.2 

-F 

0.7 

-10.8 

-F 

2.3 

-13.3 

±  6.3 

A/342 

2.0 

-t- 

0.4 

85.6 

+ 

2.1 

12.4 

-F 

2.0 

-13.8 

-F 

7.0 

-11.4 

±  4.2 

E22D  A/340 

22.9 

-t- 

8.9 

64.9 

-F 

7.6 

12.2 

-F 

2.0 

-3.9 

-F 

2.6 

-13.1 

±  6.2 

A/342 

1.9 

+ 

1.1 

62.9 

-F 

6.6 

35.2 

-F 

6.6 

-5.9 

-F 

3.4 

-0.4 

±  10.7 

A2S  A/340 

17.6 

-F 

9.1 

49.3 

-F 

6.4 

33.1 

-F 

5.0 

-1.2 

-F 

9.9 

-3.2 

±  9.7 

A/342 

2.2 

-F 

1.3 

60.4 

-F 

3.4 

37.4 

Hr 

3.3 

2.4 

-F 

3.3 

-5.3 

±  8.7 

F19S  A/340 

1.6 

-F 

2.4 

5.3 

-F 

4.4 

93.1 

-F 

4.9 

4.9 

-F 

10.3 

57.9 

±  16.1 

A/342 

0.1 

-F 

0.1 

10.8 

-F 

5.2 

89.1 

-F 

5.2 

7.8 

Hr 

13.6 

43.6 

±  14.9 

I32S  A/340 

14.2 

-F 

7.0 

34.7 

-F 

8.7 

51.1 

-F 

10.7 

-1.7 

-F 

3.2 

17.0 

±  12.6 

A/342 

1.5 

+ 

0.7 

49.5 

-F 

13.9 

49.1 

Hr 

14.1 

-1.1 

-F 

4.5 

10.2 

±  11.5 

13  2 V  A/340 

16.0 

-t- 

3.9 

58.8 

-F 

2.8 

25.3 

HF 

1.3 

4.4 

-F 

1.6 

-10.8 

±  6.8 

A/342 

2.5 

+ 

1.2 

63.2 

-F 

2.1 

34.4 

-F 

2.0 

-0.1 

-F 

3.7 

-5.7 

±  5.6 

V36E  A/340 

83.0 

-t- 

13.9 

10.6 

-F 

9.3 

6.4 

-F 

6.2 

-12.7 

-F 

3.0 

-5.3 

±  13.3 

A/342 

16.7 

-t- 

11.9 

36.2 

-F 

19.5 

47.1 

-F 

10.5 

1.6 

-F 

5.1 

5.6 

±  11.8 

H6R  A/340 

17.3 

-t- 

4.8 

47.7 

-F 

7.7 

35.0 

-F 

8.4 

2.2 

-F 

2.5 

-1.2 

±  8.8 

A/342 

2.0 

-t- 

0.6 

52.9 

-F 

8.2 

45.1 

-F 

8.4 

1.8 

-F 

3.7 

2.5 

±  7.3 

Shown  are  results  for  artificial  A/3  mutants  mentioned  in  the  literature:  percentage  transbilayer,  percentages  in  the  fully  inserted  (Fig.  2  b)  and  partially  inserted 
(Fig.  2  c)  conformations  as  calculated  from  the  three-peak  analysis  on  residue  30.  The  last  two  columns  measure  the  change  in  percentage  of  the  conformations  of 
Fig.  3,  a  and  b,  relative  to  wild-type  (from  the  four-peak  analysis).  The  A/340  E22QD23N  mutant  is  more  toxic  to  HCSM  cells  than  either  mutant  alone  (Van 
Nostrand  et  al.,  2001),  whereas  the  E22D  mutant  of  A/340  is  not  toxic  to  HCSM  cells  and  the  E22A  mutant  is  (Melchor  et  al.,  2000).  The  A2S,  F19S,  I32V,  I32S, 
and  V36E  mutants  are  known  to  reduce  aggregation  of  A/342  (Wurth  et  al.,  2002),  and  the  H6R  mutant  has  been  suggested  as  an  FAD  mutant  (Janssen  et  al.,  2003). 


residue  22  to  strongly  prefer  the  upper  head  region,  as  we  find. 
The  E22K  mutant  is  quite  similar,  but  involves  an  even  larger 
increase  in  polarity,  and  should  also  cause  residue  22  to  prefer 
the  upper  head  region,  as  we  find.  And  the  D23N  mutant 
involves  a  similar  increase  in  polarity  and  slight  increase  in 
hydrophobicity,  so  again  the  polarity  should  dominate  and 
increase  this  conformation,  as  we  find.  On  the  other  hand,  the 
E22G  mutant  replaces  a  polar,  charged  residue  with  a  neutral, 
nonpolar  residue  and  thus  residue  22  is  now  mostly  indif¬ 
ferent.  This  results  in  the  observed  huge  decrease  in  the  par¬ 
tially  inserted  conformation. 

On  one  hand,  the  different  behavior  of  the  E22G  mutant 
could  be  taken  as  evidence  that  there  is  no  consistent  change 
in  insertion  behavior,  and  thus  suggest  that  this  may  not  be 
the  toxicity  mechanism,  or  at  least  that  oligomerization,  not 
insertion,  is  important.  That  is  certainly  a  possible  meaning 
of  these  results.  However,  we  would  suggest  that  the  reader 
compare  the  conformation  of  Fig.  2  c  with  the  channel 
structure  of  Durell  et  al.  (1994),  specifically  in  Fig.  4  of  that 
article.  Note  that  there  are  two  helices  predicted,  from 
residues  15-24,  running  near  the  upper  head  region,  and 
from  25^10  or  25-42,  inserted.  In  Fig.  2  c,  we  have  a  nearly 
helical  region  from  around  residues  14-24  or  25  that  is  near 
the  upper  head  region,  and  one  from  25-40  which  hangs 


down  inserted.  Thus  it  resembles  the  predicted  structure 
more  than  either  of  the  other  conformations.  In  the  structure 
prediction,  the  part  of  the  chain  before  residue  15  folds  down 
inside  the  oligomer  and  forms  a  /3-barrel,  which  we  obviously 
do  not  capture  here  because  we  do  not  have  oligomers  neces¬ 
sary  to  stabilize  such  a  structure,  and  because  the  model  is 
not  designed  to  capture  /3-structure.  With  the  similarity  of 
Fig.  2  c  to  the  predicted  channel  structure  in  mind,  we  form 
the  hypothesis  that  the  partially  inserted  conformations 
(Fig.  2  c,  or  Fig.  3,  a  and  b)  are  more  likely  to  form  channels, 
simply  because  of  their  resemblance  to  the  structure  of  the 
monomers  making  up  the  channels.  Below,  in  the  subsection 
“Insertion  of  other  mutant  A/3  peptides”,  we  test  this 
hypothesis  against  some  experimental  data  for  artificial  A/3 
mutants,  and  present  predictions  for  further  tests.  One  might 
argue  that  the  observed  changes  in  the  prevalence  of  different 
conformations  are  relatively  small  and  thus  would  be 
unlikely  to  result  in  a  large  difference  in  the  toxicity  of  the 
different  mutants.  However,  these  channels  consist  of  eight 
or  twelve  monomers.  If  one  thinks  of  a  large  number  of 
monomers  inserted  in  a  bilayer,  some  will  be  in  each 
conformation  at  any  given  time,  so  we  can  think  of  the 
concentration  of  each  conformation.  Consider,  then,  the 
concentration  of  the  conformation  that  can  form  channels. 
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The  time  it  takes  for  12  monomers  within  a  bilayer  to  find 
each  other,  or  the  probability,  should  scale  as  c12,  where  c  is 
this  concentration.  Thus  a  small  difference  in  the  amount  in 
the  proper  conformation  could  make  a  big  difference  in  the 
likelihood  of  forming  channels. 

Insertion  of  other  mutant  A p  peptides 

A  variety  of  other  data  is  available  on  mutant  A/3  peptides. 
Some  reduce  aggregation  in  vitro  (Wurth  et  al.,  2002),  some 
might  cause  FAD  (Janssen  et  ah,  2003),  and  some  have 
various  effects  on  cultured  cerebrovascular  smooth  muscle 
cells  (Melchor  et  ah,  2000;  Van  Nostrand  et  ah,  2001).  Here, 
we  examine  some  of  the  mutations  which  have  been  men¬ 
tioned  in  the  literature  as  a  test  of  our  hypothesis — if  it  is 
right,  and  toxicity  depends  in  part  on  the  relative  prevalence 
of  the  partially  inserted  conformation,  our  results  ought  to 
correlate  with  the  experimental  toxicity  measurements. 

First,  we  address  the  E22Q.D23N  A/340  double  mutant 
created  by  Van  Nostrand  et  al.  (2001).  This  mutation’s  effect 
has  been  examined  on  human  cerebrovascular  smooth  muscle 
cells  (HCSM  cells)  because  the  E22Q  mutation  causes 
especially  pronounced  cerebral  amyloid  angiopathy  and 
patients  with  this  mutation  typically  die  of  hemorrhage,  as 
discussed  above.  HCSM  cells  are  known  to  degenerate  in 
CAA  in  a  manner  that  is  associated  with  A/3  deposition 
(Melchor  et  al.,  2000).  Van  Nostrand  et  al.  (2001)  found  that 
the  E22Q.D23N  double  mutant  is  even  more  toxic  to  HCSM 
cells  than  E22Q  or  D23N  alone.  Our  idea  was  that  the 
mechanism  for  this  toxicity  also  involves  insertion  of  the 
peptide  and  formation  of  channels,  so  we  modeled  this 
mutant  as  well.  Results  for  this  mutant  are  shown  in  Table  2; 
we  find  that  the  prevalence  of  the  partially  inserted  con¬ 
formation  does  increase  relative  to  wild-type.  Relative  to 
the  D23N  mutant,  which  we  would  predict  would  be  the  more 
toxic  of  the  E22Q  and  D23N  mutants,  we  observe  an  increase 
in  the  partially  inserted  conformation  only  for  A/342.  Again, 
however,  we  would  argue  simply  based  on  the  residue 
properties  that  since  both  the  E22Q  and  D23N  mutants 
increase  preference  of  those  residues  for  the  upper  head 
region,  the  double  mutant  should  have  a  stronger  effect  on 
this  than  either  alone,  and  thus  we  expect  that,  given  better 
statistics,  we  would  agree  with  Van  Nostrand  et  al.  (2001)  and 
predict  that  the  double  mutant  is  more  toxic  than  either  alone. 

Melchor  et  al.  (2000)  have  found  that  an  artificial  E22D 
mutant  of  A/340  does  not  effect  HCSM  cells,  in  contrast  to 
biological  E22Q  and  E22K  mutants.  They  also  observed  that 
the  A/340  E22A  mutant  is  toxic  to  HCSM  cells.  Therefore  we 
model  insertion  of  E22D  and  find  (Table  2)  that  the  E22D 
mutant  results  in  a  decrease  in  the  partially  inserted  con¬ 
formation  (as  one  would  expect  due  to  the  decrease  in  polar¬ 
ity),  especially  for  A/340.  The  E22A  mutant  results  in  a  large 
decrease  in  the  partially  inserted  conformation  for  both 
A/340  and  A/342  (again,  as  one  would  expect  due  to  the  de¬ 
crease  in  polarity).  Thus,  based  on  our  hypothesis,  we  would 


agree  that  the  E22D  A/340  mutant  would  not  be  toxic,  but 
disagree  that  the  E22A  mutant  would  be  toxic.  This  could  be 
taken  as  evidence  against  our  hypothesis  or  evidence  that  the 
toxicity  mechanism  is  different  for  HCSM  cells. 

Janssen  et  al.  (2003)  recently  identified  a  previously 
unknown  mutation  in  the  A/3  peptide  in  two  early-onset  AD 
patients  in  the  same  family.  This  mutation,  H6R,  produced 
ages  of  onset  ~55.  We  have  here  tried  this  mutation  (results 
in  Table  2)  and  find  that  it  produces  insertion  behavior  that  is 
within  error  bars  of  wild-type.  Given  the  position  of  the 
mutation,  this  is  what  we  would  expect,  as  it  is  within  the 
range  of  amino  acids  (1-14  or  more)  that  are  firmly  anchored 
in  or  near  the  upper  head  region,  where  there  is  some  water 
content.  Replacing  histidine  with  arginine,  which  is  even 
more  polar,  does  not  have  a  strong  effect  on  this  as  both  try  to 
remain  where  there  is  water,  i.e.,  the  surface  of  the  upper 
head  region.  Thus,  if  toxicity  depends  only  on  insertion 
conformation,  we  would  suggest  that  this  is  indeed  not  a  FAD 
mutation.  Thus  we  would  suggest  testing  toxicity  of  this 
mutant  in  cell  culture,  particularly  as  already  done  for  A/342 
by  Lin  et  al.  (2001). 

Some  in  vitro  work  has  been  done  to  find  artificial  mutants 
that  can  reduce  aggregation  of  A/342  (Wurth  et  al.,  2002). 
Unfortunately,  this  has  not  yet  been  extended  to  include 
A/340.  However,  we  selected  some  of  the  point  mutations 
which  are  known  to  reduce  aggregation  of  A/342  and 
modeled  the  insertion  of  these.  We  tried  A2S,  F19S,  I32V, 
I32S,  and  V36E.  Results  are  shown  in  Table  2.  We  find  that 
F19S  and  I32S  strongly  increase  the  prevalence  of  the 
partially  inserted  conformation  for  both  A/340  and  A/342. 
For  F19S,  this  is  due  to  the  substitution  of  polar  serine  for 
strongly  hydrophobic  phenylalanine,  causing  residue  29  to 
prefer  the  upper  head  region;  for  I32S,  the  reason  is  similar. 
In  this  case,  the  conformation  is  actually  different  from 
normal  in  that  residue  32  also  sticks  in  the  upper  head  region. 
In  contrast,  the  I32V  mutant  does  not  result  in  a  large  change 
relative  to  wild-type,  consistent  with  the  relatively  small 
change  between  isoleucine  and  valine.  Unsurprisingly,  the 
A2S  mutant  makes  no  change  to  insertion  behavior  (residue 
2  is  firmly  in  the  upper  head  region,  so  changing  it  to  polar 
makes  little  difference).  The  V36E  mutant,  however,  as  one 
might  expect,  drastically  increases  the  number  of  steps  that 
the  peptide  is  transbilayer,  as  the  polar  and  charged  glutamic 
acid  strongly  prefers  to  be  in  an  environment  with  more 
water.  Thus  it  decreases  the  prevalence  of  the  partially 
inserted  and  fully  inserted  conformations.  Thus  we  would 
predict  that,  if  our  hypothesis  is  correct,  F19S  and  I32S 
should  be  the  most  toxic,  I32V  and  A2S  should  be  com¬ 
parable  to  wild-type,  and  V36E  might  be  much  less  toxic 
than  wild-type.  We  say  might  because  this  probably  depends 
on  the  thickness  of  the  bilayer;  if  the  bilayer  is  too  thick  it 
will  probably  behave  just  like  wild-type  since  the  E36  could 
not  reach  the  lower  head  region.  Thus  a  testable  experimental 
prediction  of  our  hypothesis  is  that  toxicity  of  these  mutants 
would  be  related  to  their  insertion  behavior  as  just  described. 
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We  can  make  a  second  prediction  which  is  simply  based 
on  the  observed  insertion  behavior  of  A/342.  Looking  at  the 
insertion  behavior  of  the  A/342  versions  of  the  Wurth  et  al. 
(2000)  mutants,  we  find  that  I32V,  A2S,  and  possibly  V36E 
insert  more  like  natural  A/342,  whereas  F19S,  132S,  and 
possibly  V36E  insert  differently.  The  reduction  in  aggre¬ 
gation  splits  the  group  differently — I32S  aggregates  most, 
then  A2S,  V36E,  and  F19S  are  similar  to  one  another 
and  intermediate,  and  I32V  aggregates  least  (Wurth  et  al., 
2002).  So  if  fibrillar  aggregation  primarily  causes  toxicity, 
experiments  looking  at  toxicity  should  see  the  latter  group¬ 
ing,  whereas  if  insertion  behavior  is  of  much  more  impor¬ 
tance,  toxicity  experiments  should  see  the  former  grouping. 

DISCUSSION  AND  CONCLUSIONS 

We  have  here  presented  work  applying  a  model  of  peptide 
insertion  to  A/3,  a  peptide  implicated  in  Alzheimer’s  disease. 
Specifically,  we  have  examined  the  effect  of  FAD  mutations 
on  the  peptide’s  insertion  behavior,  with  the  idea  that  any 
successful  hypothetical  toxicity  mechanism  should  be  able  to 
explain  why  FAD  mutations  are  toxic.  Thus  if  FAD  mutants 
do  not  affect  peptide  insertion  into,  or  oligomerization 
within,  membranes,  the  ion  channel  toxicity  mechanism 
proposed  previously  (Arispe  et  al.,  1993;  Durell  et  al.,  1994; 
Lin  et  al.,  2001)  is  probably  not  relevant  biologically. 

What  we  find  is  that  the  FAD  mutations  do  affect  peptide 
insertion.  Four  of  these  five  mutations  involve  an  increase  in 
polarity  or  decrease  in  hydrophobicity  and  thus  cause  the 
peptide  to  prefer  (relative  to  wild-type)  a  conformation 
where  those  residues  are  in  the  upper  lipid  head  region  (i.e.. 
Fig.  2  c).  It  is  interesting  to  note  that  a  channel  stmcture 
suggested  previously  (Durell  et  al.  (1994),  see  Fig.  4)  has 
these  residues  laying  along  the  surface  of  the  bilayer.  Thus 
we  find  that  four  of  the  five  FAD  mutations  increase  the 
resemblance  to  this  configuration. 

Based  on  this  similarity,  we  develop  the  hypothesis  that 
causing  the  peptide  to  hang  up  in  the  upper  leaflet  (the 
partially  inserted  conformation)  facilitates  formation  of 
harmful  channels.  We  test  this  hypothesis  on  several  artificial 
mutations  examined  in  vitro  and  find  that  it  can  explain  the 
change  in  toxicity  of  two  and  is  wrong  on  a  third.  As  a  fur¬ 
ther  test  of  our  hypothesis,  we  can  offer  some  testable 
predictions.  For  example,  if  the  hypothesis  is  right,  we  would 
suggest  that  the  F19S  mutant  would  increase  toxicity  of  A/3 
relative  to  wild-type.  Additionally,  there  are  several 
mutations  known  to  reduce  aggregation  of  A/342  that  we 
predict  would  promote  the  insertion  behavior  that  would 
facilitate  channel  formation.  Thus,  it  would  be  simple  to 
distinguish  between  the  channel  formation  toxicity  mecha¬ 
nism  and  the  aggregation  toxicity  mechanism  by  looking  at 
the  toxicity  of  these  mutants.  Even  more  evidence  could  be 
provided  by  replicating  the  work  of  Lin  et  al.  (2001)  but 
using  various  FAD  mutants  and  looking  at  how  these  effect 
the  abundance  of  channels.  Additional  information  could 


also  be  gained  from  theoretical  work  along  the  lines  of  that 
by  Durell  et  al.  ( 1994)  to  see  what  effect  these  FAD  mutants 
would  have  on  channel  structures. 

Overall,  this  approach  of  modeling  peptide  insertion 
provides  a  simple  way  of  making  concrete  predictions  to 
distinguish  the  proposed  mechanism  of  channel  formation 
from  others.  The  limitation  of  this  approach,  however,  is  that 
we  can  only  look  at  insertion  of  single  peptides,  and  not 
interaction  between  these.  This  is  tine  if  formation  of  channels 
depends  on  having  peptides  initially  in  the  correct  conforma¬ 
tions,  which  seems  reasonable.  However,  it  is  also  possible 
that  interaction  between  inserted  peptides  causes  them  to 
adopt  conformations  appropriate  for  channel  formation.  If  this 
is  the  case,  then  insight  into  channel  formation  would  require 
a  more  sophisticated  model  that  can  include  interaction 
between  peptides.  Even  if  this  is  not  the  case,  including 
interaction  between  peptides  will  certainly  answer  the  ques¬ 
tion  of  whether  FAD  mutants  affect  channel  formation  in 
a  consistent  way  much  more  thoroughly  than  we  are  able  to  do 
here.  Therefore,  this  would  be  a  logical  continuation  of  this 
work  and  is  something  we  hope  to  do  in  the  future. 

Even  if  our  hypothesis  proves  to  be  wrong,  we  have 
shown  that  the  FAD  mutants  affect  insertion  behavior  of  the 
A/3  peptide  into  lipid  bilayers,  and  provided  an  understand¬ 
ing  as  to  why  the  FAD  mutants  affect  insertion  in  the  way 
they  do. 

In  conclusion,  our  work  shows  that  most  FAD  mutations 
have  a  significant  effect  on  the  insertion  of  the  A/3  peptide  in 
lipid  bilayers  in  this  theoretical  model,  and  this  effect  can 
easily  be  understood  by  looking  at  the  change  in  polarity  and 
hydrophobicity  accompanying  the  mutations.  The  effect  of 
FAD  mutations  on  insertion  has  not  been  studied  previously, 
and  may  be  significant.  Additionally,  we  offer  a  hypothesis 
based  on  promoting  channel  formation  by  causing  peptides 
to  insert  less  fully  that  can  help  explain  toxicity  of  A/3  FAD 
mutants,  as  well  as  several  artificial  mutants  studied  in  vitro. 
While  this  hypothesis  is  unproven,  it  is  based  on  the 
observation  that  these  peptides  do  insert  into  cell  membranes 
and  form  ion  channels  (Lin  et  al.,  2001),  and  similarities  to 
theoretically  predicted  channel  structures  (Durell  et  al., 
1994).  We  provide  testable  predictions  based  on  this 
hypothesis.  It  should  be  simple  for  experimentalists  to 
disprove  this  hypothesis,  if  it  is  false,  or  to  offer  additional 
evidence  for  it,  if  they  follow  the  experimental  suggestions 
we  offer  above.  Additionally,  our  work  suggests  the  value  of 
further  modeling  work  to  describe  the  full  formation  of  these 
channels,  rather  than  just  single-peptide  insertion. 

APPENDIX:  DATA  ANALYSIS 

The  basic  output  of  the  simulation  is  the  binned  number  of  steps — essen¬ 
tially  the  frequency  with  which  each  residue  is  found  at  each  z-coordinate. 
This  can  be  plotted  across  all  residues  and  illustrates,  as  in  Fig.  4,  that  there 
are  multiple  conformations  that  the  peptide  can  switch  between. 

In  terms  of  data  analysis,  it  is  more  useful  to  plot  z-distributions  of 
specific  residues  which  can  be  used  to  distinguish  between  the  conforma- 
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tions  described  above.  For  example,  residue  40  has  a  well-defined  peak  in 
the  lower  head  region  that  can  be  used  to  distinguish  the  transbilayer  con¬ 
formation  from  the  other  conformations.  A  plot  of  residue  40  z-distributions 
for  various  mutations  is  shown  in  Fig.  5. 

It  is  somewhat  more  difficult  to  distinguish  exactly  how  many  steps  the 
peptide  spends  in  the  fully  inserted  and  partially  inserted  conformations  (Fig. 
2,  b  and  c).  To  understand  this,  it  is  important  to  recognize  that  each 
conformation  results  in  a  peak  for  a  given  residue.  That  is,  from  Fig.  2,  one 
can  easily  see  that  the  position  of  residue  28  has  three  significantly  different 
locations  depending  on  which  conformation  the  peptide  is  in.  As  it  turns  out, 
it  also  has  a  fourth,  which  is  not  very  different  from  the  third,  as  can  be  seen 
in  Fig.  3.  Unfortunately,  these  peaks  tend  to  overlap  with  each  other  quite 
significantly — that  is,  residue  28  can  sometimes  have  similar  locations 
whether  it  is  in  the  partially  inserted  conformation  or  the  fully  inserted 
conformation.  This  means  that,  to  extract  the  number  of  steps  in  a  given 
conformation  from  the  data,  it  becomes  necessary  to  fit  some  function  to  the 
peaks  and  then  calculate  the  number  of  steps  from  that. 

To  do  this,  we  select  residues  where  the  peaks  appear  to  be  particularly 
well  separated.  For  one,  we  choose  residue  30  (for  A/340),  because  at  that 
residue,  the  peaks  corresponding  to  the  two  partially  inserted  conformations 
actually  overlap — that  is,  residue  30  has  the  same  average  location  no  matter 
whether  it  is  residues  22-23  alone  that  are  stuck  in  the  upper  head  region, 
or  22—23  and  26-28.  This  means  that  we  can  extract  the  weight  of 
conformation  from  Fig.  2  c — the  conformation  grouping  these  two 
together — by  looking  at  this  residue.  For  A/342,  the  two  peaks  of  the 
partially  inserted  conformation  do  not  quite  overlap  for  residue  30,  but  do  for 
residue  35,  so  we  do  the  same  analysis,  but  for  residue  35.  A  plot  of  the 
number  of  steps  at  each  z-value  for  these  residues  is  shown  in  Fig.  6. 

To  separate  the  two  partially  inserted  conformations  of  Fig.  3,  we  choose 
three  residues  that  have  four  peaks  that  are  the  most  well-separated  (unlike 
residue  30  or  35,  where  these  two  peaks  merge  into  one).  These  are  residues 
26,  28,  and  31.  Plots  of  the  number  of  steps  at  each  z-value  for  different 
mutations  for  one  of  these  residues,  residue  26,  are  shown  in  Fig.  7. 

Our  goal  was  to  find  the  number  of  steps  in  each  conformation.  Having 
chosen  the  residues  where  the  peaks  are  the  best  separated,  it  is  necessary  to 
find  the  number  of  steps  under  each  peak.  Since  they  overlap,  we  find  that 
the  best  way  to  do  this  is  to  use  the  least-squares  method  to  find  fits  to  each 
peak  that  best  describe  the  whole  function.  This  is  probably  best  understood 
using  a  concrete  example,  for  which  we  choose  residue  26. 

For  residue  26,  we  know  there  are  four  peaks,  and,  as  seen  in  Fig.  7,  the 
locations  of  these  are  fairly  clearly  visible.  To  fit  these  peaks,  we  begin  by 
averaging  all  of  our  results  for  the  wild-type  form  of  A/340.  We  assume  that 
the  peaks  are  Gaussian.  Since  the  peaks  are  reasonably  well  separated,  we 
assume  initially  that  the  center  of  each  Gaussian  is  at  the  maximum.  Then, 
we  perform  a  least-squares  fit  of  the  standard  deviation  and  amplitude  of 
each  Gaussian.  Particularly,  we  perform  our  fit  by  looping  through  the  peaks 
and  suggesting  changes  first  to  the  standard  deviation  and  then  to  the 
amplitude.  Each  time,  we  try  both  increasing  and  decreasing  the  standard 
deviation  by  a  specified  step  size,  and  check  whether  this  improves  the 
quality  of  the  fit.  We  then  move  on  to  the  next  peak  and  do  the  same  thing, 
then  repeat  the  process  for  the  amplitude.  Then  we  reduce  the  step  size  and 
repeat  the  whole  process.  We  do  this  until  the  fit  can  no  longer  be  improved 
by  further  iterations  of  the  process.  Having  done  all  that,  we  also  then  try 
altering  the  locations  of  the  centers  of  the  Gaussians  slightly  to  see  if  it 
improves  the  fit.  Having  then  selected  optimal  standard  deviations  and  center 
locations,  we  store  those,  and  assume  that  the  shape  and  location  of  each 
peak  will  remain  the  same  for  other  simulations,  and  that  simply  the 
amplitudes  will  vary.  A  sample  such  fit  is  shown  in  Fig.  8.  As  is  apparent 
from  the  figure,  the  fit  is  good.  Therefore,  we  concluded  that  assuming  the 
peaks  are  Gaussian  is  sufficient,  at  least  to  get  a  reasonable  estimate  of  the 
number  of  steps  in  each  conformation. 

Having  done  this,  we  are  able  to  then  perform  a  four  parameter  fit  for 
every  other  data  set  for  every  A/340  mutation  data  set,  simply  by  fitting  the 
amplitudes  while  keeping  the  standard  deviations  and  peak  locations  fixed. 
We  have  also  tried  allowing  the  standard  deviations  to  vary  when  doing 
these  subsequent  fits  but  this  does  not  substantially  improve  the  quality  of 


(a)  A-Beta  40 

Residue  26 


(b)  A-Beta  42 

Residue  26 


FIGURE  8  Binned  number  of  steps  spent  at  each  z-coordinate  for  residue 
26  of  ( a )  A/340  and  ( b )  A/342,  along  with  fits  used  for  data  analysis.  Shown 
here  are  the  data  for  wild-type,  compared  to  fits  of  Gaussians,  as  described  in 
Appendix:  Data  Analysis.  The  Gaussians  appear  to  provide  fairly  good  fits 
and  thus  we  use  these  in  this  work  to  calculate  the  number  of  steps  in  each 
conformation. 

the  fit,  and  it  seems  reasonable  that  the  shape  of  a  given  peak  for  a  given 
residue  should  be  constant. 

We  apply  a  similar  technique  to  fitting  residues  28  and  3 1 ,  and  residue  30, 
but  we  fit  residue  30  using  only  three  peaks,  and  we  do  the  same  for  A/342. 
Additionally,  using  this  technique,  we  are  able  to  calculate  standard 
deviations  for  a  given  peptide  or  mutant.  For  example,  for  the  D23N  mutant, 
where  residue  23  is  changed  from  aspartate  (D)  to  asparagine  (N),  we 
perform  a  fit  separately  for  each  of  our  five  MC  runs  and  calculate  the  area 
under  each  peak.  Then,  we  average  the  results  over  all  five  trials  and 
calculate  the  standard  deviation. 

For  the  four-peak  case,  it  is  somewhat  difficult  to  accurately  separate  the 
third  and  fourth  peaks,  thus  it  helps  that  we  are  able  to  separately  perform  fits 
to  four  peaks  on  residues  26,  28,  and  31.  We  are  thus  able  to  take  the 
apparent  change  of  steps  under  each  of  these  three  peaks  relative  to  the  wild- 
type  and  average  over  all  three  residues.  This  significantly  reduces  our  error. 

There  is  one  more  factor  which  complicates  issues.  Going  back  to  our 
earlier  example  of  the  location  of  residue  28  in  Fig.  2,  it  appears  in  very 
different  positions  along  the  z  axis  depending  on  which  conformation  the 
peptide  is  in.  However,  it  is  also  possible  to  have  conformations  where 
residue  28  may,  for  example,  be  quite  low,  like  it  is  in  the  transbilayer 
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conformation  (Fig.  2  a)  but  residues  38 — 40  are  not  in  the  lower  leaflet.  This 
effect  gets  worse  the  further  the  residue  being  examined  is  from  residue  40, 
and  is  apparent  in  Figs.  6  and  7  as  the  area  under  the  transbilayer  peak  on 
those  residues  is  obviously  much  larger  than  the  area  under  the  same  peak  on 
residue  40,  simply  because  more  conformations  look  similar  to  the 
transbilayer  conformation. 

This  means  that  the  calculated  number  of  steps  under  each  peak  for 
residues  26,  28,  30,  and  31,  while  accurate,  are  not  really  to  be  taken  as 
a  measure  of  how  many  simulation  steps  the  peptide  spends  in  each  of  those 
conformations.  Since  the  peak  shapes  remain  constant  over  the  different 
mutations,  however,  these  do  accurately  describe  the  change  relative  to  the 
wild-type.  That  is,  an  increase  in  the  prevalence  of  the  conformation  of  Fig.  2 
c  relative  to  wild-type  for  a  given  mutation  is  correct,  while  the  absolute 
number  of  steps  in  that  conformation  may  not  be  accurate,  except  for  the 
conformation  of  Fig.  2  a ,  which  we  can  extract  accurately  from  residue  40. 

As  we  discussed  above,  for  many  of  the  peptides  we  have  studied  results 
are  averaged  over  five  trials.  However,  for  some  we  have  used  significantly 
more  trials.  Specifically,  for  wild-type  of  A/340  and  A/342,  we  used  more 
trials,  as  it  was  particularly  important  to  have  a  good  average  for  those  results 
since  we  compare  all  of  our  other  results  to  those.  For  A/340  we  used  25 
trials;  for  A/342  we  used  20.  Additionally,  for  mutants  with  particularly 
small  changes  relative  to  wild-type,  or  particularly  large  standard  deviations, 
we  also  used  more  runs.  We  did  20  trials  for  A21G  A/340;  10  for  A21G 
A/342;  and  10  each  for  A/342  E22Q  and  D23N.  Of  the  artificial  mutants,  we 
did  10  trials  for  A/340  and  A/342  of  both  I32S  and  H6R. 
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Structure  of  infectious  prions:  stabilization 
by  domain  swapping 
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ABSTRACT  A  candidate  structure  for  the  minimal 
prion  infectious  unit  is  a  recendy  discovered  protein 
oligomer  modeled  as  a  (3-helical  prion  trimer  (BPT); 
BPTs  can  stack  to  form  cross-p  fibrils  and  may  provide 
insight  into  protein  aggregates  of  other  amyloid  dis¬ 
eases.  However,  the  BPT  lacks  a  clear  intermonomer 
binding  mechanism.  Here  we  propose  an  alternative 
domain-swapped  trimeric  prion  (DSTP)  model  and 
show  with  molecular  dynamics  (MD)  that  the  DSTP  has 
more  favorable  intermonomer  hydrogen  bonding  and 
proline  dihedral  strain  energy  than  the  BPT.  This  new 
structural  proposal  may  be  tested  by  lysine  and  N 
terminus  fluorescent  resonance  energy  transfer  (FRET) 
either  direcdy  on  recombinant  prion  protein  amyloid 
aggregates  or  on  synthetic  constructs  that  contain  the 
proline/lysine-rich  hinge  region  critical  for  domains  to 
swap.  In  addition,  the  domain  swapping  may  provide  1 ) 
intrinsic  entanglement,  which  can  contribute  to  the 
remarkable  temperature  stability  of  the  infectious 
prion  structure  and  help  explain  the  absence  of  PrPSc 
monomers,  2)  insight  into  why  specific  prolines  are 
potentially  relevant  to  three  inherited  forms  of  prion 
disease,  and  3)  a  simple  explanation  of  prion  strains 
assuming  the  strain  is  encoded  in  the  monomer  number 
of  the  oligomers. — Yang,  S.,  Levine,  H.,  Onuchic,  J.  N., 
Cox,  D.  L.  Structure  of  infectious  prions:  stabilization 
by  domain  swapping.  FASEB  J.  19,  1778-1782  (2005) 

Key  Words:  prion  •  domain  swapping  ■  (3  -helix  ■  amyloid 


MODEL  (5-HELICAL  PRION  TRIMER 

The  protein-only  model  for  prion  diseases  remains 
the  leading  candidate  to  explain  the  infectious  form  of 
the  disease  (1),  and  this  has  been  significantly  strength¬ 
ened  by  the  successful  production  of  disease  from 
synthetic  prion  protein  aggregates  (2).  Clearly,  the 
understanding  of  prion  disease  pathology,  propagation 
within  the  body,  and  therapeutic  approaches  will  be 
greatly  enhanced  by  a  detailed  knowledge  of  the  struc¬ 
ture  of  the  minimal  infectious  protein  unit. 

The  discovery  of  areal  aggregates  of  infectious  prion 
protein  oligomers  in  purified  brain  extracts  from  mice 
has  inspired  a  leading  candidate  model  for  the  minimal 
infectious  unit,  a  trimer  in  which  large  portions  of  the 
N  terminus  are  converted  to  left-handed  p-helices  (3, 


4).  We  call  this  the  |3-helical  prion  trimer  (BPT)  model. 
Figure  1A  shows  a  top  view  of  the  BPT  in  the  (3-helical 
region.  This  BPT  displays  a  3-fold  symmetry  consistent 
with  maps  of  the  electron  density  difference  between  the 
PrP  27-30  oligomers  and  the  “miniprion”  PrPSc106  oli¬ 
gomers.  The  BPT  can  easily  be  stacked  to  fonn  filaments, 
which  can  then  be  wound  into  amyloid  fibrils.  The  model 
trimer  also  shares  detailed  structural  features  with  the 
bacterial  protein  carbonic  anhydrase  trimer  (5).  Inspired 
in  part  by  this  model,  others  have  explored  the  possibility 
that  the  [3-helix  plays  a  role  in  other  amyloid  structures 
such  as  aggregates  of  polyglutamine  peptides  (6),  (3-amy¬ 
loid  peptides  (7),  and  yeast  prion-like  proteins  (8). 

However,  it  is  not  clear  in  the  BPT  model  what  holds 
the  trimer  together.  The  carbonic  anhydrase  trimer, 
which  was  the  template  for  the  BPT  model  is  in  fact 
held  together  by  intermonomeric  covalent  bonding  of 
zinc  ions  (5),  which  are  not  present  in  PrPSc.  It  was 
suggested  that  a  hydrogen  bonding  network  is  respon¬ 
sible  for  this  binding  (4) ,  as  in  the  case  of  the  bacterial 
(3-helical  trimer  serine  acyltransferase  (SAT)  (9),  but 
no  effort  was  made  to  test  this  idea  quantitatively. 
Moreover,  the  yellow  loops  of  the  BPT  model  (which 
are  highly  conserved  (10))  shown  in  Fig.  1  A,  B  and 
expanded  in  Fig.  1 C  contain  two  prolines,  and  prolines 
tend  to  disfavor  this  type  of  bending. 


HYPOTHESIS:  A  DOMAIN-SWAPPED  PRION 
TRIMER  IS  THE  MINIMAL  INFECTIOUS  UNIT 

We  hypothesize  that  a  stable  oligomeric  structure  for 
the  infectious  prion  is  the  domain-swapped  trimeric 
prion  (DSTP)  of  Fig.  IB,  which  can  increase  hydrogen 
bonding  and  reduce  elastic  energy  relative  to  the  BPT 
structure  proposed  earlier.  Further,  we  conjecture  that 
the  DSTP  structure  may  be  relevant  to  1)  the  unusual 
temperature  stability  of  the  prion  protein  and  the  lack 
of  observation  of  PrPSc  monomers  (both  due  to  the 
entanglement  of  the  DSTP  structure),  2)  the  kinetics  of 
at  least  three  inherited  forms  of  human  prion  disease, 
and  3)  the  encoding  of  prion  strains  in  conformation. 


1  Correspondence:  CTBP,  UC  San  Diego,  MC  0374,  9500 
Gilman  Dr.,  La  Jolla,  CA  92093-0374,  USA.  E-mail:  cox 
@  physics  .ucdavis.edu 
doi:  10.1096/fj.05-4067hyp 
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Figure  1.  A)  Proposed  fi-helical  structure  of  the  truncated 
PrPSc106  prion  trimer  (BPT)  taken  from  ref  4  (residues 
G89-F140).  B)  Structural  model  of  domain-swapped  trimeric 
prion  (DSTP)  of  the  present  paper.  C)  The  proposed  hinge 
region  expanded  from  the  yellow  loops  in  panels  A.  B. 

DOMAIN  SWAPPING:  MOTIVATION  AND 
RESULTS 

Concern  with  stability  of  the  BPT  model  led  us  to 
modify  it  by  domain  swapping  (11-13),  especially  since 
prolines  are  present  in  the  hinge  or  bridging  regions  of 
several  proteins  that  domain  swap  (14).  For  example,  in 
pl3sucl,  the  strain  on  two  prolines  in  the  hinge  mod¬ 
ulates  the  equilibrium  between  monomers  and  domain- 
swapped  dimers  (15,  16). 

We  constructed  our  DSTP  model  of  Fig.  IB  from  the 
proposed  BPT  structure  (we  thank  C.  Govaerts  for 
providing  his  coordinate  hie)  by  reorienting  the  mono¬ 
mers  to  place  the  loop  regions  near  the  trimer  center 
and  swapping  residues  89-99  to  a  neighboring  (3-helix 
in  a  cyclical  fashion,  with  the  BPT  loops  serving  as 
hinges.  This  DSTP  model  straightens  these  loops,  adds 
more  possible  hydrogen  bonds  in  the  loop  region,  and 
preserves  the  putative  monomer-monomer  spacing  and 
3-fold  symmetry  of  the  electron  density  maps  (3).  We 
ran  all-atom,  explicit  solvent  MD  (see  Appendix  for 
details)  on  energy  minimized  DSTP  and  BPT  structures 
using  AMBER8  (http://amber.scripps.edu).  As  a  con¬ 
trol,  we  applied  the  same  protocol  to  the  serine  acyl- 
transferase  (SAT)  (3-helical  trimer,  which  has  no 
strained  prolines  and  for  which  there  is  direct  crystal¬ 
lographic  evidence  for  intermonomer  hydrogen  bond- 
mediated  cohesion  (9).  Wee  also  carried  out  two  other 
simulation  tests:  1)  we  examined  the  stability  of  the 
DSTP  model  with  the  inclusion  of  the  C  terminus 
a-helices  proposed  in  ref  4,  and  2)  we  carried  out  three 
distinct  1-ns  runs  for  both  the  BPT  and  DSTP  structures 
to  make  sure  there  is  no  special  circumstance  with 
regard  to  individual  simulations. 


We  have  concluded  that  the  DSTP  is  significantly  more 
stable  than  the  BPT  and  therefore  a  more  promising 
candidate  structure.  First,  we  have  enumerated  direct  and 
water-bridged  hydrogen  bonds  between  monomers  in  die 
BPT  model,  our  DSTP  model,  and  the  SAT  structure  (see 
Appendix),  and  carried  this  out  for  three  separate  MD 
runs  on  each  structure.  We  show  the  direct  hydrogen 
bond  counts  for  three  different  1-ns  runs  of  the  BPT  and 
DSTP  structures  in  Fig.  2.  For  the  BPT,  we  find  2  ±  1, 
5  ±  1,  and  2  ±  1  direct  hydrogen  bonds  for  the  three  runs 
(averages  over  the  hist  100-ps  of  simulation  time),  and  at 
the  last  simulation  time  0,  0,  and  1  water-mediated  hydro¬ 
gen  bonds.  For  the  DSTP,  we  find  for  all  three  runs  9  ±  1 
direct  hydrogen  bonds  for  the  hist  100-ps  of  simulation 
time  and  2,  1,  and  1  water-mediated  hydrogen  bonds  in 
the  last  simulation  snapshot.  In  the  one  BPT  case  with  a 
relatively  high  number  of  direct  hydrogen  bonds,  the  N 
terminus  from  one  monomer  has  contacted  with  the  N 
terminus  region  of  another.  In  general,  the  BPT  hydro¬ 
gen  bonds  form  around  Q89  only,  while  the  DSTP 
hydrogen  bonds  form  in  the  hinge  region.  Note  that  Fig. 
2  shows  direct  hydrogen  bonds  formed  in  the  monomer- 
monomer  interface  for  the  DSTP  to  genetically  increase 
with  time  while  those  from  the  BPT  generically  decrease 
with  MD  time.  In  contrast  to  the  BPT,  the  simulated 
(3-helical  portion  (F140-A237)  of  the  SAT  trimer  has  11 
direct  and  3  water-bridged  hydrogen  bonds  (after  two 
stage  equilibration  with  no  restraints) . 

Second,  we  remark  that  the  domain  swapping  also 
relaxes  the  P101  and  P104  dihedral  energy  by  ~1.0 
kcal/mol  after  1-ns  of  MD.  Hence,  the  domain  swap¬ 
ping  helps  relieve  elastic  strain  as  anticipated. 

Third,  we  have  verified  that  the  DSTP  structure  for 
the  PrPs<  1  06  can  accommodate  the  C  terminus  a-h eli¬ 
tes  without  affecting  stability.  This  is  a  nontrivial  point, 
since  a  rigid  rotation  of  the  monomers  within  the  BPT 
is  not  possible  due  to  steric  constraints.  As  shown  in  Fig. 
3,  we  find  that  by  “peeling  off’  one  triangular  edge 
(residues  133-140)  of  the  lower  (C  terminus)  end  of 
each  (3-helix  to  form  a  random  coil  segment  that  we  can 
link  to  the  a-helices  that  begin  at  residue  177,  with  no 
resultant  change  in  direct  hydrogen  bonding  structure. 

Finally,  high  temperature  MD  provides  direct  evi¬ 
dence  for  BPT  instability;  at  500K,  the  DSTP  remains 
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Figure  2.  MD  simulations  show  that  the  nonswapped  BPT  (A) 
has  less  direct  hydrogen  bonds  in  the  interface  region  than 
the  domain-swapped  DSTP  (B) .  In  both  cases  we  have  simu¬ 
lated  residues  G89-F140.  Different  colors  correspond  to  three 
different  1  ns  MD  runs  in  each  case  initialized  at  different 
initial  temperatures  (see  Appendix) . 
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Figure  3.  A)  Model  for  DSTP  PrPSc106  (with  residues  G141- 
H176  deleted)  and  a-helices  included.  Residues  133-140  have 
been  removed  from  the  lower  rung  of  the  (3-helices  for  each 
monomer,  converted  to  random  coil,  and  attached  to  residue 
177  of  the  a-helices.  B)  Number  of  direct  hydrogen  bonds  for 
the  DSTP  structure  with  intact  a-helices.  The  a-helices  have  no 
effect  on  the  stabilization  induced  by  the  domain  swapping. 


compact  and  entangled,  while  the  BPT  monomer- 
monomer  separation  nearly  doubles.  Principal  compo¬ 
nent  analysis,  which  identifies  the  amplitude  and  char¬ 
acter  of  the  dominant  fluctuations,  shows  similar  BPT 
unbinding  tendencies  at  300K  as  in  the  500K  simula¬ 
tion,  but  no  unbinding  tendencies  for  the  DSTP.  Figure 
4  shows  the  radius  of  gyration  in  the  500K  runs  as  a 
function  of  simulation  time.  Within  the  1-ns  window  of 
simulation,  the  BPT  radius  appears  to  keep  growing, 
reaching  about  twice  that  of  the  apparently  more  stable 
DSTP  structure.  This  stabilization  through  entangle¬ 
ment  may  be  relevant  to  the  unusual  temperature 
stability  of  infectious  prion  proteins,  although  we  can¬ 
not  rule  out  that  such  stabilization  is  related  to  larger 
scale  aggregates  of  prion  oligomers.  The  entanglement 
may  also  be  relevant  to  the  inability  to  isolate  mono¬ 
meric  forms  of  PrPSc,  since  the  scrapies  form  is  stabi¬ 
lized  in  this  picture  only  in  aggregate. 


POSSIBLE  RELEVANCE  TO  INHERITED 
DISEASES  AND  PRION  STRAINS 

Finally,  we  discuss  the  possible  relevance  of  the  DSTP 
model  to  inherited  diseases  and  strains.  First,  proline 
point  mutations  yield  heritable  Creutzfeldt-Jakob  dis¬ 
ease  (P105T  for  humans,  P104T  for  mice)  or  Gerst- 
mann-Straussler-Scheinker  disease  (P102L,  P105L  for 


humans,  P101L,  P104L  for  mice;  see  chapter  14,  ref  1). 
Transgenic  mice  with  P101I  have  even  greater  disease 
susceptibility  whereas  those  with  PlOlYhave  somewhat 
less  susceptibility  than  PI  OIL  (K.  Nazor  and  G.  Telling, 
private  communication).  We  note  that  the  centrality  of 
these  sites  in  the  DSTP  hinges  affords  possible  burial  of 
the  relatively  hydrophobic  I,  L,  T,  and  Y  side  chains  and 
a  potential  speed  up  of  trimer  formation  through 
increased  hinge  flexibility.  Moreover,  the  insertion  of 
relatively  hydrophobic  residues  into  the  unstructured 
region  of  the  normal  PrPc  protein  will  raise  the  energy 
of  this  structure.  Indeed,  it  is  already  known  that  the 
P102L  mutation  is  marginally  less  stable  than  the 
wild-type  protein  (see  p.  694,  ref  1  ).  A  decrease  in  the 
stability  of  the  cellular  protein  form  relative  to  wild-type 
monomers  coupled  to  an  increase  in  stability  of  the 
scrapies  form  (associated  with  more  effective  hydro- 
phobic  burial  of  the  mutant  residue)  should  lead  to 
more  rapid  conversion  kinetics  via  ramping  of  the 
overall  free  energy  surface  slope  in  favor  of  the  scrapies 
form. 

In  Fig.  5  we  illustrate  two  ways  in  which  strains  might 
breed  true  if  the  monomer  content  within  an  oligomer 
encodes  the  strain  information.  With  this  assumption, 
strains  can  breed  true  either  by  1)  oligomer  conversion 
and  crystallization  on  membranes  (Fig.  5A),  with  local 
shape  matching  favoring  homogeneous  growth  in  a 
manner  analogous  to  the  oriented  aggregation  of  nano¬ 
particles  (17),  or  2)  via  interneuronal,  epitaxially  tem- 
plated  conversion  across  the  synaptic  cleft  or  interneu¬ 
ronal  gaps  (Fig.  5 B).  The  plausibility  of  the  latter 
conjecture  rests  on  the  length  of  unstructured  N  ter¬ 
minus  PrPc,  which  is  ~100  residues,  or  ~30  nm 
stretched  out,  easily  sufficient  to  span  the  interneuro¬ 
nal  gaps  of  10-20  nm.  Both  of  these  models  would 
enjoy  enhanced  conversion  kinetics  in  the  presence  of 
unconverted  or  “soft”  oligomers  of  PrPc  (18);  the 
interneuronal  mechanism  in  this  case  would  allow  for 
exponential  growth  via  oligomeric  autocatalytic  conver¬ 
sion.  We  note  that  a  different  domain  swapping  model 
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Figure  4.  The  entanglement  due  to  domain  swapping  intro¬ 
duces  additional  stability  for  the  DSTP  relative  to  the  BPT  as 
revealed  by  high  temperature  (T=500K)  simulations,  where 
the  radius  of  gyration  of  the  BPT  increases  more  rapidly  than 
that  of  the  DSTP,  which  reflects  unbinding  of  the  BPT. 
Principal  component  analysis  at  300K  reveals  a  similar  un¬ 
binding  tendency. 
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Figure  5.  Schematic  models  for  strains  breeding  true  under 
the  assumption  that  strain  is  encoded  in  monomer  number. 
For  simplicity,  we  assume  one  strain  is  encoded  as  trimers, 
one  as  tetramers.  A)  Areal  aggregation.  Given  local  size  or 
orientation  mismatches  between  trimers  (upper  image)  and 
tetramers  (lower  image),  energetics  will,  for  example,  disfa¬ 
vor  growth  of  trimers  off  of  tetramer  seeds  (lower  image).  B) 
Interneuronal  templating.  Assuming  that  seeds  template  con¬ 
version  across  interneuronal  gaps  or  the  synaptic  cleft,  trim¬ 
ers  will  seed  trimer  growth  (upper  image) ,  tetramers  will  seed 
tetramer  growth  (middle  image),  and  there  will  be  no  cross- 
seeding  between  trimers  and  tetramers. 

for  strains  has  been  conjectured  previously  (19);  sepa¬ 
rately,  domain-swapped  prion  dimers  have  been  crystal¬ 
lized  (20).  The  latter  result  is  consistent  with  prion 
oligomerization  via  intermolecular  disulfide  bonding  as 
has  been  shown  by  redox  induced  prion  hbrillization  in 
vitro  (21),  although  only  intramolecular  disulfide 
bonding  has  so  far  been  observed  in  vivo  (22).  Other 
specific  domain  swapping  mechanisms  than  those  we 
have  discussed  here  may  play  a  critical  role  in  amyloid 
fibril  formation  for  prions  as  has  been  observed  for 
other  proteins  (24-26).  For  the  PrP  27-30  scrapies 
form,  the  large  loop  (residues  E145-R163)  of  ref  4) 
made  from  some  of  the  deleted  residues  of  the 
PrPSc106  structure  could  serve  as  a  hinge  region  for 
alternative  domain-swapped  structures. 


PROPOSED  EXPERIMENTAL  TESTS 

While  the  DSTP  structure  may  ultimately  be  testable 
direcdy  on  recombinant  infectious  prion  oligomers,  as 


have  been  synthesized  by  Prusiner  and  collaborators 
recently  (2),  we  first  propose  testing  an  analog  system 
by  FRET  labeling.  Poly-Q  peptides  are  of  considerable 
interest  because  of  the  polyglutamine  diseases  (espe¬ 
cially  Huntington’s).  The  “critical”  length  near  36 
glutamine  repeats  for  disease  is  suggestive  of  a  left- 
handed  (3-helix  structure,  as  discussed  elsewhere  (4,  6, 
27).  The  key  experiment  we  propose  is  to  construct 
synthetic  peptides  of  the  form  Q12(KPSKPK)Qm, 
where  m  =  24-36.  The  goal  is  to  label  the  N  terminus 
Q’s  in  half  the  peptides  with  a  donor,  the  K’s  with  a 
different  donor,  and  in  the  other  half  put  in  N  termi¬ 
nus  Q  acceptors  and  K  acceptors.  If  the  argument  about 
stress  relaxation  and  hydrogen  bond  stabilization  works 
for  the  prion,  it  should  also  work  here.  There  should  be 
intense  K-K  FRET  given  intermonomer  separations  of 
the  side  chains  by  ~4  A  if  the  domain  swapping  arises. 
There  will  also  be  weak  N  terminus-Q-  N  terminus-Q 
FRET.  Presumably  there  could  be  an  N  terminus-Q-to- 
N  terminus-Q  FRET  signal  arising  with  time  if  the 
trimers  are  first  formed  and  then  they  fibrillize,  since 
the  corresponding  distances  would  get  reduced  by 
~60%  case  assuming  stacking  of  the  [3-helices.  If  this 
works  on  the  synthetic  peptide  it  should  be  attempted 
on  recombinant  prion  protein  prior  to  in  vitro  aggre¬ 
gation  as  per  ref  2.  We  would  anticipate  little  interfer¬ 
ing  contributions  from  other  lysines  in  the  prion  pro¬ 
tein  sequence  that  are  far  from  the  hinge  region. 


SUMMARY 

We  propose  a  domain-swapped  trimer  model  for  the 
minimal  infectious  prion  oligomer,  and  note  that  the 
domain  swapping  stabilizes  the  oligomer  relative  to  the 
earlier  proposed  model  by  reducing  stress  in  proline 
containing  loops,  increasing  intermonomer  hydrogen 
bonding,  and  promoting  entanglement.  We  conjecture 
that  domain  swapping  may  play  a  role  in  explaining  the 
temperature  stability  of  the  scrapies  form  of  the  prion 
protein,  the  absence  of  PrPSc  monomers,  inherited 
prion  disease  phenotypes  associated  with  mutating  pro¬ 
lines  in  hinges  of  the  DSTP,  and  that  conformation  may 
be  encoded  in  the  monomer  number  in  a  given  oli¬ 
gomer.  Finally,  we  propose  experimental  tests  based  on 
FRET  labeling  of  principally  lysine  residues  near  the 
crucial  pralines,  which  should  produce  intense  FRET 
response  if  the  domain  swapping  is  present.  pj| 


Appendix 

The  simulations  were  carried  out  with  the  AMBER8  molecu¬ 
lar  dynamics  package  (28)  with  the  AMBER  parm99  force 
field  and  explicit  TIP3P  waters,  following  these  procedures.  1 ) 
Before  we  start  a  molecular  dynamics  simulation,  we  perform 
an  energy  minimization  of  2000  steps  to  partially  relax  the 
entire  molecular  system.  2)  We  perform  a  two-stage  equilibra¬ 
tion  to  further  relax  the  protein  and  the  surrounding  solvent. 
In  the  first  stage,  we  start  the  system  from  a  low  temperature 
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of  approximately  100K  and  gradually  heat  up  to  300K  over 
20-ps  of  simulation  time.  We  perform  this  stage  of  equilibra¬ 
tion  with  the  volume  held  constant.  In  the  second  stage,  we 
equilibrate  the  system  using  pressure  and  temperature  con¬ 
trol  to  adjust  the  density  of  water  to  experimental  values.  We 
have  equilibrated  the  system  for  a  total  of  40-ps.  During  the 
first  two  steps,  the  backbone  heavy  atoms  of  the  (3-helical  (and 
also  the  a-helical  if  applicable)  portions  are  positionally 
restrained  using  a  harmonic  potential.  The  hinge  loops  and 
the  regions  linking  the  (3-helical  and  a-helical  domains  are 
allowed  to  freely  move.  To  achieve  the  three  different  MD 
runs  for  the  BPT  and  DSTP  structures,  we  start  with  different 
initial  temperatures  that  seed  a  different  run  through  a 
different  initial  thermal  distribution  of  atom  velocities.  Fi¬ 
nally,  we  carry  out  a  1-ns  production  run  at  constant  pressure 
and  temperature  with  no  positional  restraints.  For  the  high 
temperature  (T=500K)  simulations,  the  production  runs  are 
performed  at  constant  volume  with  a  reduced  time  step  of 
integration  to  prevent  any  possible  system  blowup. 

Hydrogen  bonds  are  calculated  based  on  the  following 
criteria:  1)  direct  hydrogen  bonds  in  interface  are  counted 
with  a  donor-acceptor  distance  of  3.5  A  and  a  hydrogen- 
donor-acceptor  angle  of  60°,  and  2)  water-bridged  hydrogen 
bonding  in  interface  is  counted  if  a  water  molecule  couples  by 
criterion  1 )  to  an  atom  from  each  of  two  monomers. 
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Prion  Disease:  Exponential  Growth  Requires  Membrane  Binding 
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ABSTRACT  A  hallmark  feature  of  prions,  whether  in  mammals  or  yeast  and  fungi,  is  exponential  growth  associated  with 
fission  or  autocatalysis  of  protein  aggregates.  We  have  employed  a  rigorous  kinetic  analysis  to  recent  data  from  transgenic 
mice  lacking  a  glycosylphosphatidylinositol  membrane  anchor  to  the  normal  cellular  PrPc  protein,  which  show  that  toxicity 
requires  the  membrane  binding.  We  find  as  well  that  the  membrane  is  necessary  for  exponential  growth  of  prion  aggregates; 
without  it,  the  kinetics  is  simply  the  quadratic-in-time  growth  characteristic  of  linear  elongation  as  observed  frequently  in  in  vitro 
amyloid  growth  experiments  with  other  proteins.  This  requires  both:  i),  a  substantial  intercellular  concentration  of  anchorless 
PrPc,  and  ii),  a  concentration  of  small  scrapies  seeding  aggregates  from  the  inoculum,  which  remains  relatively  constant  with 
time  and  exceeds  the  concentration  of  large  polymeric  aggregates.  We  also  can  explain  via  this  analysis  why  mice 
heterozygous  for  the  anchor-full/anchor-free  PrPc  proteins  have  more  rapid  incubation  than  mice  heterozygous  for  anchor-full/ 
null  PrPc,  and  contrast  the  mammalian  membrane  associated  fission  or  autocatalysis  with  the  membrane  free  fission  of  yeast 
and  fungal  prions. 
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Prions  are  distinguished  from  other  amyloid  diseases  both 
by  their  infectious  character  and  the  observed  exponential 
growth  of  infectious  material  in  vivo  (1).  There  is  a  corre¬ 
spondence  to  this  of  prion-like  proteins  in  yeast  and  fungi, 
for  which  spontaneous  fission  is  reported  in  vitro  (2).  In  fact, 
it  has  been  argued  that  the  replication  necessary  for  infection 
in  mammals  and  non-Mendelian  inheritance  in  yeast/fungi 
requires  the  fission  or  autocatalysis  that  drives  the  exponen¬ 
tial  growth  (1).  It  has  been  established  for  yeast  prions  that 
additional  chaperone  proteins  most  likely  facilitate  the  fis¬ 
sion  of  aggregates  in  living  cells  ( 1 ).  It  is  an  open  question 
what  mechanism  drives  the  exponential  growth  in  mammals. 
Here  we  show  by  a  rigorous  kinetic  analysis  of  recent  disease 
time  course  data  that  the  exponential  growth  is  tied  to  mem¬ 
brane  anchoring  of  the  prion  protein,  suggesting  that  either 
mechanical  fission  of  areal  prion  aggregates  or  oligomeric 
autocatalysis  of  membrane  bound  prions  explain  the  ob¬ 
served  behavior. 

Chesebro  et  al.  (3)  recently  studied  transgenic  (Tg)  mice 
lacking  a  GPI  membrane  anchor  to  the  normal  cellular  PrPc 
protein  and  discovered  that  these  mice  grew  infectious  prions 
without  suffering  neuronal  death.  We  denote  these  anchor¬ 
less  cellular  prions  as  PrPjg  and  anchor-full  wild-type  (WT) 
cellular  prions  by  PrP^,T.  When  inoculated  with  infectious 
scrapies  prions  (PrPSc)  at  a  dose  that  induces  clinical  symp¬ 
toms  within  140-160  days  for  WT  mice,  the  Tg  mice  were 
symptomless  up  to  400-600  days,  even  though  proteinase 
resistant  PrP-res,  an  indicator  of  infectivity,  accumulated  and 
surpassed  the  maximal  WT  levels. 

In  Fig.  1  the  Tg  mice  PrP-res  concentration  ( crosses )  of 
Chesebro  et  al.  (3)  are  plotted  versus  the  square  of  time,  together 
with  a  linear  regression  fit  (line)  with  a  high  regression 


coefficient  ( R  =  0.97).  This  time  dependence  is  consistent 
with  short  time  kinetics  described  by  linear  polymer  elon¬ 
gation  via  monomer  addition  without  fission  or  autocatalysis 
(4),  illustrated  schematically  in  Fig.  2.  Assuming  the  PrP-res 
concentration  to  be  a  proxy  for  the  total  protein  content  in 
aggregate,  simple  kinetic  arguments  predict  a  behavior,  before 
monomer  depletion  and  seed  nucleus  depletion,  of 

[PrP-res]  (t)  «  (1/2)0+  [PrP?g])2[PrPf  ]  t\  (1) 

with  p+  [PrP!f  ]  the  elongation  rate  at  intercellular  monomer 
Tg  prion  concentration  |  PrPc[  g  ] ,  and  [PrP[jc]  is  the  intercel¬ 
lular  concentration  of  seeding  nuclei  from  inoculated  scra¬ 
pies  protein  after  initial  hydrodynamic  clearance.  The 
validity  of  Eq.  1  at  long  times  suggests  that:  i),  there  is  a 
substantial  homeostatic  concentration  of  intercellular  PrP(|g 
presumably  due  to  slow  clearance,  and  ii),  [PrP[]c]  is  hardly 
changed  implying  either  that  only  a  small  fraction  of  seeds 
grow  into  large  polymers  or  [PrP^c]  is  maintained  by  steady 
proteolytic  degradation  of  large  remnant  aggregates  from  the 
dose.  Given  a  similar  de  novo  production  rate  of  PrP^,T  in 
inoculated  WT  mice,  we  speculate  that  the  associated  sat¬ 
uration  of  [PrP-res]  arises  from  loss  of  PrP(NV  I  after  cell  death. 

This  elongation  hypothesis  is  testable  by:  i),  genetically 
engineering  mice  to  overexpress  PrP(|  g ,  which  will  quadrat- 
ically  modulate  the  PrP-res  concentration  (4),  and  ii),  by 
varying  the  initial  dose  of  PrPSc,  which  will  linearly  mod¬ 
ulate  the  PrP-res  concentration. 
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FIGURE  1  Quadratic-in-time  fit  to  infectious  prion  time  course 
data  of  Chesebro  et  al.  (3). 

Another  striking  observation  of  Chesebro  et  al.  (3)  was 
that  mice  heterozygous  for  expression  of  PrP^g  and  PrP^v| 
display  shorter  incubation  time  upon  inoculation  than  mice 
with  one  PrP(NVT  copy  and  one  inactive  fusion  construct.  We 
argue  that  this  is  due  to  an  enhanced  concentration  for  PrP(l  g 
relative  to  PrP(NV  |  and  that  PrP-res  obtained  from  PrP^g  also 
templates  PrP^j.  conversion.  The  latter  is  speculative  but 
should  be  testable. 

To  establish  the  expectation  that  the  PrP^,T  concentration 
is  lower  than  PrP^g  in  the  WT/Tg  heterozygotes,  it  is  suf¬ 
ficient  to  establish  that  PrP(NV  I  concentration  in  WT  mice  is 
lower  than  PrP.],  concentration  in  the  anchorless  Tg  mice.  To 
make  this  clear,  we  compare  estimates  for  cellular  prion 
concentrations  in  homozygous  WT  mice  with  homozygous 
Tg  mice.  Before  the  postulated  cell-death  driven  saturation 
of  infectious  material,  WT  mice  inoculated  with  a  concen¬ 
tration  [PrP^c]  of  scrapies  seeds  will  have  a  time  course  of  the 
general  form  (4) 

[PrP-res]  (f)  °c  A  +  (cosh[Ar]  —  1 ) ,  (2) 

where  A  =  ln(2)/f2  is  the  percentage  growth  rate  and  t2  is  the 
doubling  time.  The  coefficient  A+  can  be  determined  from 
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FIGURE  2  Schematic  model  for  linear  elongation  driven  growth 
of  infectious  prion  material  from  inoculated  seeds  and  anchor¬ 
less  cellular  prion  proteins. 


the  short  time  behavior  of  Eq.  2,  which,  like  Eq.  1,  is 
described  by  linear  elongation  given  by 

1/2A  +  (At)2  =  l/2(p+  [PrPwT])2[PrPnCk2,  (3) 

where  [PrP^v|  ]  is  the  homeostatic  concentration  of  mem¬ 
brane  bound  WT  PrPc  and  we  have  assumed  that  the  WT  and 
Tg  mice  have  the  same  PrP-res  elongation  coefficient  p+. 

On  the  other  hand,  at  long  times  (but  before  saturation) 
Eq.  2  gives 


[PrP-res]  (t)  x  1/2A+  exp(Ar).  (4) 

Now,  we  set  t+  =  A'd'  f2  as  the  time  it  takes  WT  mice  to 
reach  clinically  detectable  levels  of  PrP-res  concentration  at 
the  inoculum  level  generating  a  seed  nuclei  concentration 
[PrpSc],  wjlere  /\r  is  the  number  of  doublings  experienced  in 
that  process,  and  t_  is  the  time  it  takes  Tg  mice  to  reach  the 
same  clinical  concentration  of  PrP-res  for  the  same  initial  in¬ 
oculation  dose.  By  taking  suitable  ratios  to  eliminate  A+, 
[PrP®c],  and  p+,  the  ratio  of  homeostatic  concentrations  of 
cellular  prions  from  the  WT  mice  to  the  Tg  mice  is  given  by 

[PrPw-r]  /  [PrPig]  =  (l/2)ln(2)  Nd+  (f_/f  + ) 
Xexp(-ln(2)(/Vd+  /2)). 

From  Chesebro  et  al.  (3),  t+  =  150  days,  and  t_  =  400 
days.  A  reasonable  estimate  (5)  for  the  number  of  doublings 
is  N A  =  20  for  the  dose  of  Chien  et  al.  (1).  With  these 
numbers,  Eq.  5  gives  a  concentration  ratio  (and  hence 
elongation  rate  ratio)  for  WT/Tg  mice  of  0.036.  This  is 
reasonable  given  that  likely  slower  PrPc  clearance  in  the  Tg 
case  will  lead  to  a  higher  extracellular  concentration  of 
cellular  prion  protein.  By  employing  the  arguments  of  Chien 
et  al.  (1)  we  obtain  elongation  rate  values  of  0.13/day(Wt) 
and  3.5/day  (Tg).  The  former  is  in  good  agreement  with  es¬ 
timates  made  elsewhere  for  linear  elongation  based  upon 
analysis  of  dose-incubation  curves  (6). 

Hence,  in  the  WT/Tg  and  WT/null  heterozygotes  explored 
in  Chesebro  et  al.  (3),  we  anticipate  in  each  case  the  mem¬ 
brane  bound  PrP^j.  concentration  to  be  about  half  that  of  the 
homozygous  WT  mice,  whereas  the  intercellular  PrPjg 
concentration  should  be  about  half  that  of  the  homozygous 
Tg  mice.  Because,  as  shown  with  infectious  prions  bound 
to  electrodes  (7),  templating  and  conversion  can  be  driven 
by  scrapies  material  not  bound  to  the  membrane  surface,  we 
expect  the  incubation  time  of  the  WT/Tg  heterozygotes  to  be 
significantly  accelerated  relative  to  the  WT/null  heterozy¬ 
gotes  as  is  observed. 

We  note  that  membrane  associated  exponential  growth 
might  be  due  to:  i),  as  yet  undiscovered  membrane  specific 
enzymes  splitting  aggregates,  in  analogy  to  the  role  of 
Hp  104a  for  yeast  prions  ( 1 );  ii),  mechanical  breakage  of  aggre¬ 
gates  due  to  membrane  curvature  or  membrane  undulations 
(8);  iii),  oligomeric  autocatalysis  arising  from  interneuronal 
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templating  by  infectious  oligomeric  seeds  bound  to  one  or  the 
other  membrane  (9). 
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ABSTRACT  We  employ  ab  initio  electronic  structure  calculations  to  obtain  two  structural  models  for  copper  bound  in  the 
strongest  binding  site  of  the  noninfectious  form  of  the  prion  protein.  The  models  are  compatible  with  available  experimental 
constraints  from  electron  spin  resonance  data.  The  bending  of  the  peptide  backbone  attendant  with  the  copper  binding  is  not 
compatible  with  the  requisite  straight  beta-strand  backbone  structure  for  the  same  sequence  contained  in  two  recently 
proposed  models  of  the  prion  protein  structure  in  its  infectious  form.  We  hypothesize  that  copper  binding  at  this  site  is 
protective  against  conversion  to  the  infectious  form,  discuss  experimental  data  which  appear  to  support  and  conflict  with  our 
hypothesis,  and  propose  tests  using  recombinant  prion  protein,  genetically  modified  cultured  neurons,  and  transgenic  mice. 
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The  functional  role  of  the  normal  cellular  prion  protein 
(PrPc)  may  well  be  related  to  copper  binding,  potentially 
protecting  against  oxidative  damage  in  the  synaptic  region 
by  sequestering  divalent  copper1.  This  suggestion  is 
supported  by  data  for  transgenic  knockout  mice  devoid  of 
the  gene  for  expressing  PrPc;  these  mice  appear  to  suffer 
late  stage  oxidative  degradation  in  the  neuronal  synaptic 
regions  where  surface  bound  prion  protein  is  preferentially 
concentrated  in  wild  type  mice2'1.  The  strongest  copper 
binding  site  is  in  the  protein  region  which  converts  to  beta 
sheet  structure  in  the  infectious  (PrPSt)  form1.  We  use 
electronic  structure  calculations  to  study  two  possible 
geometries  for  the  bound  copper-PrPc  complex;  we  find 
that  these  geometries  are  incompatible  with  recently 
proposed  models45  for  PrPSc  oligomers,  and  we  thus 
hypothesize  a  mechanism  for  inhibition  of  PrPc-to-PrPSc 
conversion  via  copper  binding. 

Copper  binds  principally  to  mammalian  PrPc  at 
octarepeat  sites  of  highly  conserved  form  (though  not 
number)  in  a  now  well  understood  pyramidal  geometry 
with  binding  to  the  peptide  backbone,  histidine  side  chains, 
and  axial  waters6.  These  octarepeats  (residues  60-91  in 
humans)  are  not  essential  to  PrPSc  ,  which  remains 
infectious  even  after  proteinase  exposure  that  leaves 
residues  92-230  intact.  One  strong  copper  binding  site  is 
present  in  the  PrPSc  core  region,  containing  the  sequence 
92-96  GGGTH  for  humans.  Electron  spin  resonance  data 
suggest  there  is  binding  to  the  H96  (H95  in  mice)  side 
chain  and  the  G94  amide  group1.  It  is  not  known  whether 
the  primary  peptide  coordination  is  to  four  nitrogens 
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(NNNN)  or  to  three  nitrogens  and  an  oxygen  (NNON), 
although  the  former  structure  has  been  conjectured  to  be 
more  likely1. 

We  have  studied  this  problem  computationally  with  the 
SIESTA  local  orbital  based  density  functional  theory  code 
using  conjugate  gradient  based  energy  minimization  to 
examine  possible  GGGTH-Cu(II)  geometries7'8.  We  have 
considered  the  NNNN  structure  and  an  NNON  structure 
compatible  with  the  ESR  data1  and  our  results  are  shown  in 
the  upper  two  panels  of  the  figure  for  the  mouse  prion. 
We  built  initial  candidate  structures  of  the  form  Ac- 
Cu(H20)6(GGGTH)-NH2  using  ChemSite  Pro  (ChemSW, 
Fairfield  CA)  and  VMD4.  For  each  geometry,  the  lowest 
energy  sampled  provides  our  model  bound  copper-prion 
complex.  We  used  double-zeta  (DZ)  basis  sets  for  light 
atoms  (H,  C,  N,  O)  and  double-zeta  polarized  orbitals  for 
copper.  We  used  the  Troullier-Martins  norm-conserving 
pseudopotential.  We  employed  the  Perdew-Becke- 
Ernzerhof  based  generalized  gradient  approximation 
exchange  correlation  energy  functional.  Our  energy  cutoff 
for  matrix  element  integration  was  120  Rydberg.  All 
calculations  were  carried  out  in  a  periodic  unit  cell  size: 
30.7  A  x  30.7  A  x  30.7  A.  Geometry  minimization  was 
carried  out  to  a  force  tolerance  of  32pN. 

Our  main  result  is  inferred  from  comparing  the  two 
possible  GGGTH-Cu(II)  geometries  to  recently  proposed 
(3-helical  PrPSc  trimer  models4  5  We  see  that  the  backbone 
bending  induced  by  copper  binding  is  not  compatible  with 
formation  of  (3-strands  in  the  left  handed  helices.  (It  is 
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unclear  if  another  recently  proposed  model  has  the  same 
characteristics10.  )  Hence,  copper  in  the  non-octarepeat 
sequence  GGGTH  can  protect  against  PrPc-to-PrPSc 
conversion;  removal  of  this  copper  by  pH  reduction 
associated  with  synaptic  fluctuations  or  endocytosis  may  be 
a  key  step  in  the  conversion  pathway.  This  picture  is 
supported  by  the  observation  that  only  the  GGTH  sequence 
is  required  for  binding,  and  eutherian  prion  sequences 
reveal  high  conservation  of  either  GGTH  or  the  very 
similar  GGSH  sequence  in  the  converting  region11. 

Removal  of  this  copper  by  pH  reduction  may  be  a  key 
step  in  the  PrPc-to-PrPSc  conversion  pathway  since  it  is 
known  that  acidic  conditions  favor  detachment  of  the 
copper12.  We  note  two  routes  for  such  a  pH  reduction  with 
accompanying  conversion  of  PrPc.  First,  upon  endocytosis 
into  the  cell,  the  pH  is  reduced.  Assuming  that  some  PrPSc 
is  endocytosed  as  well  as  PrPc,  then  conversion  can  take 
place  within  the  endocytosed  compartments  after  the 
copper  is  expelled.  Second,  given  that  PrPc  is  preferentially 
concentrated  near  the  synapse,  then  from  signal  bearing 
vesicles  released  near  or  at  the  synapse13  can  reduce  the  pH 
locally  and  facilitate  copper  expulsion. 


Figure:  Potential  copper  binding  motifs  in  the  converting 
region  of  the  normal  (PrPc)  mouse  prion  protein  which  are 
consistent  with  ESR  data1  are  shown  in  panels  (A,B).  The 
corresponding  copper  free  stretch  of  the  left  handed  beta 
helix  model  of  the  infectious  (PrPSc)  protein  from  Ref.  (4) 
is  shown  in  (C). 

Three  experimental  results  support  our  proposed 
mechanism  indirectly:  (1)  Post-mortem  studies  show 
copper  depletion  in  regions  of  infection,  with  80-90% 


reduction  of  copper  content  in  the  corresponding  prion 
proteins14.  This  result  is  consistent  with  our  mechanism, 
since  this  GGGTH  binds  Cu(II)  first  with  potential 
cooperative  enhancement  of  subsequent  binding  to  other 
sites.  If  infection  promotes  removal  of  copper  from  this 
site,  then  it  is  more  likely  that  copper  will  not  be  retained  at 
the  other  sites  assuming  near-equilibrium  conditions.  (2) 
Copper  uptake  in  infected  cultured  neurons  is  suppressed 
approximately  10  fold  compared  to  control  cells13.  Again, 
this  is  consistent  with  the  hypothesis,  although  this  requires 
an  understanding  of  how  the  prion,  which  is  easily  shown 
to  contain  less  than  1%  of  cellular  copper,  can  shutdown 
copper  transport.  (3)  Copper  in  solution  inhibits  in  vitro 
growth  of  amyloid  fibrils1617.  This  is  a  more  direct 
corollary  to  our  hypothesized  protective  role  of  copper. 

Other  experiments  provide  mixed  support  for  our 
mechanism.  The  copper  chelator  D-penicillamine 
inoculated  in  vivo  delays  disease  onset,  which  is  at  apparent 
odds  with  our  hypothesis18.  The  chelator  cuprizone, 
however,  is  known  to  induce  a  spongiosis  (vacuolation  of 
brain  tissue)  very  similar  to  that  of  prion  disease  although  it 
also  induces  demyelination  of  neurons  and  no 
transmissibility19'20.  Moreover,  it  is  also  known  that  copper 
confers  proteinase  resistance  to  the  cellular  prion  protein  in 
the  weaker  binding  octarepeat  regions21.  This  is  significant 
for  two  reasons.  First,  proteinase  resistance  is  one  of  the 
generic  hallmarks  of  the  infectious  form  of  the  disease 
(although  not  all  infectious  prion  protein  is  in  fact 
proteinase  resistant).  Second,  the  octarepeat  deletion  does 
not  remove  infectivity  from  PrPSc.  However,  in  transgenic 
mice  without  octarepeats,  disease  incubation  upon  direct 
intercerebral  inoculation  is  slowed"".  Hence,  these 
octarepeats  are  not  necessary  for  disease,  but  (i)  they  can 
impart  structure  in  the  presence  of  copper  that  confers 
proteinase  resistance  to  the  prion  protein,  and  (ii)  their 
presence  enhances  disease  incubation.  Indeed,  it  is  known 
that  the  octarepeat  region  adopts  a  unique  structural  motif 
when  copper  is  bound23.  We  can  rationalize  these  data  by 
assuming  the  following  affinity  ranking:  cuprizone  >  non- 
octarepeat  site  >  D-PEN  >  octarepeat  site.  This  ranking  is 
plausible  given  the  hexadentate  chelation  of  cuprizone 
versus  at  most  tetradentate  coordination  of  D-PEN.  With 
this  assumption,  D-PEN  lifts  the  proteinase  resistant 
structure  associated  with  copper  binding  to  octarepeats,  and 
this  slows  conversion  to  the  infectious  form.  These 
assumptions  may  be  tested  by  performing  competitive 
binding  experiments  in  vitro  with  D-PEN18,  cuprizone19 
and  the  relevant  copper  binding  segments  of  the  prion 
protein.  We  note  that  the  high  affinity  of  cuprizone  for 
divalent  copper  apparently  leads  to  stripping  of  copper 
from  other  proteins  so  that  the  in  vivo  phenomena  need  not 
be  limited  to  prions;  if  our  assumptions  are  correct,  these 
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non-prion  related  effects  must  lead  to  more  rapid 
degeneration  of  neurons  than  conversion  to  PrPSt. 

To  further  test  the  mechanism,  it  is  desirable  to  perform 
mutation  studies  on  the  stretch  GGSH  or  GGTH.  The  key 
idea  is  to  mutate  out  the  H  for,  e.g.,  A,Y,G  or  other  amino 
acids  which  should  severely  mitigate  copper  binding.  We 
propose  first  that  recombinant  mutant  PrPc  from  residues 
90-230  (with,  e.g.,  H96A  in  the  human  form)  be  allowed  to 
aggregate  in  copper-ful  and  copper-free  environments. 
Using  the  reduced  length  prion  will  avoid  the  confounding 
effects  from  the  octarepeat  regions  discussed  above.  If  the 
mechanism  is  correct,  the  fibril  formation  rate  should  be 
approximately  the  same  as  wild  type  fibrils  grown  in 
copper  free  solutions.  Second,  one  can  transfect  cultured 
neurons  with  the  mutated  DNA  construct  and  look  for 
increased  susceptibility  to  infection.  Third,  assuming 
positive  results  from  the  fibril  growth  and  cultured  neuron 
experiments,  transgenic  mice  with  the  H95  mutated  away; 
these  mice  should  prove  to  have  shorter  incubation  times 
for  a  given  dose  than  wild  type  mice.  Finally,  fibrillization 
experiments  on  recombinant  marsupial  PrPc  with  and 
without  copper  should  be  informative  since  the  H  is 
replaced  with  Y24;  our  prediction  based  upon  the  present 
hypothesis  would  be  that  the  copper  should  have  reduced 
inhibition  of  fibril  formation  relative  to  growth  of 
recombinant  human  prion  protein. 
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Mammalian  prion  proteins  (PrP)  are  of  significant  public  health  interest.  Yeasts  have  proteins,  which  can 
undergo  similar  reconformation  and  aggregation  processes  to  PrP,  without  posing  a  threat  to  the  organism. 
These  yeast  “prions,”  such  as  SUP35,  are  simpler  to  experimentally  study  and  model.  Recent  in  vitro  studies 
of  the  SUP35  protein  found  long  aggregates,  pure  exponential  growth  of  the  misfolded  form,  and  a  lag  time 
which  depended  weakly  on  the  monomer  concentration.  To  explain  this  data,  we  have  extended  a  previous 
model  of  aggregation  kinetics  along  with  a  stochastic  approach.  We  assume  reconformation  only  upon  aggre¬ 
gation  and  include  aggregate  fissioning  and  an  initial  nucleation  barrier.  We  find  that  for  sufficiently  small 
nucleation  rates  or  seeding  by  a  small  number  of  preformed  nuclei,  the  models  achieve  the  requisite  exponen¬ 
tial  growth,  long  aggregates,  and  a  lag  time  which  depends  weakly  on  monomer  concentration.  The  spread  in 
aggregate  sizes  is  well  described  by  the  Weibull  distribution.  All  these  properties  point  to  the  preeminent  role 
of  fissioning  in  the  growth  of  misfolded  proteins. 

DOI:  10.1 103/PhysRevE. 72. 051915  PACS  number(s):  82.39.Fk 


I.  INTRODUCTION 

Bovine  spongiform  encephalopathy  (BSE)  in  cows, 
Scrapies  in  sheep,  Creutzfeldt-Jakob  Disease  (CJD)  and 
Kuru  in  humans  are  all  diseases  caused  by  a  specific  mis¬ 
folded  protein  residing  on  neurons  [1,2].  This  prion  protein 
(denoted  PrP'  in  its  normal  form  and  PiP"  in  its  misfolded 
form)  is  present  in  all  mammals  and  its  full  function,  to  date, 
is  still  unknown.  A  growing  body  of  evidence  strongly  im¬ 
plies  that  this  disease  propagates  not  by  nucleic  acids,  such 
as  DNA  or  RNA,  but  by  misfolded  proteins  (PrPvr)  [3,4]. 
In  our  current  understanding  of  prion  diseases,  the  autocata- 
lytic  misfolding  of  the  prion  proteins  plays  a  central  role. 
The  misfolded  form  PrP"  entices  normal  versions  of  the  cel¬ 
lular  prion  protein  (PrP'  )  to  change  conformation  to  the  mis¬ 
folded  or  disease  causing  form.  Although  PrPe  and  PrPIC 
have  the  same  amino  acid  sequence  (230  amino  acids  or 
residues),  PrPsc  has  a  higher  beta  sheet  content  than  its  nor¬ 
mal  form  [1].  In  humans,  the  disease  is  mostly  sporadic, 
perhaps  caused  by  a  rare  spontaneous  misfold  of  the  protein. 
But,  infectious  forms  are  also  known  to  occur,  such  as  in  the 
case  of  new-variant  CJD  caused  by  eating  BSE  infected 
meat.  In  the  latter  case,  infectious  agents  are  presumably 
misfolded  “seeds”  that  have  entered  the  body  from  outside. 
Developing  an  understanding  of  the  misfolding  process 
should  provide  insights  to  prevention  and/or  cure  of  these 
diseases. 

A  class  of  proteins  in  yeast  (SUP35)  undergo  a  conforma¬ 
tion  change  similar  to  the  mammalian  PrP.  However,  in  case 
of  yeast,  this  does  not  kill  the  organism.  In  fact,  reconfigured, 
aggregated  forms  lead  to  a  new  phenotype  [5].  Because  they 
are  not  toxic,  the  misfolding  is  much  easier  to  study  for  these 
proteins.  Besides,  it  has  been  shown  that  aggregates  pro¬ 
duced  in  vitro  can  lead  to  the  same  phenotype,  when  added 
to  a  yeast  cell  [6-8].  Thus,  it  is  possible  to  study  the  molecu¬ 
lar  aggregation  processes  for  yeast  prions  in  a  test  tube  in 
tremendous  detail,  and  a  large  quantity  of  experimental  data 
is  indeed  available  [9,10]. 


In  this  paper,  we  explore,  using  conventional  kinetic  as 
well  as  stochastic  models,  the  aggregation  dynamics  for 
model  yeast  prions  with  a  particular  goal  of  explaining 
in  vitro  data  [9,10]  including  (i)  linear  aggregates  which 
grow  by  monomer  addition,  (ii)  a  modest  concentration  de¬ 
pendence  of  initial  aggregate  growth  time  varying  roughly 
inversely  as  the  square  root  of  the  monomer  concentration, 
(iii)  an  observed  distribution  of  aggregate  sizes  with  a  mode 
that  is  large  (lOs-lOOs  of  monomers),  and  (iv)  a  sigmoidal 
growth  that  is  nearly  exponential.  We  extend  the  nucleation 
and  fissioning  model  of  Collins  et  al.  [9,10]  from  a 
moments-only  model  to  monitor  individual  polymer  lengths 
and  show  that  we  can  readily  obtain  the  above  listed  features 
in  the  data.  In  particular,  we  argue  against  a  model  based 
upon  rapid  equilibration  with  micelles,  which  can  also  pro¬ 
vide  a  weak  concentration  dependence  to  initial  aggregation 
times  [11,12],  but  not  the  exponential  growth  implied  by 
including  explicit  fission  of  aggregates.  The  fission  rate  plays 
a  central  role  in  determining  the  time  scales  as  well  as  the 
mean  aggregate  length.  The  range  of  aggregate  sizes  assumes 
a  Weibull  distribution,  also  well  known  from  other  processes 
where  fissioning  is  important.  Thus  the  agreement  with  ex¬ 
periments  highlights  the  importance  of  fissioning  in  the 
growth  dynamics. 

We  study  our  kinetics  both  with  continuous  time  model¬ 
ing  and  via  the  discrete  time  stochastic  Gillespie  algorithm. 
The  motivation  for  applying  the  latter  is  to  determine  where 
stochastic  effects  from  the  relatively  rare  processes  of  large 
length  polymer  generation  may  be  observable  in  experiment. 

In  considering  the  relevance  of  these  studies  to  mamma¬ 
lian  prions,  there  is  an  important  caveat  that  arises  from  the 
fact  that  mammalian  prions  are  glycosyl  phosphatidyl- 
inositol  (GPI)  anchored  on  neuron  cells  and,  thus,  move  in  a 
two-dimensional  plane  and,  hence,  can  lead  to  areal  aggre¬ 
gates  [13,14].  The  aggregation  and  fissioning  processes  are 
likely  to  be  very  different  there.  Nevertheless,  the  dominance 
of  fissioning  can  lead  to  long-time  growth  kinetics,  which  are 
in  many  ways  closely  related. 
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FIG.  1.  Aggregation  model  where  circles  are  normal  proteins 
(psi-)  and  squares  are  the  misfolded  type  [PSI+], 

II.  MODELS 

We  consider  a  one-dimensional  aggregation  model  that 
includes  a  nucleation  barrier,  elongation  by  monomer  addi¬ 
tion,  and  bssion  (Fig.  1).  The  same  kind  of  kinetics  was 
considered  by  Collins  et  al.  in  the  form  of  moments  of  the 
full  polymer  kinetic  equations,  but  not  in  the  explicit  length 
resolved  form  considered  here  [9,10].  The  assumption  that 
monomer  (rather  than  oligomer)  addition  dominates  and  the 
need  for  fission  are  underscored  by  the  data  we  seek  to  de¬ 
scribe  on  yeast  prions  [9,10].  We  shall  distinguish  the  model 
quantitatively  from  a  nonfission  based  model  later  in  the 
paper. 

A  nucleation  event  occurs  very  rarely  and  is  composed  of 
two  normal  [psi-]  proteins  interacting  in  such  a  way  that  a 
misfold  occurs  and  they  bind  together;  this  two  prion  aggre¬ 
gate  is  called  a  [PSI+]  dimer.  The  actual  size  of  the  nucleus  is 
not  a  critical  parameter  in  our  model,  as  long  as  the  nucle¬ 
ation  process  is  much  slower  than  aggregation.  Once  an  ag¬ 
gregate  has  been  created,  it  can  elongate  by  monomer  addi¬ 
tion.  Our  model  does  not  consider  fusion  between  two 
oligomers.  The  rate  of  fissioning  becomes  important  only 
when  an  oligomer  of  large  enough  size  arises.  After  an  oli¬ 
gomer  has  fissioned  into  smaller  pieces,  the  individual  pieces 
can  grow  and  fission.  When  a  fission  leads  to  a  monomeric 
product,  we  assume  that  it  returns  to  the  properly  folded 
state.  This  process  leads  to  a  steady  state  distribution  of 
oligomer  sizes  and  a  pure  exponential  growth.  It  continues 
until  the  monomer  concentration  begins  to  dwindle 
significantly. 

We  note  that  the  lag  time  for  the  in  vitro  growth  of  yeast 
prion-like  proteins  is  reproducible  in  experiments,  and  shows 
a  weak  dependence  on  monomer  concentration.  Hence,  it  is 
not  likely  that  it  can  be  associated  with  pure  nucleation  (at 
least  at  low  initial  monomer  concentration),  which  being  a 
rare  event  will  lead  to  a  distribution  of  time  scales  as  well  as 
a  rather  strong  dependence  on  monomer  concentration,  de¬ 
pending  on  the  size  of  the  nucleus.  We  will  assume  that  the 
lag  time  arises  from  the  growth  of  misfolded  material  to  an 
observability  threshold.  This  assumption  is  valid  if  the  nucle¬ 
ation  time  scales,  which  are  long  compared  to  other  micro¬ 


scopic  time  scales,  are  still  shorter  than  this  latter  time  scale. 
It  is  also  valid  when  a  small  number  of  nucleated  seeds  exist 
at  the  time  of  the  preparation  of  the  experimental  samples. 
As  shown  elsewhere,  a  weak  concentration  dependence  to 
lag  time  can  derive  from  rapid  equilibration  of  monomer 
concentration  with  nonfibrillar  oligomers  [11,12];  this  model 
is  not  supported  in  the  current  context,  and  we  shall  discuss 
this  later  in  this  section. 


A.  Kinetic  model 

Our  first  approach  is  to  use  kinetic  equations  to  describe 
the  aggregation  process  [15,16].  In  order  to  use  the  rate  equa¬ 
tions,  one  must  assume  a  large  number  of  monomers.  The 
kinetic  rate  equations  for  our  model  are  as  follows 


dnx 

dt 


=  -  2  p, 


nucn  1 


N-\  N-\ 

2p+«l2  ni  +  2pf^j  nh 
i=2  i= 2 


(1) 


dn2  2 
~  Pnuc11 1 


N- 1 

2p+n  |  n2  -  pfi2  +  2pf^  nh 
i'=3 


(2) 


dnk 

dt 


N- 1 

=  2p+nk_tn]  -  2 p+nknx  -{k-  1  )pfik  +  2pf  2)  «„ 

i-k+2 


(3) 


dnN 

—  =2p+nN_]ni  - pf(N  —  \)nN_u  (4) 

where  nx  is  the  monomer  concentration,  nk  is  the  concentra¬ 
tion  of  A'-mers,  and  N  is  the  longest  oligomer  kept  in  the 
calculation.  The  parameters  pnuc,  p+,  and  pj-  can  be  adjusted 
to  change  the  rate  of  nucleation,  fusion,  and  fission,  respec¬ 
tively.  Since  one  of  these  parameters  can  be  absorbed  into 
the  definition  of  the  time  f,  we  will  set  p+=  1.  Furthermore, 
we  will  set  the  initial  monomer  concentration  to  unity.  The 
equations  listed  above  form  a  system  of  A'-coupled  differen¬ 
tial  equations  which  cannot  be  solved  analytically  because  of 
the  nonlinear  terms. 

An  analytical  solution  is  possible  if  nx  can  be  assumed  to 
be  a  constant  and  n2  and  higher  are  much,  much  smaller  than 
the  monomer  concentration.  Then  the  set  of  equations  be¬ 
come  linear  and  can  be  solved  by  several  techniques,  such  as 
a  Laplace  transformation.  This  is  not  a  bad  approximation, 
for,  as  we  will  see,  during  much  of  the  growth  process  the 
monomer  concentration  is  nearly  constant.  Only  at  the  end, 
the  monomer  concentration  begins  to  dwindle  and  the 
amount  of  misfolded  monomers  saturates.  It  is  useful  to  de¬ 
fine  the  zeroth  and  first  moments  of  the  aggregate  size  dis¬ 
tribution  as 


N 

A  =  2  nh 

i=2 


(5) 
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M  =  2  *  *  nn  (6) 

1=2 

where  A  is  the  number  of  polymer  aggregates  and  M  is  the 
total  number  of  monomers  in  aggregate.  Since  the  overall 
number  of  proteins  is  conserved,  the  rate  of  change  of  M  is 
given  by 

dM  dn  | 

- =  -— 7 

dt  dt 

whereas  the  rate  of  change  of  A  is  given  (ignoring  the  cutoff 
N)  by  the  equation 

dA  -)  ,  ■ 

—  =Pnucn\  -Pfll+Pflj  (;-3)«/.  (8) 

i= 3 


If  we  ignore  the  terms  depending  on  n2  and  assume  « ,  is  a 
constant,  we  get  two  linear  coupled  equations, 

dA 

—  =  const  +  pfM  -3pfA,  (9) 


dM 

- =  const  +  (2  p+nx 

dt 


-2 Pf)A. 


(10) 


We  set  p+=  1  to  set  the  unit  of  time.  Anticipating  M  > A , 
which  means  that  the  mean  aggregate  size  is  much  bigger 
than  unity.  The  above  equations  mean  that  the  M  and  A  both 
grow  exponentially  as  ekl,  with  X  given  by 

X  =  V 2pf{nx  -pf)~  i2pfn x .  (11) 

During  the  exponential  growth,  the  mean  aggregate  size  is 
given  by 


If  we  further  assume  that  the  lag  time  is  given  by  the  time  it 
takes  for  the  misfolded  proteins  to  reach  a  detectable  thresh¬ 
old  Mj-,  then  the  lag  time  is  given  by 


In 


Mr 


‘lag  ' 


C(n  i)2 
V2  pfnl 


(13) 


where  C  is  a  constant  and  the  n  x  dependence  inside  the  loga¬ 
rithm  comes  from  the  prefactors.  We  will  see  that  these  re¬ 
sults  are  confirmed  by  a  complete  numerical  integration  of 
the  differential  equations.  Thus,  this  model  demonstrates 
the  requisite  /-lag  dependence  seen  in  the  experimental 
data,  namely  the  inverse  square  root  dependence  on  the 
monomer  concentration,  up  to  logarithms,  shown  in  Fig.  2. 
We  note  that  the  fit  in  Fig.  2  is  not  sensitive  to  the  logarith¬ 
mic  terms. 


Kinetic  model  results 

A  numerical  integration  approach  was  used  to  study  the 
full  set  of  equations.  We  used  the  fourth  order  Runge-Kutta 
algorithm  for  solving  differential  equations  and  to  obtain 


FIG.  2.  Lag  time  ( ttag )  vs  initial  monomer  concentration  [/? x (0)]. 
The  crosses  are  data  from  Collins  et  al.  in  Ref.  [9].  The  continous 
curve  is  a  fit  to  the  equation  f/a£=«o[oi(0)]_I/2  log  axhix(t)),  with 
ci0=9.22  min* (pM)m  and  fl[  =  1.80X  10 6pM,  following  Eq.  (13). 

length  distributions  and  growth  curves  for  the  system. 

In  Fig.  3,  we  show  plots  of  log(M)  as  a  function  of  time. 
We  have  chosen  a  set  of  parameters  yielding  mean  lengths  in 
the  10.?- 100s  and  displaying  manifest  exponential  growth  in 
M{t)  vs  t  plots.  It  is  evident  that  there  is  a  regime  where 
log (M)  varies  linearly  with  time,  implying  a  pure  exponen¬ 
tial  growth.  This  region  of  exponential  growth  is  not  limited 
to  just  a  few  chosen  parameter  values,  but  can  be  achieved 
over  a  wide  range  of  pf  values.  In  Fig.  3,  we  plot  the  total 
number  in  aggregate  vs  time,  now  on  a  linear  scale.  The 
sigmoidal  growth  is  now  evident  as  the  exponential  increase 
tapers  off  when  the  monomers  begin  to  deplete.  This  is  also 
consistent  with  the  experimental  results  [9]. 

As  discussed  earlier,  the  mean  length  of  the  aggregates 
has  a  simple  inverse  square-root  dependence  on  the  rate  of 
fission.  Thus  by  tuning  the  rate  of  fission,  it  is  possible  to 
obtain  mean  aggregate  lengths  (see  Fig.  4),  which  appear 


t>0 

O 


FIG.  3.  Log[M(t)]  vs  time,  with  parameter  values  ni(0)  =  l,  p+ 
=  1,  p,mc=  10-8.  The  linear  regions  illustrate  regions  of  pure  expo¬ 
nential  growth.  M(t)  vs  time  illustrating  sigmoidal  growth  is  shown 
in  the  inset.  Note  that  time  is  in  units  of  l/fn^Ojp+j. 
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FIG.  4.  Length  distribution,  nt  vs  length  L, 
from  the  kinetic  equations  during  the  exponential 
growth  phase,  compared  with  the  Weibull  distri¬ 
bution,  Eq.  (14). 


qualitatively  consistent  with  those  observed  experimentally 
displaying  lengths  in  the  tens  of  hundreds  (see  especially 
slide  9  at  Ref.  [10]).  In  Fig.  5,  this  mean  length  is  plotted  vs 
time  for  a  different  fission  rate  pf.  We  see  that  the  mean 
aggregate  length  is  constant  during  the  exponential  growth. 
In  Fig.  6,  we  show  the  variation  of  the  exponential  growth 
rate  with  compared  with  the  analytic  calculations.  One  can 
see  that  they  agree  closely.  Note  also  that  pf  values  of  order 
10~4  will  give  the  lengths  seen  in  vitro. 

We  can  also  determine  the  distribution  of  aggregate  sizes 
in  our  calculations.  The  concentration  of  aggregates  of  dif¬ 
ferent  sizes  relative  to  the  initial  concentration  of  monomers 
are  shown  in  Fig.  4  for  a  given  value  of  the  fission  parameter 
(j>f=  4X  1(T4)  at  three  different  times  during  the  exponential 
growth  phase  of  the  simulation.  Clearly,  the  distribution  is 
fairly  stable  during  the  exponential  growth  phase,  further 
strengthening  the  argument  for  a  simple  steady  state.  These 
length  distributions  are  well  fitted  by  a  Weibull  distribution 
given  by  the  relation. 


FIG.  5.  Mean  length  vs  time  for  different  pf,  «j(0)=  1,  p+=  1, 
Pnuc=  10-8.  The  plateau  corresponds  to  the  constant  mean  length 
during  the  exponential  growth  phase.  Note  that  time  is  in  units  of 

l/[«i(0)p+]. 


fix)  = 


Oq 
Cl  i 


(14) 


Apart  from  an  overall  normalization,  the  two  key  parameters 
of  the  distribution  are  the  mean  ai  and  the  power  law  at 
small  x  set  by  ciq  -  1 .  This  distribution  is  found  to  be  ubiqui¬ 
tous  in  nature  [17].  For  example,  a  dropped  piece  of  coal  will 
have  shattered  pieces  that  follow  an  asymmetrical  distribu¬ 
tion  with  more  pieces  on  the  smaller  end.  Our  model  follows 
a  similar  idea  (i.e.,  taking  a  larger  length  and  shattering  it 
into  smaller  lengths)  and  a  similar  distribution.  In  our  simu¬ 
lations,  the  quantity  a0  is  close  to  2  for  a  range  of  parameters 
studied.  During  subsequent  times,  the  peak  of  this  curve 
shifts  to  the  left  as  saturation  occurs  and  the  only  process  left 
is  to  fission.  (The  shift  in  peak  at  later  times  are  not  shown 
here.) 


FIG.  6.  Exponential  growth  rate  \  vs  pt  compared  with  simple 
square-root  dependence  predicted  by  Eq.  (11)  for  different  param¬ 
eter  values,  nj(0)  =  l,  p+=  1,  pnuc=  10  . 
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B.  Comparison  to  nonfissioning  micelle  model 

The  coupling  of  the  observation  of  weak  concentration 
dependence  with  exponential  growth  strongly  supports  a 
model  dominated  by  fissioning  of  aggregates  constructed  by 
monomer  addition.  Indeed,  centrifugation  studies  confirm 
that  monomers  are  the  dominant  species  in  the  in  vitro  me¬ 
dium,  and  examination  of  the  aggregation  kinetics  suggests  a 
critical  nucleus  of  small  size  [9], 

It  is  thus  important  to  contrast  with  other  models  for 
which  a  weak  concentration  dependence  can  be  derived.  No¬ 
tably,  the  work  of  Lomakin  et  al.  [11,12]  considers  kinetics 
with  monomer  growth  of  polymers  as  well  as  micellar  oligo¬ 
mers  composed  of  A4  monomers  of  concentration  nM  which 
may  or  may  not  be  off  pathway.  This  model  may  be  of  rel¬ 
evance  to  other  in  vitro  studies  of  yeast  prion  aggregation 
[18]  for  which  there  is  evidence  of  oligomer  formation  and 
growth  and  a  weak  concentration  dependence  to  the  lag 
time.  If  the  oligomers  are  off  pathway,  then  assuming  an 
equilibration  with  monomers  leads  to  a  buffering  of  the 
monomer  concentration  assuming  an  initially  high  input  con¬ 
centration  and  sufficiently  rapid  equilibration.  There  is  a 
crossover  or  critical  concentration  n*  for  this  dictated  by  the 
equilibrium  constant  KM  between  monomers  and  micelles 
(■ KM=nM/n given  by 

n=(MKM)m-M.  (15) 

Lomakin  et  al.  argue  that  for  a  total  equilibrated  protein  con¬ 
centration  nT  obtained  prior  to  fibril  growth, 


=  nl  +  MnM 

(16) 

-•cr 

(17) 

From  these  equations,  it  is  clear  that  there  is  an  approximate 
two-state  behavior  provided  J\Ai>  1.  For  nT<n*,  then  nT 
—  Hi  and  aggregation  kinetics  can  be  strongly  concentration 
dependent.  However,  for  nT^n*,  then  n ,  ~  n* ,  that  is,  the 
micelles  buffer  the  monomer  concentration.  This  would  lead 
to  a  weak  dependence  of  fibrillar  growth  upon  subsequent 
increase  of  initial  monomer  concentrations. 

There  are  two  routes  by  which  this  model  can  be  applied 
to  the  data  of  Collins  et  al.  [9,10].  First,  if  the  fibrils  are 
formed  from  monomer  addition  to  the  micelles  without  re¬ 
sulting  fission  of  fibrils,  then  one  will  expect  a  quadratic 
growth  in  time  [19,20]  to  the  total  number  of  aggregated 
proteins,  namely 


M{t)~-p2+n\nMt2,  (18) 

~  ^^pln]+M(n*)'~MnMt2,  (19) 

which  is  very  weakly  dependent  upon  «i(0)  provided  «i(0) 
~nT^n*.  but  varies  strongly  as  n1(0)2+A(  in  the  opposite 
limit.  Solving  for  where  M(tlay)  =  M c  characteristic  of  the  ob¬ 
servability  threshold  yields 


i 

,ag  p+nx  V  (m)M 


(20) 


For  ii\{ 0)<§«*,  t!ag  varies  as  nl(0)^1+M,2\  at  odds  with  ex¬ 
periment  [9],  while  for  «i(0)  and  A4  >  1 ,  tIag  is  approxi¬ 
mately  independent  of  /q(0)  [it  will  decrease  weakly  with 
increasing  «i(0)]. 

Second,  if  the  micelles  are  off  pathway,  then  their  princi¬ 
pal  impact  on  the  fibrillar  aggregation  is  on  buffering.  As¬ 
suming  a  critical  fibrillar  nucleus  of  length  p,  which  equili¬ 
brates  with  the  monomers  via  a  constant  Kp=nplnpx ,  the  same 
considerations  of  aggregation  kinetics  as  in  the  preceding 
paragraph  give  a  lag  time  given  by 


1  /  2  Mc 

P+n\  V  Kp{nx)p 


(21) 


For  n1(0)<^n*,  this  varies  as  n,(0)  ,1+A2)  js  at  odds  with 
Collins  et  al.  [9],  while  it  is  approximately  independent  of 
n^O)  for  /i](0) and  p> 2.  For  small  nuclei  (p  —  2-5, 
say),  the  dependence  can  be  relatively  strong  on  7q(0)  even 
in  this  limit. 

Hence,  via  either  on-pathway  or  off-pathway  micelles,  in 
this  model  there  is  a  route  to  weak  dependence  of  lag  time 
(defined  as  the  threshold  for  observation)  upon  initial  mono¬ 
mer  concentration  requiring  nj(0)  and  A4  ,p  >  1 .  How¬ 
ever,  there  are  several  problems  in  connecting  this  to  the 
experiments  of  Collins  et  al.  First,  exponential  growth  is 
unambiguously  observed,  indicating  that  fission  of  fibrils  oc¬ 
curs.  The  micelle  model  of  this  subsection  has  no  fission. 
This  highlights  the  important  role  of  fission  for  the  prion 
phenomenon.  Second,  in  these  experiments,  the  starting  so¬ 
lution  contains  almost  entirely  monomers  and  not  micelles. 
This  implies  that  even  if  micelles  form,  one  is  always  in  the 
limit  of  ni(0)<n",  leading  to  strong  concentration  depen¬ 
dence.  Third,  the  width  of  the  fibrils  obtained  in  this  experi¬ 
ment  is  apparently  monomeric,  which  is  not  readily  compat¬ 
ible  with  growth  by  monomer  addition  to  micelles,  and 
rather  must  be  obtained  by  addition  to  a  critical  fibril 
nucleus.  The  data  suggest  that  this  critical  fibril  nucleus  is  of 
length  p^6,  not  consistent  with  the  assumption  p  >  \ . 

Hence,  we  do  not  believe  that  the  micelle  based  model 
can  explain  the  data  of  Collins  et  al.  [9,10]. 


C.  Stochastic  model 

In  this  section,  we  develop  a  stochastic  treatment  of  the 
model,  similar  to  a  model  proposed  by  Poschel  et  al.  [21]. 
Our  main  motivation  is  to  be  able  to  treat  species  with  small 
numbers,  in  which  case  the  continuous  deterministic  ap¬ 
proach  will  break  down.  As  we  will  see,  the  stochastic  ap¬ 
proach  largely  agrees  with  the  deterministic  model  for  the 
parameter  values  considered  here  for  aggregate  quantities 
[e.g.,  M(t)].  However,  for  distributions  there  can  be  substan¬ 
tial  stochastic  noise  at  short  times.  Another  advantage  of  the 
stochastic  approach  is  that  it  can  be  readily  extended  to  study 
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TABLE  I.  Propensities  of  stochastic  processes. 


Fusion  propensity  for 

monomer  addition 

Fission  propensity  for  monomer 

Higher  order 
fission  term 

propensity 

Chain  length  2 

Chain  length  3 

fusion  term 

Chain  length  2 

Chain  length  3 

gnucNMl-l) 

2gJfiN2 

2g+NjN3 

gfN2 

2g/N3 

two-dimensional  aggregates  as  well  as  the  problem  of  mul¬ 
tiple  prion  “strains,”  in  which  case  the  treatment  of  a  small 
number  of  heteroaggregates  would  be  important  to  model  the 
extent  to  which  strains  breed  true. 

The  Gillespie  algorithm  provides  an  exact  way  to  treat  the 
stochastic  problem  of  chemical  reactions  [22].  In  our  case, 
the  polymers  of  different  lengths  are  the  different  chemical 
species.  The  processes  of  nucleation,  monomer  addition,  and 
fission  are  assumed  to  be  stochastic.  In  other  words,  the 
number  of  polymers  of  different  length  at  time  t,  only  define 
the  propensities  (or  normalized  probability)  for  the  different 
reactions  to  happen  at  that  time.  Once  a  reaction  takes  place, 
the  number  of  polymers  is  altered  and  the  propensities  are 
changed.  The  Gillespie  algorithm  is  a  Monte  Carlo  treatment 
that  deals  with  the  stochastic  process  by  using  a  pair  of  ran¬ 
dom  numbers  at  each  step,  one  to  decide  which  event  will 
occur  next  and  another  to  decide  how  long  will  it  be  until  the 
next  event  takes  place.  This  process  can  be  repeated  to  fol¬ 
low  the  dynamical  behavior  of  the  system. 

Our  model  has  the  corresponding  propensities 


gnuMNx  -  1)  +  2g+N Ni  +  gjZ  (i  -  l)Nj,  (22) 


i=2 


where  gnuc,  g+,  and  gf  are,  respectively,  the  nucleation,  fu¬ 
sion,  and  fission  parameters  in  the  stochastic  model.  In  our 
model,  each  propensity  was  assigned  a  bin  in  an  array.  A 
random  number  generator  decides  which  bin  is  selected. 

Table  I  illustrates  physically  the  propensities  at  a  given 
time.  A  normalized  random  number  generator  selects  which 
element  in  the  array  will  occur  and  the  time  it  took  to  create 
that  event.  After  an  event  occurs,  the  corresponding  propen¬ 
sities  are  updated  and  the  process  is  repeated. 

In  order  to  compare  the  results  achieved  from  the  continu¬ 
ous  to  stochastic  model,  a  mapping  of  the  Gillespie  param¬ 
eters  to  the  kinetic  parameters  is  required.  For  this,  we  need 
to  develop  an  approximate  equation  satisfied  by  the  mono¬ 
mer  concentration  in  the  stochastic  simulations  and  compare 
it  with  the  rate  equations.  In  the  stochastic  case,  one  has  a 
master  equation  that  relates  the  probability  distribution  asso¬ 
ciated  with  a  different  number  of  polymers  at  time  t+dt  to 
those  at  time  t.  We  will  make  a  mean-field  approximation  for 
different  kind  of  k-mers,  namely  ( V; * Nk) = (A,)  *(Nk).  Thus, 
we  can  consider  the  approximate  master  equation,  which 
only  tracks  changes  in  the  monomer  number.  We  arrive  at  the 
equation. 


Pl{Nl,t  +  dt)  =  Pl{Nl,t) 


1  ~  8, mcNiin^  -  l)dt-g+N Njdt 


i=2 


■  2 gf^j  Njdt 


i=2 


+  2P,(A,  -  1  ,t)g£  Njdt 


i= 3 


+  Pt(Nt  -  2,t)gfN2dt  +  H|(A'|  +  1  .rWVj 


x2  Nidt  +  P1(Nl  +  2,t)gnucN1(Ni  -  1  )dt, 

/=2 

(23) 

where  A]  is  the  number  of  monomers,  P\(N\ , t+dt)  is  the 
probability  to  have  Nx  monomers  at  a  later  time  t+dt. 
P|(V|  ,t)  is  the  probability  to  have  N1  monomers  at  time  t.  It 
is  multiplied  by  the  probability  that  no  reaction  occurs  in 
time  dt  that  changes  A).  The  rest  of  the  terms  represent  the 
probability  to  have  a  different  number  of  monomers  at  time  t 
but  then  a  reaction  happens  in  time  dt  leading  to  Ai  mono¬ 
mers  at  time  t+dt. 

Now  using  the  definition  of  the  derivative,  one  can  rear¬ 
range  the  equation  and  multiply  both  sides  by  Aj  to  get  the 
average  rate  of  change  of  Nt .  After  shifting  some  indices  to 
get  every  probability  to  have  the  form  P\(Nj  ,t)  and  using  the 
appropriate  volume  element  (denoted  by  V)  to  normalize  the 
counts  to  a  concentration,  one  arrives  at  the  equation  for  the 
mean  number  of  monomers. 


1  dnx 
Vg+dtgil 


=  ~  2gnucn\  -  2n^  n,  +  2^2 


i=2 


V 


n-j- 


(24) 


By  comparing  with  the  rate  equations  discussed  before,  we 
obtain  the  parameter  mapping 


P+=  1> 

(25) 

Pnuc  ~  Snuct 

(26) 

ii 

(27) 

t  kinetic  8  +  gillespie  • 

(28) 

Results  stochastic  approach 

After  achieving  the  parameter  mapping,  we  studied  the 
stochastic  models  with  106  initial  monomers.  The  numerical 
results  are  shown  in  the  next  few  figures.  In  Fig.  7,  the 
growth  of  the  aggregate  material  is  shown  on  a  logarithmic 
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FIG.  7.  Log[M(f)]  and  M(t),  shown  in  the  inset,  vs  time  ob¬ 
tained  by  the  Gillespie  algorithm  with  IVi(0)  =  10  ,  g+=l,  gnuc 
=  10“8  for  different  fission  rate  gt.  Note  that  time  is  in  units  of  1  /g+. 

and  on  a  linear  scale.  It  is  evident  that  they  are  in  excellent 
agreement  with  the  corresponding  calculations  for  the  kinetic 
model.  In  Fig.  8,  we  show  plots  of  the  mean  length  as  a 
function  of  time,  and  in  Fig.  9,  the  corresponding  distribution 
of  length  scales  during  the  exponential  growth  phase.  It  is 
again  evident  that  there  is  a  steady  state  during  the  exponen¬ 
tial  growth  and  the  mean  lengths  reach  a  plateau  value.  The 
stochastic  effects  are  much  larger  in  the  distribution  and  in 
the  mean  length,  but  the  overall  results  agree  well  with  the 
kinetic  model.  The  fits  to  the  Weibull  distribution  are  again 
good  with  comparable  parameters.  The  differences  between 
the  kinetic  and  stochastic  Weibull  parameters  ( a0  and  «,)  are 
less  than  10%.  This  shows  that  for  106  monomers  and  for 
mean  aggregate  lengths  up  to  a  few  hundred,  the  dominant 
species  still  occur  in  large  enough  numbers  so  that  stochastic 
effects  do  not  change  the  results  in  a  significant  way. 

III.  DISCUSSION 


In  this  paper,  we  have  used  a  deterministic  kinetic  and  a 
stochastic  Monte  Carlo  approach  to  model  the  dynamics  of 


FIG.  8.  Mean  length  vs  time  in  units  of  1  / g+  for  different  fission 
rates  gt  obtained  by  the  Gillespie  algorithm.  The  parameter  values 
are  A1(0)  =  106,  g+=  L  gnUc=  FT8. 


FIG.  9.  Length  distributions  at  different  times  during  the  expo¬ 
nential  growth  phase  obtained  by  the  Gillespie  algorithm  for  g, 
=  400  and  fitted  by  the  Weibull  distribution. 

yeast  prion  growth.  For  the  parameters  relevant  to  in  vitro 
experiments,  where  mean  aggregate  size  is  of  order  lOOs  of 
monomers,  and  the  number  of  monomers  is  larger  than  106, 
the  two  approaches  agree  well  with  each  other,  showing  that 
stochastic  effects  are  not  dominant.  The  studies  confirm  that 
the  key  experimental  features  of  the  growth  of  yeast  prions 
can  be  well  captured  by  these  nucleation-growth-hssioning 
models.  They  lead  to  an  exponential  and/or  sigmoidal 
growth,  an  inverse  square-root  dependence  of  the  lag  time  on 
the  monomer  concentration,  and  aggregate  sizes  that  depend 
on  the  rate  of  fission.  We  predict  that  the  aggregate  size 
distributions  should  be  Weibullian,  reflecting  the  importance 
of  fissioning  in  the  growth  process.  We  also  predict,  for  suf¬ 
ficiently  large  times,  that  the  steady  state  aggregate  length 
will  drop  as  fissioning  converts  large  aggregates  to  those  of 
minimal  size.  We  have  argued  that  models  in  which  mono¬ 
mers  equilibrate  with  micelles  but  do  not  fission  [11,12], 
which  are  capable  of  producing  lag  times  having  a  weak 
concentration  dependence  for  a  sufficiently  high  initial 
monomer  concentration,  are  in  fact  not  appropriate  for  the 
data  of  Collins  et  al.  [9,10]. 

We  note  that  the  in  vitro  aggregation  work  of  Serio  et  al. 
[18]  appears  to  come  to  different  conclusions  than  that  of 
Collins  et  al.  [9]  and  may  be  ripe  for  a  discussion  in  terms  of 
the  micelle  models.  We  do  not  understand  the  discrepancies 
between  these  two  sets  of  experiments. 

We  note  that  the  results  on  yeast  prions  discussed  here 
and  recent  work  on  mammalian  prions  stripped  of  their  GPI 
membrane  anchors  [23,24]  suggest  that  a  fundamentally  dif¬ 
ferent  mechanism  for  fission  is  at  play  for  the  latter  proteins 
vis-a-vis  yeast.  In  the  recent  experiments,  transgenic  mice 
expressing  the  cellular  prion  protein  without  the  residues 
necessary  for  the  GPI  anchor  are  exposed  to  anchorless  in¬ 
fectious  prions;  in  time  course  experiments,  prion  aggregates 
are  produced  which  retain  infectivity.  However,  the  infected 
mice  with  anchorless  PrPc  do  not  show  clinical  symptoms. 
Moreover,  the  time  course  data  [see  Fig.  1(c)  of  Chesebro 
et  al.  [24]}  show  two  remarkable  characteristics:  (i)  they  are 
nonsigmoidal  in  shape,  with  the  infectious  prion  content  at 
long  times  significantly  exceeding  that  of  infected  wild  type 
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mice  at  death,  (ii)  the  time  scales  to  reach  the  levels  of  in- 
fectivity  characteristic  of  wild  type  mice  at  the  correspond¬ 
ing  infectious  dose  are  quite  long.  Indeed,  the  time  course 
data  can  be  fit  roughly  by  a  quadratic  in  time  curve  charac¬ 
teristic  of  nonfissioning  aggregation  to  an  initial  concentra¬ 
tion  of  infectious  seeds.  Evidently,  the  binding  of  cellular 
mammalian  prions  to  the  membrane  is  critical  to  the  fission¬ 
ing  process,  while  for  yeast  prions  fissioning  of  aggregates 
in  vitro  is  observed  [9]. 

In  the  future,  we  hope  to  extend  this  model  to  study  yeast 
prion  strains,  where  we  expect  the  stochastic  treatment  to  be 
the  key  to  dealing  with  rare  heteroaggregates.  An  under¬ 


standing  of  this  process  should  lead  to  insights  into  the  im¬ 
portant  problem  of  strain  dynamics  in  mammalian  prions. 
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Abstract 

Numerous  human  diseases  are  associated  with  conformational  change  and 
aggregation  of  proteins,  including  Alzheimer’s,  Parkinson’s,  prion  diseases  (such  as  mad 
cow  disease),  familial  amyotrophic  lateral  sclerosis  (ALS,  or  Lou  Gehrig’s  disease), 
Huntington’s,  and  type  II  (mature  onset)  diabetes.  In  many  cases,  it  has  been 
demonstrated  that  conformational  change  and  aggregation  can  occur  outside  living  cells 
and  complex  biochemical  networks.  Hence,  approaches  from  materials  and  physical 
science  have  enhanced  our  understanding  of  the  role  of  protein  aggregation  in  these 
diseases  at  the  molecular  and  nanoscale  levels.  In  this  article,  we  will  review  what  is 
known  about  these  protein  structures  from  the  perspective  of  materials  science,  focusing 
on  the  details  of  emergent  oligomeric  and  nanotube-like  structures,  their  interactions  with 
model  lipid  bilayers,  how  the  structures  relate  to  observed  biological  phenomena,  and 
how  protein  aggregation  and  amyloid  formation  can  be  employed  for  the  good  in  biology 
and  materials  science. 

Keywords:  amyloid  diseases,  complex  adaptive  matter,  emergent  behavior,  prion 
diseases,  protein  aggregation,  protofibrils,  nanotubes,  nanowires. 


Amyloid  Diseases:  A  Truly 
Emergent  Phenomenon 

The  primary  goal  of  the  study  of  com¬ 
plex  adaptive  matter  is  to  identify  the  key 
organizing  principles  that  govern  mate¬ 
rials  phenomena  at  given  length  and  time 
scales.  In  the  case  of  the  folding  of  indi¬ 
vidual  globular  proteins  into  their  func¬ 
tional  structures,  we  know  a  great  deal. 
For  example,  the  local  rules  governing  a- 
helical  secondary  structure  have  been  un¬ 
derstood  in  detail  for  some  40  years;1,2 
indeed,  the  a-helix  was  predicted  by  Paul¬ 
ing  and  collaborators  prior  to  discovery.3 
For  the  global  tertiary  structure  of  the  pro¬ 
tein,  extensive  experiments  and  simula¬ 
tion  studies  show  that  proteins  engineered 
by  evolution  to  experience  minimal  frus¬ 
tration  in  the  interactions  between  closely 
contacted  amino  acids  (or  residues)  de¬ 
velop  funneled  energy  landscapes  and  rel¬ 
atively  rapid  (and  multiple)  pathways  to 
folding  after  synthesis.4  As  discussed  in 
the  article  by  Ramirez  in  this  issue,  "frus¬ 
tration"  refers  to  the  effect  of  competing 
interactions  that  make  it  impossible  to  fa¬ 
vorably  lower  the  interaction  energy. 


However,  we  have  comparatively  little 
understanding  of  the  organizing  princi¬ 
ples  governing  structure  formation  for 
proteins  interacting  with  other  proteins  or 
membranes.  This  lack  of  knowledge  is 
problematic,  because  proteins  left  alone 
tend  to  spontaneously  aggregate,  often  by 
formation  of  (3-sheet  structures,  which  are 
not  governed  by  local  formation  rules  like 
a-helices  (the  distances  between  hydro¬ 
gen  bonding  residues  along  the  backbone 
can  be  great).  |3-sheets  are  formed  from 
approximately  linear  stretches  of  peptide 
that  hydrogen-bond  from  line  to  line. 
These  structures  are  especially  prone  to 
protein  aggregation  due  to  favorable 
edge-to-edge  hydrogen  bonding  between 
sheets.5  This  aggregation  tendency  is  an 
obstacle  in  high-throughput  proteomics, 
where  one  is  interested  in  measuring  the 
properties  of  individual  proteins.6 

Purposeful  biological  aggregation  of 
monomeric  proteins,  as  in  the  assembly  of 
actin  filaments  or  microtubules,7  is  usually 
highly  regulated  and  energetically  con¬ 


trolled  (we  will  discuss  tightly  controlled 
and  biologically  useful  (3-sheet  self- 
assembly  later  in  the  article).  Living 
organisms  have  evolved  an  effective 
quality-control  system  to  prevent  protein 
misfolding  and  aggregation,  where  chap¬ 
erone  proteins  provide  "safe  houses"  for 
folding  proteins,8  and  the  ubiquitin/pro- 
teosome  system  ensures  rapid  degrada¬ 
tion  or  disposal  of  misfolded  proteins. 

Table  I  conveys  the  tragic  side  of  uncon¬ 
trolled  (3-sheet  self-assembly:  it  summa¬ 
rizes  the  key  aspects  of  seven  (out  of 
dozens)  prominent  human  amyloid  dis¬ 
eases.  Amyloid  means  "starch-like" — the 
aggregates  stain  like  starch.  These  dis¬ 
eases  typically  arise  in  old-  or  middle- 
aged  populations,  and  frequently  arise 
spontaneously  or  sporadically  rather  than 
from  genetic  predisposition.  Indeed,  for 
spontaneous  Alzheimer's,  Parkinson's, 
prion,  and  immunoglobulin  light  chain 
diseases,  incidence  varies  little  between 
countries.  Moreover,  prion  diseases,  the 
lone  infectious  type  of  amyloid  disease, 
show  highly  reproducible  dose-versus- 
incubation-time  distributions  for  inter- 
cerebrally  inoculated  animals.9  Remarkably, 
infectious  protein-only  prion  aggregates 
have  been  grown  in  vitro, 10  proving  that 
the  protein  interaction  properties  alone, 
without  additional  biochemical  guidance, 
account  for  the  disease. 

These  observations  suggest  that  these 
diseases  can  be  studied  from  the  perspec¬ 
tive  of  materials  growth,  without  exten¬ 
sive  biological  modulation.  Indeed,  it 
appears  that  the  growth  of  amyloid  struc¬ 
ture  has  much  in  common  with  the  ori¬ 
ented  aggregation  of  inorganic,  nearly 
monodispersed  nanoparticles.11  We  seek 
here  to  portray  the  growing  scientific 
movement  toward  the  use  of  concepts  and 
tools  from  materials  science  in  the  study  of 
amyloidogenic  proteins  to  elucidate  the 
mechanisms  of  disease  and  design  new 
materials. 

Amyloid  Structures:  Plaques, 
Protein  Nanotubes,  and  Oligomers 
Plaques 

The  extracellular  and/or  intracellular 
accumulation  of  amyloid  fibrils  in  the 
form  of  plaques  or  inclusions  (Lewy  bod¬ 
ies)  in  the  brain  is  a  defining  hallmark  for 
several  neurodegenerative  diseases  (Fig¬ 
ure  1).  For  example,  Alzheimer's  disease 
patients  have  large  quantities  of  postmortem 
brain  plaques,  predominantly  composed 
of  40-42-amino-acid-long  Ap  peptides 
that  are  cleaved  by  protease  proteins  from 
the  Alzheimer's  precursor  protein  (BAPP).13 
These  micrometer-scale  plaques  are  com¬ 
posed  of  multi-polymeric  strands  of  the 
peptide,  called  fibrils  (Figure  le),  which 
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Table  I:  Aspects  of  Some  Human  Amyloid  Diseases  with  Associated  Protein  Aggregation. 


Disease 

Alzheimer’s51 

Protein/Peptide 

p-42  (from  BAPP)  and  x 

Function 

? 

Heritable 

Component 

~25% 

Incidence 

—50%  of  post-85-yr-old 
population 

Onset  Age* 
(Years) 

>65 

Parkinson’s52 

a-synuclein 

? 

5-10% 

-1%  of  post-50-yr-old 
population 

55-60  yrs 

Huntington’s23 

huntingtin 

? 

100% 

1  in  20,000  (Caucasian) 

35-40 

Familial  ALS 
(Lou  Gehrig’s 
disease)53 

Superoxide  dismutase 
(20%  of  cases) 

Lowers  oxidative 

stress 

100%  (5-10%  of 
all  ALS) 

2  in  106 

46 

Type  II  diabetes54 

IAPP 

? 

High  (obesity 
trigger) 

14  x  106per  year  (U.S.) 

>40 

Immunoglobulin 
light  chain55 

IG  light  chain 

Immune  response 

Small/unknown 

1  in  1 0s 

64 

Prion  diseases35 

Prpc 

Lowers  oxidative 

stress? 

10-15% 

-1  in  106 

63 

*lncidence/onset  age  is  for  non-heritable  sporadic  disease,  unless  otherwise  noted. 


have  high  quantifies  of  P-sheet  structure 
as  revealed  by  crystallography,  circular 
dichroism  (which  detects  the  different 
light  polarization  rotation  tendencies  of  a- 
helices  and  P-sheets),  and  more  recently, 
solid-state  magnetic  resonance  studies. 
These  P-strands  are  aligned  perpendicular 
to  the  fibrillar  axis,  in  a  so-called  "cross-P" 
structure,  shown  in  Figure  If.  Postmortem 
plaques  and  inclusions  from  a  variety  of 


diseases  are  shown  in  Figure  l.13  We  note 
that  the  amyloid  fibrils  are  protein  nano¬ 
tubes,  hollow  in  the  middle,  with  diame¬ 
ters  of  the  order  of  10-20  nm. 

Despite  a  dominant  research  focus  on 
both  structure  measurements  and  com¬ 
puter  modeling  on  plaques,  an  emerging 
perspective  is  that  the  plaques  may  repre¬ 
sent  disease  end  points  having  little  to  do 
with  toxicity.14  This  view  is  supported  by 


observations  such  as  (1)  the  abundant  AP 
plaques  observed  in  the  brains  of  individ¬ 
uals  displaying  no  symptoms  of  Alzheimer's 
disease  and  (2)  the  non-uniformity  of 
plaque  observation  in  prion  diseases:  for  ex¬ 
ample,  victims  of  kuru  (a  disease  among 
the  Fore  people  of  New  Guinea,  arising 
from  ritual  cannibalism  of  deceased  tribe 
members)  exhibit  them,  while  victims  of 
spontaneous  Creutzfeldt-Jakob  disease 
usually  do  not.15  In  this  view,  toxicity 
is  engendered  by  small  p-sheet  aggre¬ 
gates,  possibly  on  the  pathway  to  amyloid 
fibrils. 

P-, Helices 

Accordingly,  attention  has  turned  to 
small  P-sheet  motifs  with  multiple  assem¬ 
bly  outcomes,  especially  the  left-handed 
p-helix  (LHBH)  structures  shown  in  Fig¬ 
ure  2, 16  recently  proposed  as  the  P-sheet 
unit  for  infectious  mammalian  prion 
trimers  on  the  basis  of  cryogenic  electron 
microscopic  data.17  The  LHBH  P-sheet 
structure  has  also  been  proposed  for  yeast 
prion-like  proteins,18  Alzheimer's  disease,19 
and  Huntington's  disease.16'20  This  LHBH 
motif  was  first  observed  in  several  bacter¬ 
ial  enzymes  and  the  "antifreeze"  protein 
of  the  spruce  budworm;  to  date,  there  are 
11  structures  in  the  protein  data  bank 
(http://www.rcsb.org/pdb/)  confirmed 
to  have  LHBHs.  LHBHs  are  usually  pre¬ 
sumed  (or,  in  two  cases,  observed)21'22  to  be 
in  protein  trimers.  The  LHBH  motif  has  a 
fundamental  repeat  unit  of  triangular 
cross  section,  consisting  of  18  amino  acids 
with  two  per  bend  region  and  four  per 


Figure  1.  Plaques  and  fibrils,  (a)-(d)  Postmortem  tissue  plaques  from  human  amyloid 
disease  patients  (from  Reference  13):  (a)  Ap42  (Alzheimer's)  plaque,  (b)  a-synuclein 
(Parkinson's)  plaque,  (c)  huntingtin  plaque  (Huntington's  disease),  and  (d)  PrPSc  (kuru) 
plaque,  (e)  Transmission  electron  microscope  images  of  amyloid  fibrils  (H.  Lashuel, 
unpublished  data),  (f)  Model  of  a  hollow-core  SH3  domain  fibril:  upper  panel  shows  a 
density  map  cross  section,  lower  panel  shows  the  cross-p  structure.56 
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Figure  2.  (a)  Peptide  backbone  of  measured  left-handed  p-helix  (LHBH)  structure  of  1 FWY 
protein,  (b)  Proposed  LFiBFi  structure  of  infectious  prion  protein  (after  Reference  16). 
Parallel,  hydrogen-bonded  p-sheets  are  shown  as  strips,  (c)  Superposed  structure  of  six 
LHBFI  backbones,  illustrating  the  uniform  cross  section  of  the  helical  structure  (red,  1  FWY 
protein;  blue,  proposed  infectious  prion  protein;  yellow,  1 HMO  protein;  white,  1J2Z  protein; 
green,  1L0S  protein;  fuchsia,  1 T3D  protein),  (d)  LHBFI  “nanotube”  formed  from  a  stack  of 
three  1FWY  helices.  All  images  generated  with  the  RASTOP  molecular  viewing  program 
(see  www.geneinfinity.org/rastop/);  all  files  for  the  proteins  1FWY,  1HMO,  1J2Z,  1  LOS,  and 
1 T3D  input  to  RASTOP  are  from  the  Protein  Data  Bank  (http://www.rcsb.org),  except  for 
the  PrPSc  model,  which  is  courtesy  of  C.  Govaerts. 


(3-strand,  alternating  between  hydropho¬ 
bic  and  hydrophilic  in  the  strand.  The 
triangular  cross  section  is  very  uniform 
across  observed  LHBHs  (Figure  2c); 
each  edge  is  ~1.9  nm  in  length.  In  Fig¬ 
ure  2d,  we  show  how  multiple  copies 
of  the  LHBH  of  an  enzyme  from  E.  coli 
(labeled  1FWY  on  the  protein  data  base 
at  www.rcsb.org/pdb)  can  be  assembled 
into  a  kind  of  (3-nanotube.  This  raises 
the  natural  and  intriguing  question: 
are  amyloid  fibrils  composed  of  such 
"nanofilaments"?19 

The  LHBH  has  a  tantalizing  connection 
to  Huntington's  and  other  polyglutamine 
repeat  diseases,  where  pathology  derives 
from  an  inherited  excess  number  of  re¬ 
peats  of  the  amino  acid  glutamine  on  one 
end  (the  N-terminus)  of  the  huntingtin 
protein.23  A  polyglutamine  repeat  number 
p  of  less  than  ~24  is  normal;  for  p  >  36, 
disease  is  certain.  Note  that  36  is  the  num¬ 
ber  of  amino  acids  or  residues  in  two  turns 
of  a  LHBH  that  has  fully  saturated  inter¬ 
nal  hydrogen  bonds.19 

The  LHBH  presents  a  challenge  to 
theory  and  simulation.  At  present,  this 
motif  has  not  emerged  from  any  molecu¬ 
lar  dynamics  simulations  or  semi-analytic 
(Hamiltonian  or  cellular  automata)  ap¬ 


proaches.  The  complicating  factor  is  the 
long-range  coupling  along  the  backbone 
(bonded  amino  acids  are  separated  by  18 
residues).  On  the  other  hand,  the  remark¬ 
able  conservation  of  the  shape  and  helical 
cross  section  suggest  that  this  is  a  motif 
ripe  for  study. 

Nanoscale  Oligomers 

(3-sheet  oligomers  obtained  from  in  vitro 
growth  of  aggregates  have  been  studied 
extensively.13  One  common  form  seems  to 
be  a  spherical  micelle-like  aggregate  of  (3- 
converted  proteins,  which  have  been  im¬ 
plicated  as  precursors  for  the  chain-like 
(protofibril)  and  annular  oligomers  that 
are  also  seen  frequently  during  in  vitro 
amyloid  formation  by  most  amyloido- 
genic  proteins.24  Figure  3  shows  a  com¬ 
pendium  of  oligomers  and  protofibrils 
from  amyloid  disease  proteins.25,26  The  an¬ 
nular  oligomers  underlie  the  proposed 
toxicity  mechanism  discussed  later  in  this 
article.  We  note  that  all  atomic  force  micro¬ 
scope  (AFM)  measurements  of  spherical, 
chain-like,  and  annular  oligomers  share 
the  same  height  and  diameter,  suggesting 
that  the  spherical  aggregates  are  the  pre¬ 
cursors  to  the  larger  chain  and  annular 
pore-like  structures. 


The  diversity  of  pre-fibrillar  oligomer 
structures,  the  formation  of  which  are 
highly  dependent  upon  protein  sequence 
and  environmental  conditions  (e.g.,  pH, 
salt  concentration,  and  levels  of  molecular 
crowding),  may  explain  the  mystery  of 
prion  disease  strains.  Strains,  for  a  given 
mammal,  have  unique  incubation-time 
versus  dose  distributions,  tissue  lesion  pro¬ 
files,  and  distributions  of  post-translationally 
attached  sugars:  the  prion  protein  can 
have  0,  1,  or  2  sugars  attached.27,28  More¬ 
over,  strains  breed  true  upon  multiple 
passage  in  animals.  This  means  that  upon 
passage  from  a  diseased  animal  to  a 
healthy  animal,  the  same  properties  (incu¬ 
bation  time,  lesion  profile,  sugar  binding) 
are  preserved.  There  is  considerable  evi¬ 
dence  that  strain  is  encoded  in  prion  con¬ 
formation  but  no  detailed  understanding 
of  the  underlying  mechanisms.27,28  Prion 
oligomers17  might  have  a  spectrum  of  dif¬ 
ferent  shapes,  and  oriented  aggregation  of 
such  protein  "nanoparticles"  might  "breed" 
the  conformation  true  (nonmatching  shapes 
will  be  energetically  unfavorable). 

Domain  Swapping 

A  separate  way  to  generate  fibrils  and 
oligomers  is  by  "domain  swapping,"29,30 
in  which,  say,  two  identical  copies  (A,B)  of 
a  protein  exchange  a  domain  (a  well-defined 
protein  region  attached  to  a  flexible  sec¬ 
tion  of  peptide).  The  domain  of  monomer 
A  binds  to  the  corresponding  region  of  B 
and  vice  versa.  The  swapping  is  not  lim¬ 
ited  to  dimers;  it  can  lead  to  filamentary 
structures  in  which  the  zth  protein  swaps 
with  the  (i  +  l)th  protein,  for  example, 
closed-chain  oligomers  (where  filament 
ends  are  brought  together  to  domain- 
swap)  or  two-dimensional  structures. 

Domain  swapping  has  been  proposed 
to  play  a  role  in  prion  strains.31  Domain 
swapping  in  protein-protein  interactions 
generically  and  amyloid  diseases  in  par¬ 
ticular  seems  certain  to  emerge  as  a  critical 
theme  in  the  coming  years.  At  the  concep¬ 
tual  level,  the  study  of  coarse-grained  pro¬ 
tein  models  with  molecular  dynamics  has 
shown  that  the  formation  of  domain- 
swapped  dimers  will  proceed  down  a 
funneled  landscape  if  the  dimer  enjoys 
minimal  frustration. 

Aggregation  Pathways  and 
Kinetics 

AFM  and  kinetic  modeling,  well  known 
to  materials  scientists,  have  been  instru¬ 
mental  in  advancing  our  understanding 
of  the  structural  properties  of  the  protein 
aggregates  linked  to  disease  and  their 
growth  kinetics.  Figure  4  schematically 
shows  accepted  models  on  amyloid  con¬ 
version/aggregation  kinetics.  We  note 
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Figure  3.  Amyloid  oligomers,  (a)  Afi  peptide  implicated  in  Alzheimer’s  disease,  observed  by 
electron  microscopy  (scale,  50  nm/box)  (from  Reference  24).  (b)  Nine  different  oligomers  of 
mutant  a-synuclein,  implicated  in  Parkinson's  disease,  observed  by  electron  microscopy 
(scale,  40.5  nm/box)  (from  Reference  25).  (c)  Six  oligomers  of  mutant  S0D1,  implicated  in 
familial  ALS  (from  Reference  56).57  (d)  Nine  mutant  a-synuclein  oligomers  (scale,  40.5  nm/box) 
(from  Reference  25).  (e)  Atomic  force  microscope  (AFM)  images  of  Afi  oligomers  inserted 
in  a  supported  lipid  bilayer  (from  Reference  26).  Clearly  resolved  annular  oligomers  have 
outer  diameters  of  —16  nm.  (f)  AFM  image  of  a-synuclein  oligomer  on  supported  lipid 
bilayer  (from  Reference  25).  (g)  Schematic  illustration  of  oligomerization  pathways.  Flere, 
the  subscript  n  in  the  second  stage  represents  the  number  of  monomers  in  an  oligomer. 


that  P-sheet  conversion  is  typically  con¬ 
comitant  with  aggregation  (i.e.,  conforma¬ 
tional  changes  and  protein  assembly  are 
linked),32-33  with  the  possible  exception  of 
polyglutamine  proteins,34  and  requires 
templating  either  by  spontaneously  formed 
(and  rare)  nuclei  or  by  external  seeding  of 
aggregates.33  Also,  fission  of  mammalian 
prions  is  necessary  for  disease  propagation 
and  exponentially  growing  aggregates  in 
mammals;  prion-like  proteins  in  yeast 
have  fission  effected  by  another  protein35 
(although  the  fission  can  occur  sponta¬ 
neously  in  vitro).36  Why  mammalian  prions 
fission  and  other  amyloidogenic  proteins 
do  not  remains  a  mystery. 

In  vitro  and  theoretical  studies  of  con¬ 
version/aggregation  suggest  that  the  late- 


age  onset  in  amyloid  diseases  derives  from 
slow  underlying  molecular  processes.  For 
example,  in  vitro  kinetics  experiments  for 
polyglutamine  peptides  extrapolated  to  in 
vivo  concentrations  of  huntingtin  protein 
suggest  that  within  the  "sampling  win¬ 
dow"  of  a  human  lifespan,  toxic  aggregate 
concentrations  should  arise  only  for  gluta¬ 
mine  number  p  >  36,  consistent  with  clin¬ 
ical  observations.34 

For  prion  diseases,  theoretical  modeling 
of  two-dimensional  aggregation  and  fis¬ 
sion  (prions  live  mostly  on  neuronal 
membranes)  yielded  a  sporadic  incuba¬ 
tion  time  distribution  that  peaked  at  ~  100 
times  that  obtained  from  dilute  seeding 
for  physiological  concentrations  of  the 
normally  expressed  protein  called  PrPc 


that  misfolds  in  the  diseased  form.37 
Given  a  mean  incubation  time  for 
kuru  of  12  years,38  this  suggests  that  en¬ 
demic  sporadic  prion  disease  requires 
~  1000-year  life  spans!  Meanwhile,  the  1  in 
106  sporadic  disease  background  inci¬ 
dence  may  reflect  the  low-amplitude, 
pre-peak  tail  in  the  incubation  time 
distribution.38 

Amyloid-Membrane  Interaction 
and  Toxicity 

Many  amyloidogenic  proteins  associate 
with  lipid  membranes.  AFM  studies  on 
supported  bilayers  and  molecular  model¬ 
ing  have  helped  support  a  potential  unify¬ 
ing  hypothesis  for  amyloid  disease 
toxicity:  that  small  oligomers  pierce  cell 
membranes,  triggering  cell  death  through 
superfluous  ion  pores.  It  was  proposed  in 
detail  that  a  pair  of  P-coupled  annular 
tetramers  of  Ap  peptides  can  insert  into 
the  neuronal  membrane  leaflets  and  create 
an  ion  pore.39  In  vitro- grown  oligomers 
qualitatively  consistent  with  this  hypothe¬ 
sis  have  been  found  (Figure  3)  and  studied 
with  AFM  on  supported  bilayers  (Fig¬ 
ure  3f);  these  are  evidently  composed  of 
4  nm  spherical  oligomers,  much  larger 
than  the  peptide  tetramers  of  Refer¬ 
ence  39.  These  oligomers  permit  excess 
calcium  flux,  which  is  toxic  to  cultured 
neurons.26  Intracerebral  inoculation  of  rats 
and  mice  with  a  solution  presumed  rich  in 
oligomeric  Ap  particles  leads  to  a  re¬ 
versible  short-term  memory  deficit.40  A 
recent  theoretical  study  found  strong  cor¬ 
relations  between  the  membrane  insertion 
configuration  and  the  pore  model  for  4 
of  5  mutations  leading  to  early-onset 
Alzheimer's  disease.41  While  the  pore 
model  is  not  universally  accepted,  and 
may  not  apply  to  all  the  diseases  (the 
prion  trimer  model,  for  example,  will  not 
allow  ion  passage),  it  remains  a  vibrant 
area  of  research. 

Amyloids  for  Good  in  Biology  and 
Materials  Science 
Biologically  Useful  Amyloids 

Emerging  evidence  suggests  that  amy¬ 
loid  structures  can  provide  useful  biologi¬ 
cal  functions.  Some  examples  are 

1.  Heritable  amyloid  structure  in  yeast.  As 
alluded  to  previously,  prion-like  proteins 
in  yeast  form  aggregates  that  fission 
upon  cell  division,  can  actively  confer 
phenotype,  and  may  provide  some  stress 
protection.42 

2.  Spider  silk.  Spiders  produce  insoluble  fil¬ 
aments  of  fibroin  protein  that  possess  30x 
greater  extensibility  and  toughness  than 
steel.  Recent  circular  dichroism  studies 
show  that  significant  amyloid-like  cross-P 
structure  develops  in  a  region  of  reduced 
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Figure  4.  Schematic  illustration  of  amyloid  conversion  and  aggregation.  Native  monomers  of 
the  relevant  peptides  or  proteins  are  shown  as  circles  and  high-p-content  converted 
proteins  as  squares.  Peptides  can  either  spontaneously  convert  or  be  seeded.  Subsequent 
aggregation  and  conversion  can  generate  either  proto-fibrils  via  the  elongation  step,  which 
then  form  amyloid  fibrils  and  plaques,  or  pass  to  “off-pathway”  oligomers.  In  the  case  of 
prions,  the  oligomers  and/or  fibrils  can  fission,  which  then  re-seeds  the  initial 
conversion/aggregation  process.  Prions  are  the  only  known  amyloid  proteins  to 
spontaneously  fission  in  the  converted  form  either  in  vitro  or  in  vivo.  It  is  assumed  (but  not 
shown  explicitly  here)  that  both  oligomers  and  native  monomers  experience  some  level  of 
removal  by  cellular  systems  in  vivo. 


pH  downstream  from  the  initial  extrusion 
site.43 

3.  Chorion  in  fish  and  insect  egg  shells.  Trun¬ 
cated  peptides  from  the  central  regions  of 
two  chorion  proteins  self-assembled  into 
spherulites  possessing  (3-structure,  which 
then  converted  to  fibrils  upon  maturation, 
suggesting  amyloid  character  in  actual  egg 
shells.44 


4.  Amyloid-like  structure  in  synapses  of 
Aplysia  (slugs).  The  N-terminus  region  of 
the  Aplysia  synapse  protein  CPEB  is 
glutamine-rich,  like  huntingtin  protein. 
Engineered  expression  of  CPEB  in  yeast 
yields  prion-like  aggregates  similar  to  the 
native  ones  discussed  earlier;  hence,  prion¬ 
like  states  of  CPEB  might  effect  long-term 
strengthening  of  synaptic  contacts.45 


Amyloids  in  Materials  Science 

The  regular  diameter  and  periodicity  of 
the  amyloid  fibrils  make  them  good  mate¬ 
rials  templates.  One  group  employed  yeast 
prion-like  proteins  to  template  ~100-nm- 
wide  gold  nanowire  growth;  the  gold- 
coated  protein  filaments  after  initial  decor¬ 
ation  by  small  gold  nanoparticles  were 
found  to  bind  to  genetically  engineered 
cysteine  residues.46  Silver  nano  wires  of 
20  nm  width  were  grown  inside  filamen¬ 
tary  cross-(3  tubes  grown  from  a  dipheny¬ 
lalanine  peptide.47  Engineering  of  con- 
trollably  switched  (3-sheet  materials  could 
prove  valuable  for  tissue  growth  scaffold¬ 
ing,  as  one  example.48  Hybrid  molecules 
including  eight  amino  acids — four  per 
strand,  but  with  non-amino  bends — have 
been  developed  that  controllably  self- 
assemble  into  different  (3-sheet  structures 
(as  shown  in  Figure  5)  depending  upon 
the  pH,  analogous  to  spider  silk.49  Clearly, 
the  future  of  engineered  amyloid  struc¬ 
tures  in  materials  science  looks  bright.50 

Conclusion 

In  this  article,  we  have  developed 
themes  that  bring  together  the  fields  of 
amyloid  diseases  and  materials  science. 
Protein  misfolding  and  aggregation  phe¬ 
nomena  are  intimately  linked  to  many  se¬ 
rious  public  health  issues.  However,  many 
aspects  of  the  phenomena  have  close 
analogies  in  synthesized  materials,  and 
their  full  molecular  understanding  re¬ 
quires  experimental  and  modeling  tools 
more  familiar  in  the  physical  and  mate¬ 
rials  sciences.  In  addition  to  the  possible 
medical  breakthroughs  that  such  multi¬ 
disciplinary  studies  can  lead  to,  there  is 
a  growing  possibility  that  understand¬ 
ing  the  mechanisms  of  amyloid  formation 
can  have  wide  impact  in  fields  ranging 
from  basic  neurobiology  to  materials 
science. 
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Figure  5.  (a)-(c)  Images  from  electron  microscopy  showing  pH-dependent  fibril  growth  of  designed  peptidomimetic  molecules  (from 
Reference  49).  Scale  bars  are  100  nm. 
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Sort  Results  by: 

1.  California  researcher  finds  that  mad  cow  infection  spreads  rapidly 

Scottsbluff  Star-Herald  -  Jun  1 6  8:06  AM  ^ 

DAVIS,  Calif.  (AP)  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow 
so  frighteningly  fast. 

2.  Modeling  to  Understand  Mad  Cow  Disease 

SAP  Info -Jun  16  2:15  AM  ^ 

Two  physics  professors  at  the  University  of  California  Davis,  who  have  been  researching  mad  cow  disea 
use  computerized  mathematical  models  to  explain  how  the  disease  acts  so  fast. 

3.  Researchers  explore  spread  of  mad  cow  disease 

CTV.ca  -  Jun  15  4:49  PM  ^ 

DAVIS,  Calif.  —  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dis 
frighteningly  fast  using  an  unconventional  tool  mathematical  modeling. 

4.  UC  Researchers  Explore  Spread  of  Mad  Cow 

TechNewsWorld.com  -  Jun  15  8:27  AM  ^ 

Scientists  believe  the  human  body  usually  rejects  misfolded  proteins  without  harm;  it's  when  a  microscopi 
flawed  proteins  is  consumed  from  another  animal  that  the  disease  can  take  hold  in  a  person  or  animal  wh 
would  never  develop  the  disease. 

5.  The  Capital  City's  Newspaper  Online 
Independent  Record  -  Jun  14  11 :1 2  PM  ^ 

AP  photo  -  Rajiv  Singh,  a  physics  professor  at  the  University  of  California  Davis,  holds  up  a  slide  showint 
of  a  human  protein  cell  at  his  campus  office  in  Davis,  Calif.,  on  Tuesday. 

6.  Mad  cow  computer  modeling  finds  infection  spreads  rapidly 

Billings  Gazette  -  Jun  15  7:31  AM  ^ 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  disc 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

7.  Mad  cow  computer  modeling  finds  infection  spreads  rapidly 

Billings  Gazette  -  Jun  15  7:32  AM  ^ 

DAVIS,  Calif.  --  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  disc 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

8.  UC  Researchers  Explore  Spread  of  Mad  Cow 

First  Coast  News  -  Jun  15  10:04  AM  ^ 

DAVIS,  CA  (AP)  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  < 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

9.  Mad  cow  computer  modeling  finds  infection  spreads  rapidly 
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SanLuisObispo.com  -  Jun  15  3:35  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

10.  UC  researchers  explore  spread  of  mad  cow  ^ 

Kansas  City  Star  -  Jun  1 5  5:20  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

11.  UC  Researchers  Explore  Spread  of  Mad  Cow ^ 

RedNova  -  Jun  15  7:26  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

12.  UC  researchers  explore  spread  of  mad  cow  ^ 

The  Times  Leader  -  Jun  15  5:18  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

13.  UC  researchers  explore  spread  of  mad  cow  ^ 

Macon  Telegraph  -  Jun  15  5:28  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

14.  Mad  cow  computer  modeling  finds  infection  spreads  rapidly  ^ 

Monterey  County  Herald  -  Jun  15  3:29  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

15.  UC  researchers  explore  spread  of  mad  cow  ^ 

TimesLeader.com  -  Jun  15  5:08  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

16.  Mad  cow  computer  modeling  finds  infection  spreads  rapidly  ^ 

The  Times  Leader  -  Technology  -  Jun  15  5:48  AM 

DAVIS,  Calif.  (AP)  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow 
so  frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

17.  UC  Researchers  Explore  Spread  of  Mad  Cow  ^ 

WJLA-TV  Washington  D.C.  -  Jun  15  6:08  AM 

DAVIS,  Calif.  (AP)  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow 
(  news  )  acts  so  frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

18.  UC  researchers  explore  spread  of  mad  cow  ^ 

Kansas  City  Star  -  Jun  15  5:27  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 
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